Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,263 G→T 6.3% intergenic (+272/+135) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,2630GT 6.3% 238.5 / 4.8 91intergenic (+272/+135)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (59/25);  new base T (5/2);  total (64/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAG  >  NC_000913/1552114‑1552412
                                                                                                                                                     |                                                                                                                                                     
aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                                                                                       <  2:678651/150‑1 (MQ=255)
   tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                                                                                    <  1:1859320/150‑1 (MQ=255)
    ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                   <  1:1390066/150‑1 (MQ=255)
     cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTCATATTGAAGAAAACGAAGGCTTTGTTCATCCCCGATGCGGATTGAAAGc                                                                                                                                                  >  2:5695/1‑150 (MQ=18)
     cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGTCAATGTTGCTATTGTCATAAGTTTCAAATGGCGCTGTAAAAGGAGACAATTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACgg                                                                                                                                                  >  1:1119982/1‑149 (MQ=25)
      aaTAGATTGCAGTGAACTTGTATGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                 <  2:2087032/150‑1 (MQ=25)
      aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                 <  2:1037377/150‑1 (MQ=255)
               ctgttaaCGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGAGTTTGTCTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                        <  2:858539/145‑1 (MQ=21)
               ctgttaaCGTGTAGGCATGATATTCGTAGCTCTTCAGGAAATGTTGCTTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                        <  2:850854/145‑1 (MQ=18)
               cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCATTAAAAGGAGACATATTAATATAGTAGACAACGTAAACTTTGTTAATACCTGAAGCTCAGTGAACTCCatataagaa                                                                                                                                        >  2:1379161/1‑145 (MQ=9)
                aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                       <  2:1923693/150‑1 (MQ=255)
                  tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATATTCATATAGTAGACAACGTAGGCTTTGTTAATTCCGGATGCGGAGTGAACGCCTTATAAGGCAta                                                                                                                                     >  2:470993/1‑149 (MQ=21)
                    aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaa                                                                                                                                   <  1:1763913/150‑1 (MQ=255)
                      cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGATGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                                                                                                 >  2:2051798/1‑150 (MQ=255)
                      cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                                                                                                 >  2:8530/1‑150 (MQ=255)
                       gtgtAGGCCTGATAAGCGGAGGGCATGAGGCAAGGTTGCGTTTGTCATCAGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAATGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                                                                                                <  2:1110067/150‑1 (MQ=18)
                        tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAAAGACTTATACGGAATGAAAAcc                                                                                                                               >  1:2182851/1‑150 (MQ=25)
                        tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                                                                                               >  1:575127/1‑150 (MQ=255)
                        tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                                                                                               >  2:753192/1‑150 (MQ=255)
                           aGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGAAAAAGTCGTCATTTTCATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCGTATCCGGCATGCAAGcccgt                                                                                                                            >  2:119077/1‑148 (MQ=21)
                              ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                                                                                         <  2:1575914/150‑1 (MQ=255)
                                   tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAAAACGTAGGCTTTGTTCATTCCGGAAGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCa                                                                                                                    >  2:712139/1‑150 (MQ=255)
                                      gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                                                                                                 >  1:1009524/1‑150 (MQ=255)
                                            cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                                                                                           <  2:575127/150‑1 (MQ=255)
                                            cgaaTCAGGCAATGTTGCGTTTGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                                                                                           <  2:1416964/147‑1 (MQ=255)
                                             gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTATAAAACGTAGGCTTTGTACATGCCGGAAGCGGCGTGAACGACTTAACCGGAAAGAAAACCCTTCAAAACCAAAAGATTgaa                                                                                                          >  1:874179/1‑148 (MQ=37)
                                                  aGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGAGACATTTTAATATTGTATACAACGAAGGCTTTGTACATGCCCGATGCGGCGTGAAACACAAAAACGCCATGACAACACATCAAATACAATAGAATGCATTGaa                                                                                                     >  2:1172479/1‑150 (MQ=16)
                                                   ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCCCGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAAc                                                                                                    >  1:1613395/1‑150 (MQ=255)
                                                   ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTAAAATACAATATATTGCAGTTAAc                                                                                                    >  1:1157324/1‑150 (MQ=17)
                                                    ttaaTGTTGCGTTTTTCATCAGGTGAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                                                                                                   <  2:1203718/148‑1 (MQ=255)
                                                    gCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCCGAAGCGGCGTGAACGACTGATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                                                                                                   >  2:2069021/1‑150 (MQ=17)
                                                         gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCATAACCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                              >  2:2122554/1‑150 (MQ=255)
                                                          ttGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                                                                                             <  2:660884/150‑1 (MQ=255)
                                                           tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTGCAGTGAACGTGTAGGc                                                                                            >  2:2026068/1‑150 (MQ=255)
                                                            gCGTTTGTCATCAGTTCCAAATGGCGCTGAAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGATCATGCCGTACGCGGCTTTAACGCCTTATCCGTCATGAATACCCTTCTAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                           <  2:1119982/150‑1 (MQ=14)
                                                             cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGCTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                                                                                          >  2:665821/1‑150 (MQ=255)
                                                                  gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa                                                                                     >  2:1414906/1‑150 (MQ=255)
                                                                    catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                   >  1:300113/1‑150 (MQ=40)
                                                                    catcaGTTTCAAATGGAGCTGTAAATGGCGGCATTTTCATATTGTAGTCAACGTAGGCTTTGTGCATGCCGGATGCGTCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGTTAAGc                                                                                   <  1:119077/150‑1 (MQ=11)
                                                                          tttCAAATGGCGCTGTAAAAGGCGTCATTTTCATTTTGTAGACAACGTAGGCTTTGTTCATGCCGTATGCGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAAATGATTGCAGGAAACGTGTAGGCCTCATAACAGTAgag                                                                             >  1:406462/1‑148 (MQ=11)
                                                                           ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATACAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            >  2:1395580/1‑150 (MQ=21)
                                                                                 tGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                 <  2:1704725/75‑1 (MQ=34)
                                                                                      cTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATg                                                                 >  1:1002003/1‑150 (MQ=32)
                                                                                       tGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGt                                                                <  2:621534/150‑1 (MQ=32)
                                                                                          aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                             >  1:834085/1‑150 (MQ=32)
                                                                                          aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                             >  2:68178/1‑150 (MQ=32)
                                                                                          aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                             >  2:159814/1‑150 (MQ=32)
                                                                                           aaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCg                                                            >  1:241490/1‑150 (MQ=32)
                                                                                             aGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGtt                                                          >  1:1118524/1‑150 (MQ=32)
                                                                                             aGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTGCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGtt                                                          >  1:1357064/1‑150 (MQ=17)
                                                                                              ggCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCATTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTAGCGttt                                                         >  2:982414/1‑150 (MQ=17)
                                                                                                 gTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtc                                                      >  2:840672/1‑150 (MQ=31)
                                                                                                           tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTa                                            >  2:503588/1‑150 (MQ=32)
                                                                                                               gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                        >  1:1919404/1‑150 (MQ=34)
                                                                                                               gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                        >  1:964949/1‑150 (MQ=34)
                                                                                                                tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                       >  2:2100040/1‑150 (MQ=34)
                                                                                                                 aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCATGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgc                                      >  1:1401812/1‑150 (MQ=21)
                                                                                                                       cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTGGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                                <  2:1814392/150‑1 (MQ=37)
                                                                                                                          aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTGTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaa                             >  1:594153/1‑150 (MQ=255)
                                                                                                                            gCTTTGTTCATTCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGCCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                          >  2:1623548/1‑150 (MQ=255)
                                                                                                                            gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCAATAGATTGCAGTGAACGTGTCCGCCTGATAAGCGTAGCGCAACAGCCAATGTTGCGTATGTCATCAGTTCTAAAAGGCGCTTTATAAAGt                          >  2:1632332/1‑150 (MQ=14)
                                                                                                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAATCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGTCGATTTATAAAGt                          >  2:721065/1‑150 (MQ=255)
                                                                                                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAgc                          <  2:319465/150‑2 (MQ=255)
                                                                                                                              tttGTTCATGCCGGATGCGGCTTGAAAGCCTTATCAGGCATGTAAACCCTTCCAATCCAATAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTg                         <  1:470993/150‑1 (MQ=14)
                                                                                                                              tttGTTCATGCAGTATGCTTAGGTAACGCCTTATGCTGTTTTTAAACCCTTCTACTCCAATCGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTg                         <  1:1379161/150‑1 (MQ=14)
                                                                                                                                tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGcc                       >  1:1300729/1‑150 (MQ=255)
                                                                                                                                  ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAt                     >  1:2184328/1‑150 (MQ=255)
                                                                                                                                  ttCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCTTCAGTCAATGTTGCGTTTGACATCAGTTCTAAATGGCGCTTTATAAAGTGCCAt                     >  1:1994796/1‑150 (MQ=14)
                                                                                                                                    cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                   >  2:1036221/1‑150 (MQ=255)
                                                                                                                                    cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                   >  1:902082/1‑150 (MQ=255)
                                                                                                                                    cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttt                   <  2:1403093/150‑1 (MQ=255)
                                                                                                                                    cATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCGAttt                   >  1:21184/1‑150 (MQ=14)
                                                                                                                                     aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGTGTAGCGCATCAGGCAATGTTGCGTTGGTCATCAGTTATAAATGGCGCTTTATAAAGTGCCAtttt                  >  2:906597/1‑150 (MQ=255)
                                                                                                                                     aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                  >  1:611233/1‑150 (MQ=255)
                                                                                                                                     aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGGTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                  >  1:297890/1‑150 (MQ=14)
                                                                                                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                 <  2:1873392/150‑1 (MQ=255)
                                                                                                                                      tGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                 >  2:2061042/1‑150 (MQ=14)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                >  1:1361650/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCCGGCAATGTTGCGTTTGTCCTCAGTTCTAAATGGCGCTTTATTAAGTGCCAttttttt               >  1:1397057/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt               >  1:747620/1‑150 (MQ=255)
                                                                                                                                         cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCTCTTTATAAAGTGCCATTTTTTTa              >  2:777014/1‑150 (MQ=255)
                                                                                                                                         cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa              <  1:8530/150‑1 (MQ=255)
                                                                                                                                          ggataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAATTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt             >  1:762390/7‑150 (MQ=17)
                                                                                                                                              gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg         <  1:1654772/150‑1 (MQ=255)
                                                                                                                                               cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt        >  2:2087163/1‑150 (MQ=255)
                                                                                                                                               cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt        >  2:103237/1‑150 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATTGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa       >  2:1369783/1‑150 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttttttcgtta       >  1:1766743/1‑145 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa       >  2:1942224/1‑150 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa       >  1:1400151/1‑150 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa       <  1:457065/150‑1 (MQ=255)
                                                                                                                                                   gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAcc    >  1:330658/1‑150 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGAAGTGAACGTGTATGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATACTGAGACATTTTTATATTGATTAACCAg  >  2:360444/1‑150 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAG  >  NC_000913/1552114‑1552412

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