Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,216,660 A→C 9.9% intergenic (‑169/‑119) nfeF ← / → nfeR NADPH‑dependent ferric chelate reductase/DNA‑binding transcriptional repressor NfeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,6600AC 9.9% 311.7 / 26.6 111intergenic (‑169/‑119)nfeF/nfeRNADPH‑dependent ferric chelate reductase/DNA‑binding transcriptional repressor NfeR
Reads supporting (aligned to +/- strand):  ref base A (45/55);  new base C (7/4);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 115 total aligned reads displayed.

GTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAAC  >  NC_000913/3216511‑3216809
                                                                                                                                                     |                                                                                                                                                     
gTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGa                                                                                                                                                       <  2:617060/150‑1 (MQ=255)
 tAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGaa                                                                                                                                                      >  1:984824/1‑150 (MQ=255)
 tAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGaa                                                                                                                                                      >  1:1634853/1‑150 (MQ=255)
   aTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACg                                                                                                                                                    >  2:755571/1‑150 (MQ=255)
     cAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:1022951/150‑1 (MQ=255)
     cAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:695450/150‑1 (MQ=255)
      aGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 <  2:1003847/150‑1 (MQ=255)
       gATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                >  2:38289/1‑150 (MQ=255)
           tatCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                            >  1:387979/1‑150 (MQ=255)
            atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                           >  2:176218/1‑150 (MQ=255)
             tCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 <  2:1033608/143‑1 (MQ=255)
             tCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 >  1:1033608/1‑143 (MQ=255)
             tCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          >  1:878747/1‑150 (MQ=255)
               taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                        >  2:1413147/1‑150 (MQ=255)
               taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcg                                                                                                                                        <  1:1023421/150‑2 (MQ=255)
               taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTCCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                        <  1:328371/150‑1 (MQ=255)
                   taaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                    >  2:1497284/1‑150 (MQ=255)
                   taaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                    <  2:1473983/150‑1 (MQ=255)
                      aaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                                                                                                                 >  1:742984/1‑150 (MQ=255)
                      aaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                                                                                                                 <  1:766548/150‑1 (MQ=255)
                            cAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATg                                                                                                                           >  1:1307290/1‑150 (MQ=255)
                            cAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATg                                                                                                                           >  1:484611/1‑150 (MQ=255)
                             aaaTGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          <  1:239477/134‑1 (MQ=255)
                             aaaTGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          >  2:239477/1‑134 (MQ=255)
                               atgataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGt                                                                                                                        >  2:1448108/1‑150 (MQ=255)
                                tgataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTc                                                                                                                       <  2:878747/150‑1 (MQ=255)
                                 gataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTcc                                                                                                                      <  2:906843/150‑1 (MQ=255)
                                      tgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGtt                                                                                                                 >  1:1141106/1‑150 (MQ=255)
                                         ttGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaa                                                                                                              <  1:306870/150‑1 (MQ=255)
                                         ttGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaa                                                                                                              >  2:306870/1‑150 (MQ=255)
                                           gTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaa                                                                                                            >  1:633419/1‑150 (MQ=255)
                                           gTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaa                                                                                                            <  2:1602326/150‑1 (MQ=255)
                                            ttAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaaa                                                                                                           >  2:485050/1‑150 (MQ=255)
                                            ttAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaaa                                                                                                           >  2:479409/1‑150 (MQ=255)
                                             tAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAAt                                                                                                          <  2:1307290/150‑1 (MQ=255)
                                                    aaTAAATGCAACCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCAGTATGCGGCGTGAACGCCTTATCCGCCCTACGTGTGAGATGAGTCCAAGTTTAAAAATAAttttt                                                                                                   >  1:539568/2‑150 (MQ=255)
                                                     aTAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTa                                                                                                  <  2:1242360/150‑1 (MQ=255)
                                                      tAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTac                                                                                                 <  2:1634853/150‑1 (MQ=255)
                                                       aaaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTaca                                                                                                <  1:94087/150‑1 (MQ=255)
                                                        aaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTacac                                                                                               >  2:1121418/1‑150 (MQ=255)
                                                         aTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACt                                                                                              <  2:130786/150‑1 (MQ=255)
                                                                  aTTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAAc                                                                                     <  1:1006141/150‑1 (MQ=255)
                                                                  aTTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAAc                                                                                     <  2:307105/150‑1 (MQ=255)
                                                                     tGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCt                                                                                  <  2:228659/150‑1 (MQ=255)
                                                                      gTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTc                                                                                 <  2:1456423/150‑1 (MQ=255)
                                                                       tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                                >  1:1437734/1‑150 (MQ=255)
                                                                       tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                                <  1:1051489/150‑1 (MQ=255)
                                                                          aaCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCatat                                                                             <  2:90172/150‑1 (MQ=255)
                                                                             tGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAg                                                                          <  2:689826/150‑1 (MQ=255)
                                                                              gATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGa                                                                         <  2:1077416/150‑1 (MQ=255)
                                                                               atatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGat                                                                        <  2:984824/150‑1 (MQ=255)
                                                                               atatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGat                                                                        <  1:1664905/150‑1 (MQ=255)
                                                                               atatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGat                                                                        <  1:605370/150‑1 (MQ=255)
                                                                                tatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGata                                                                       >  2:339603/1‑150 (MQ=255)
                                                                                  tGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGAtata                                                                     <  2:912422/150‑1 (MQ=255)
                                                                                   gTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATaa                                                                    >  1:1201701/1‑150 (MQ=255)
                                                                                    tCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATaaa                                                                   <  1:24292/150‑1 (MQ=255)
                                                                                     cTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAAt                                                                  >  1:464390/1‑150 (MQ=255)
                                                                                       ataGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTa                                                                >  2:642039/1‑150 (MQ=255)
                                                                                        taGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAg                                                               <  1:176218/150‑1 (MQ=255)
                                                                                                  aGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTaa                                                     >  1:1465392/1‑150 (MQ=255)
                                                                                                       aaCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAg                                                >  1:855099/1‑150 (MQ=255)
                                                                                                          cGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAAtttatatc                                             <  2:633419/150‑7 (MQ=255)
                                                                                                          cGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaa                                             >  2:125307/1‑150 (MQ=255)
                                                                                                           gCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaaa                                            >  1:502973/1‑150 (MQ=255)
                                                                                                           gCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaaa                                            >  1:814584/1‑150 (MQ=255)
                                                                                                             ctctTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAag                                          <  2:1141106/150‑1 (MQ=255)
                                                                                                              tcttctCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGACTCCACTTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAagg                                         <  2:60783/150‑1 (MQ=255)
                                                                                                                ttctCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAaggag                                       <  1:180120/150‑1 (MQ=255)
                                                                                                                  ctcAGGCATCATTACTCAACGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGc                                     >  2:519950/1‑150 (MQ=255)
                                                                                                                  ctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGc                                     >  1:781659/1‑150 (MQ=255)
                                                                                                                    cAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTg                                   <  1:1413147/150‑1 (MQ=255)
                                                                                                                    cAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTg                                   >  1:431450/1‑150 (MQ=255)
                                                                                                                     aGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGa                                  <  2:949139/150‑1 (MQ=255)
                                                                                                                     aGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGa                                  >  2:774265/1‑150 (MQ=255)
                                                                                                                      ggCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGat                                 >  1:262569/1‑150 (MQ=255)
                                                                                                                      ggCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGat                                 >  1:994079/1‑150 (MQ=255)
                                                                                                                        catcatTACTCAATGCCGGATGCGTCGTGAACTCCTTATCCGGCCTACGTTTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGatat                               <  2:539568/150‑1 (MQ=255)
                                                                                                                         atcatTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATg                              >  1:87396/1‑150 (MQ=255)
                                                                                                                          tcatTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGa                             <  2:742984/150‑1 (MQ=255)
                                                                                                                           catTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAg                            <  1:38289/150‑1 (MQ=255)
                                                                                                                            atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGc                           >  1:1041370/1‑150 (MQ=255)
                                                                                                                               aCTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcat                        >  2:362499/1‑150 (MQ=255)
                                                                                                                                cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc                       <  1:86319/150‑1 (MQ=255)
                                                                                                                                acgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc                       >  1:113013/7‑150 (MQ=38)
                                                                                                                                  cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatcat                     >  1:1379452/1‑150 (MQ=255)
                                                                                                                                   aaCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCatcatc                    >  2:1287183/1‑150 (MQ=255)
                                                                                                                                    aCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCAtcatca                   <  2:781659/150‑1 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCAc                  >  2:264474/2‑150 (MQ=255)
                                                                                                                                     cGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCAc                  <  2:153172/150‑1 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGa                >  2:230156/1‑150 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGa                >  1:1518779/1‑150 (MQ=255)
                                                                                                                                         ggataaggcgTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatca                      >  1:1031467/7‑142 (MQ=255)
                                                                                                                                         ggataaggcgTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatca                      <  2:1031467/136‑1 (MQ=255)
                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGgtgt         >  2:751805/1‑150 (MQ=255)
                                                                                                                                               ggcgTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgtt        <  1:264474/150‑1 (MQ=255)
                                                                                                                                                gcgTGAACTCCTTATCCGGCATACTTGTGAGATTAGTCCACTTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgttg       <  1:285500/150‑1 (MQ=255)
                                                                                                                                                  gTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTa     <  2:387979/150‑1 (MQ=255)
                                                                                                                                                   tGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTaa    <  2:855099/150‑1 (MQ=255)
                                                                                                                                                     aaCGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAAc  <  2:746307/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAAC  >  NC_000913/3216511‑3216809

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