Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,465,424 T→G 22.9% A41A (GCA→GCC fcuA_1 ← Ferrichrome receptor FcuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,465,4240TG22.9% 94.5 / 7.0 58A41A (GCA→GCCfcuA_1Ferrichrome receptor FcuA
Reads supporting (aligned to +/- strand):  ref base T (9/32);  new base G (12/3);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTC  >  ACB122/3465294‑3465555
                                                                                                                                      |                                                                                                                                    
cAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCGCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCTAAGCGCGCCAGTACATCGGCCAACGGGCCTGCgg                                                                                                                                  <  7:185937/135‑1 (MQ=255)
cAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCgg                                                                                                                                  <  4:199431/135‑1 (MQ=255)
       ttGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAggtg                                                                                                                           <  8:456201/135‑1 (MQ=255)
                  tGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGGCAAGCGCGCCCGTTCATCGGGCCCACGGCCGGCGGCCAGGGGGGAAGTgcggg                                                                                                                >  1:70004/1‑133 (MQ=255)
                  tGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGGCAAGCGCGCCCGGTCATCGGGCCACGGGCCGGCGGCCCGGGGGGTAGTgggcg                                                                                                                >  6:181454/1‑135 (MQ=255)
                  tGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGGCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTgcgc                                                                                                                 <  7:193840/134‑1 (MQ=255)
                          cccTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGATAACGCCACCCCGGCGTTGGGGGCCAACCGCGCCCGTTCATAGGCCCACCGGCCGGCGGCCCGGGGGGTAGTGCGCGCGGGGgtg                                                                                                        >  2:30438/1‑135 (MQ=255)
                          cccTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGGCAACCGCGCCCGGACCTCGGCCCCCCGGGCGGGCGCCCCGGGGGTAGTGCGCGCGGTgggg                                                                                                        >  8:302043/1‑133 (MQ=255)
                            cgccAGGCCTTGTTTTTGTACGCCCCTAAGAAGGGCGGGGTCTAATTCCAACGCCACCCAGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTc                                                                                                      <  2:405904/133‑1 (MQ=255)
                            cTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGCCAAACGCGCCACGACCTCGGGCAACCGGCCGGGGGGCCCGTGGTTAGTGGGGGCGGAgggggg                                                                                                      >  4:430790/1‑131 (MQ=255)
                            cTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCCGCCAAGCGCGCCCCTAAATAGGCCCACCGGCCTGCGGCCACGTGGGAAATGGCCGCGGTGgtgtg                                                                                                      >  8:335143/1‑134 (MQ=255)
                               cAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTGGGCGGTAAAACGCCCCCCTACAACGGGCAACCGGCCGGCGGCCAGGTGGGTAGATCGCGCGGAGGTGTCCGc                                                                                                   >  1:449857/1‑135 (MQ=255)
                                          tttGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGCAAAGCCCGCCAGTACCTAGGGCCACCGGGCGGGGGGCCGGGGGTTAGTGCGCGGGGGGGGGTGGGCTTCCCTCgggg                                                                                       >  4:138904/1‑135 (MQ=255)
                                                   cccTGAGCCGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCa                                                                              <  8:413353/135‑1 (MQ=255)
                                                     cTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGCGGCAACCGCCCCCCGAAAAAGGCCCACCCGGCCGGCGGCACGTGGGTAGTTCGCGGGGCGGTGTCGGCAGCCCACGCCCTggggggccgg                                                                            >  4:45638/1‑130 (MQ=255)
                                                     cTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGGCAACCCCCCCCGTACATCGGGCCCCCGGCCGGCGGCCGCGTGGTTATTGCGGCCGGGGGTGTGGGCTGCCCTCGCGGGGGGGGGCAgg                                                                            >  2:122979/1‑135 (MQ=255)
                                                      tGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGCCAACCGCGCCCGTACAACGGGCCACGGGGCAGGGGCCCCGGGGGTAGTGGGCGCGGGGGGGTCGGCTGGCCTAgggggggggggcccgg                                                                           >  1:68364/1‑129 (MQ=255)
                                                      tGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGCAAAAAGCGCCATTACATACGCCAAACGGCCTGCGGCCACGTGGGTACTGCGCGCCGCGGTGTCGGCAGCCAAAGCCGGGGGGGGCAggg                                                                           >  2:296593/1‑135 (MQ=255)
                                                      tGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGGCAACCGCGCCCGGTAATCGGGCAACGGGCCAGCGGGCCGGTGGGTAATGCGCGCGGTGGGGTGGGCTGCCCTAGCGGGGGGGGGCCggg                                                                           >  4:94964/1‑135 (MQ=255)
                                                         gctgggCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGggcg                                                                         <  3:122584/131‑1 (MQ=255)
                                                         gCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGggcgg                                                                        <  6:391184/135‑1 (MQ=255)
                                                             ggCCGGGGCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCAt                                                                    <  4:280768/135‑1 (MQ=255)
                                                              gCCGGGCCGAATCACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATg                                                                   <  5:51698/135‑1 (MQ=255)
                                                                cGGGTCGAATGACAACGCCCCCCCGGCGTTGGGGGCACAGCGCGCCAGGAAATAGGGCCACGGGGCGGGGGGCCGGGGGGTAGTGGGCGCGGTGGTGGGGGGGGGCCTCGCGGGGGGGGGCCGGGCGGGCCTggg                                                                 >  8:448846/1‑134 (MQ=255)
                                                                cGGGTCGAATGACAACGCCACCCCGGCGTTGGGGGCCAACCCCGCCCGTTCAACGGGGCACGGGCCGGCGGCCCGGGGGGTATTGGGCGCGGGGGGGTGGGGTGCCCTAGCGGGGGGGGGCCGGGGGGGCCAggg                                                                  >  4:206229/1‑135 (MQ=255)
                                                                cGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCCCGACATCGGCCCACCGGCCTGCGGGCCGGGGGTTAGTGGGCGGGGAGGGGTCGGGTGCCCTCGCGGGGGGGGGCCGGGCGGcccaggg                                                                 >  6:205460/1‑130 (MQ=255)
                                                                   gccGAATCACAACGCCCCCCCGGCGTGGGCGGCAACGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  1:168668/133‑1 (MQ=255)
                                                                   gTCGAATGACAACGCCACCCCGGCTTTGGCGGCAACGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCTGGTGGTAAGTGCGCGCTGAGGTGTCGGCTGCAATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  7:56682/135‑1 (MQ=255)
                                                                   gTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  4:180628/135‑1 (MQ=255)
                                                                   gTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  4:77879/135‑1 (MQ=255)
                                                                   gTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  7:63462/135‑1 (MQ=255)
                                                                   gTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCCGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                              <  4:252718/135‑1 (MQ=255)
                                                                    tCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                             <  7:137553/135‑1 (MQ=255)
                                                                    tCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                             <  4:85042/135‑1 (MQ=255)
                                                                    tCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                             <  4:389780/135‑1 (MQ=255)
                                                                    tCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGGACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                             <  6:43722/135‑1 (MQ=255)
                                                                     cGAATGACAACGCCACCCCGGCGTTGGGGGCCCAGCGCGGCAGTGCAACGGGCCACCGGGCGGCGGGCAGGGGGGTAGTGGGGGCGGCGGTGTCGGGGGCCCTCGCGGGGGGGGGCAGGGCGGGCCAGGCGAggg                                                            >  6:441475/1‑135 (MQ=255)
                                                                                          cGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTTGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAg                                       <  3:308645/135‑1 (MQ=255)
                                                                                          cGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAg                                       <  3:428547/135‑1 (MQ=255)
                                                                                          cGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAg                                       <  4:297219/135‑1 (MQ=255)
                                                                                             tGGCGGCAAAGCGCGCCTGTACATCGGCCAACGGGCCGGCGGCCAGGTTGTATGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGACATGGCGATGGCTAGGCTTAGGGCCTTCatatgct                                    <  4:182420/135‑7 (MQ=255)
                                                                                                          gcgcCAGTACATCGGCCAACGGGCCTGCGGCCAGGTTGTAAGTGCGCGCGGAGGTGGCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCt                        <  8:323066/135‑1 (MQ=255)
                                                                                                          gcgcCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCTGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTATGGCCTTCAAAGGCAAGCGGTACGGCt                        <  5:183364/135‑1 (MQ=255)
                                                                                                            gcCAGTCCATCGGCCAACGGGCCGGCGGCCAGGTGGTAAGTGCGCGCGGAGGGGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTgg                      <  3:288687/135‑1 (MQ=255)
                                                                                                            gcCAGTACATCGGCCAACGGGCCTGCGGCCAGGGGGTAAGTGCGCGCGGAGGTGGCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTgg                      <  1:304093/135‑1 (MQ=255)
                                                                                                             ccAGTACCTCGGCCAAGGGGCGGGGGGGCTGGTGGTACGGGCGCGCGGAGGGGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGc                     <  5:18338/135‑1 (MQ=255)
                                                                                                              cAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTgg                      >  6:146096/1‑133 (MQ=255)
                                                                                                                  aCAGCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGGGGCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGc                <  1:429764/135‑1 (MQ=255)
                                                                                                                    aTCGTCCAACGGGCCTGCGGCCAGGTGGTAAGTTCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCg               <  7:192101/134‑1 (MQ=255)
                                                                                                                         ggCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCAt           <  8:202402/135‑1 (MQ=255)
                                                                                                                          gCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCAt           <  5:148347/134‑1 (MQ=255)
                                                                                                                          gCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATg          <  1:304758/135‑1 (MQ=255)
                                                                                                                              aaCGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGt       >  5:256157/1‑135 (MQ=255)
                                                                                                                                 gggCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGGGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTGAAGGGTggg    >  1:456958/1‑135 (MQ=255)
                                                                                                                                 gggCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTggg    >  5:409431/1‑135 (MQ=255)
                                                                                                                                  ggCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGt   <  5:443327/135‑1 (MQ=255)
                                                                                                                                   gCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTc  >  4:457224/1‑135 (MQ=255)
                                                                                                                                   gCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGACCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTc  >  6:154473/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGCGGCAAAGCGCGCCAGTACATCGGCCAACGGGCCTGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTC  >  ACB122/3465294‑3465555

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: