Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,216,659 G→T 5.0% intergenic (‑168/‑120) 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,6590GT 5.0% 404.9 / 6.9 119intergenic (‑168/‑120)nfeF/nfeRNADPH‑dependent ferric chelate reductase/DNA‑binding transcriptional repressor NfeR
Reads supporting (aligned to +/- strand):  ref base G (49/64);  new base T (2/4);  total (51/68)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 122 total aligned reads displayed.

GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAA  >  NC_000913/3216510‑3216808
                                                                                                                                                     |                                                                                                                                                     
ggTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACACAACGCCGCAAGCGGCGTg                                                                                                                                                       >  2:1534871/1‑150 (MQ=255)
 gTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGa                                                                                                                                                      >  1:1184917/1‑150 (MQ=255)
  tAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGaa                                                                                                                                                     >  2:38358/1‑150 (MQ=255)
   aaTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAAc                                                                                                                                                    >  2:1539994/1‑150 (MQ=255)
   aaTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAAc                                                                                                                                                    >  2:601879/1‑150 (MQ=255)
      cAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGcc                                                                                                                                                 >  1:474487/1‑150 (MQ=255)
         atatatCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGccgca                                                                                                                                              <  1:66747/150‑4 (MQ=255)
            tatCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                           >  2:1408087/1‑150 (MQ=255)
             atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACg                                                                                                                                                   >  2:1346596/1‑141 (MQ=255)
             atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACg                                                                                                                                                   <  1:1346596/141‑1 (MQ=255)
             atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACg                                                                                                                                                   <  1:1352197/141‑1 (MQ=255)
             atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACg                                                                                                                                                   >  2:1352197/1‑141 (MQ=255)
             atCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          >  2:1458738/1‑150 (MQ=255)
             atCTAAATAAAACTCCCAAATCAAAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATGAGACGAACCGCCTCTTCTCAGGCATCATTACTCAAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  1:407239/150‑1 (MQ=255)
                taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       >  1:1117043/1‑150 (MQ=255)
                taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcg                                                                                                                                       <  1:1044782/150‑2 (MQ=255)
                 aaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgt                                                                                                                                      <  2:903905/150‑3 (MQ=255)
                 aaataaaACTAGCAAATGATAATCATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGACTCTTCTCAGGCATCATTACTAAACGACGGATGCGGCGTGAACGCCTTATCCgagat                                                                                                                                      >  2:407239/1‑146 (MQ=255)
                   ataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgttg                                                                                                                                    <  2:1136032/150‑5 (MQ=255)
                    taaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                   <  2:480640/150‑1 (MQ=255)
                     aaaaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgt                                                                                                                                  >  1:1182431/1‑150 (MQ=255)
                      aaaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtg                                                                                                                                 >  2:389114/1‑150 (MQ=255)
                       aaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                   >  1:735884/1‑147 (MQ=255)
                       aaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                   <  2:735884/147‑1 (MQ=255)
                        aCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                                                                                                                               <  2:474487/150‑1 (MQ=255)
                        aCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                                                                                                                               <  1:92017/150‑1 (MQ=255)
                         cTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTga                                                                                                                              >  2:497309/1‑150 (MQ=255)
                          tCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgag                                                                                                                             <  1:601879/150‑1 (MQ=255)
                          tCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgag                                                                                                                             <  1:201199/150‑1 (MQ=255)
                           cGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgaga                                                                                                                            >  2:791001/1‑150 (MQ=255)
                             cAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATg                                                                                                                          >  1:1387524/1‑150 (MQ=255)
                                      atgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGt                                                                                                                 >  1:1426884/1‑150 (MQ=255)
                                        gatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                                                                                                               <  2:1127652/150‑1 (MQ=255)
                                         atTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTa                                                                                                              >  2:266413/1‑150 (MQ=255)
                                          ttGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaa                                                                                                             >  2:1356325/1‑150 (MQ=255)
                                            gTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaa                                                                                                           <  1:680470/150‑1 (MQ=255)
                                            aacaaTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaa                                                                                                           <  2:552421/147‑1 (MQ=255)
                                                    tGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAtttt                                                                                                   >  1:1417192/1‑150 (MQ=255)
                                                         aaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTacac                                                                                              <  1:336900/150‑1 (MQ=255)
                                                         aaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTacac                                                                                              <  1:342831/150‑1 (MQ=255)
                                                         aaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTacac                                                                                              >  1:1496429/1‑150 (MQ=255)
                                                         aaTGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGGCCACGTTTAAAAATCATTTTTacac                                                                                              >  2:1588538/1‑150 (MQ=255)
                                                            gCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTg                                                                                           <  1:319487/150‑1 (MQ=255)
                                                             cAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGc                                                                                          <  2:1230841/150‑1 (MQ=255)
                                                              aaGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCa                                                                                         <  1:1083434/150‑1 (MQ=255)
                                                                gCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAg                                                                                       <  1:497309/150‑1 (MQ=255)
                                                                    tttGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACg                                                                                   >  2:37645/1‑150 (MQ=255)
                                                                      tGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCt                                                                                 >  1:911650/1‑150 (MQ=255)
                                                                       gTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTc                                                                                <  1:1633209/150‑1 (MQ=255)
                                                                        tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                               <  1:791001/150‑1 (MQ=255)
                                                                        tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                               <  2:186824/150‑1 (MQ=255)
                                                                         aGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCat                                                                              >  2:1151288/1‑150 (MQ=255)
                                                                         aGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTAAACGCCGGATGCCGAGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCat                                                                              >  2:490839/1‑150 (MQ=255)
                                                                            aCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATc                                                                           >  2:846084/1‑150 (MQ=255)
                                                                             cTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCa                                                                          >  2:775979/1‑150 (MQ=255)
                                                                                  atGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGatat                                                                     >  1:1042047/1‑150 (MQ=255)
                                                                                       tataGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAAtt                                                                >  1:1554261/1‑150 (MQ=255)
                                                                                           gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                            <  1:1458738/150‑1 (MQ=255)
                                                                                             cTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGatat                                                          >  1:1264158/1‑150 (MQ=255)
                                                                                                  aaGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCTGCCTTCGTTTGAGATGAGTCCTCGTTTAAAATTCATTATTACACTTTCTAGAACGCTCATATCAGATATAAATTAGAtatatata                                                     <  2:1293314/150‑1 (MQ=255)
                                                                                                    gACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAATATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAAt                                                   <  2:1426884/150‑1 (MQ=255)
                                                                                                     aCGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCATGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAAtt                                                  <  1:1366407/150‑1 (MQ=255)
                                                                                                     aCGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAAtt                                                  <  1:1479413/150‑1 (MQ=255)
                                                                                                      cGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTa                                                 >  1:1070064/1‑150 (MQ=255)
                                                                                                      cGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTa                                                 >  1:1251004/1‑150 (MQ=255)
                                                                                                      cGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTa                                                 <  1:795155/150‑1 (MQ=255)
                                                                                                         aCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGc                                              >  1:349780/1‑150 (MQ=255)
                                                                                                         aCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGc                                              <  2:204907/150‑1 (MQ=255)
                                                                                                          ccGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCa                                             <  2:783965/150‑1 (MQ=255)
                                                                                                           cGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaa                                            >  2:1414191/1‑150 (MQ=255)
                                                                                                               tcttctCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAagg                                        <  2:301770/150‑1 (MQ=255)
                                                                                                               tcttctCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAagg                                        <  2:1184917/150‑1 (MQ=255)
                                                                                                                  tctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAaggagg                                     <  2:1251004/150‑1 (MQ=255)
                                                                                                                   ctcAGGCATCATTCCTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCCCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGc                                    <  2:1203282/150‑1 (MQ=255)
                                                                                                                   ctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGc                                    <  2:633539/150‑1 (MQ=255)
                                                                                                                     cAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTg                                  <  2:1264158/150‑1 (MQ=255)
                                                                                                                       ggCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGACGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGat                                <  2:1627276/150‑1 (MQ=255)
                                                                                                                         catcatTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTCTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGatat                              >  1:1229204/1‑150 (MQ=255)
                                                                                                                            catTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                       >  2:922524/1‑122 (MQ=255)
                                                                                                                            catTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                       <  1:922524/122‑1 (MQ=255)
                                                                                                                             atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAAAAAAAGGAGGCTGatataaac                          >  2:693259/1‑146 (MQ=255)
                                                                                                                              ttACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGcc                         <  1:520616/150‑1 (MQ=255)
                                                                                                                               tACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCca                        <  2:240387/150‑1 (MQ=255)
                                                                                                                                 cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc                      >  2:63035/1‑150 (MQ=255)
                                                                                                                                  tCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatca                     <  2:311885/150‑1 (MQ=255)
                                                                                                                                  tCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatca                     <  2:1512145/150‑1 (MQ=255)
                                                                                                                                   cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGata                               >  2:494438/1‑139 (MQ=255)
                                                                                                                                   cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGata                               <  1:494438/139‑1 (MQ=255)
                                                                                                                                   cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatcat                    <  1:409921/150‑1 (MQ=255)
                                                                                                                                    aaCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCatcatc                   >  2:151749/1‑150 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACg                <  2:205459/150‑1 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACg                <  1:65792/150‑1 (MQ=255)
                                                                                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACg                >  1:581693/1‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCCTCATCCCgg               >  2:1362715/1‑149 (MQ=255)
                                                                                                                                          ggATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAg             >  2:399812/1‑150 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgtt       <  1:266413/150‑1 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgtt       <  1:1539994/150‑1 (MQ=255)
                                                                                                                                                 gcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgttg      <  1:427267/150‑1 (MQ=255)
                                                                                                                                                 gcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgttg      <  2:1117043/150‑1 (MQ=255)
                                                                                                                                                   gTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTa    <  2:843102/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAA  >  NC_000913/3216510‑3216808

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