Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 551,488 A→G 5.2% intergenic (+156/+39) ybcF → / ← purK putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,4880AG 5.2% 505.6 / 11.1 151intergenic (+156/+39)ybcF/purKputative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base A (77/66);  new base G (2/6);  total (79/72)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCG  >  NC_000913/551339‑551632
                                                                                                                                                     |                                                                                                                                                
gCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATa                                                                                                                                                  >  1:1323131/1‑150 (MQ=255)
 cATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATaa                                                                                                                                                 <  2:1252864/150‑1 (MQ=255)
   ttAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAgg                                                                                                                                               >  1:1011383/1‑150 (MQ=255)
   ttAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAgg                                                                                                                                               >  1:1011366/1‑150 (MQ=255)
     aGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCg                                                                                                                                             >  1:1552803/1‑150 (MQ=255)
      gACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGt                                                                                                                                            <  2:56597/150‑1 (MQ=255)
      gACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGt                                                                                                                                            <  2:489718/150‑1 (MQ=255)
       tcATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGtt                                                                                                                                           <  1:1329312/149‑1 (MQ=255)
       aCATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGtt                                                                                                                                           <  1:930531/150‑1 (MQ=255)
           ttGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                       >  1:390188/1‑150 (MQ=255)
            tgtgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                      >  2:142979/1‑150 (MQ=255)
             gtgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                     <  1:1179674/150‑1 (MQ=255)
              tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                    <  1:1140730/150‑1 (MQ=255)
              tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                    >  2:1236701/1‑150 (MQ=255)
                ccTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGCACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                  >  2:1293582/1‑150 (MQ=255)
                 cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  1:205093/150‑5 (MQ=37)
                 cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  1:1450831/150‑5 (MQ=37)
                 cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  2:850219/150‑5 (MQ=37)
                    aTGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                              <  2:29216/150‑1 (MQ=255)
                    aTGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                              <  2:605959/150‑1 (MQ=255)
                     tGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                             <  1:1217347/150‑1 (MQ=255)
                     tGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                             <  1:408311/150‑1 (MQ=255)
                        gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCac                                                                                                                          >  2:1598937/1‑149 (MQ=255)
                        gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                                                                                                          >  2:1598953/1‑150 (MQ=255)
                        gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                                                                                                          >  2:117121/1‑150 (MQ=255)
                         aCGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAg                                                                                                                         >  2:1431701/1‑150 (MQ=255)
                          cGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGa                                                                                                                        <  1:1458041/150‑1 (MQ=255)
                           gCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGaa                                                                                                                       >  2:1115535/1‑150 (MQ=255)
                                aCGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACGAATagag                                                                                                                  >  1:1495833/1‑150 (MQ=255)
                                aCGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACGAATagag                                                                                                                  >  2:258864/1‑150 (MQ=255)
                                aCGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATagag                                                                                                                  <  2:1114053/150‑1 (MQ=255)
                                aCGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATagag                                                                                                                  >  1:638891/1‑150 (MQ=255)
                                  gcgTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCa                                                                                                                <  1:99545/150‑1 (MQ=255)
                                    gTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCAcc                                                                                                              <  2:444095/150‑1 (MQ=255)
                                            aGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACGAATAGAGCACCCGTTAAcc                                                                                                      >  2:513684/1‑150 (MQ=255)
                                            aGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAAcc                                                                                                      >  2:874024/1‑150 (MQ=255)
                                            aGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAAcc                                                                                                      >  2:306912/1‑150 (MQ=255)
                                               ccTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                   <  2:1350723/150‑1 (MQ=255)
                                               ccTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                   <  2:43246/150‑1 (MQ=255)
                                               ccTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGATTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                   <  2:1011383/150‑1 (MQ=255)
                                                cTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                  >  1:530298/1‑150 (MQ=255)
                                                cTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                  >  2:898411/1‑150 (MQ=255)
                                                cTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                  >  2:959515/1‑150 (MQ=255)
                                                   cAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTa                                                                                               <  2:1124488/150‑1 (MQ=255)
                                                   cAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTa                                                                                               <  1:172519/150‑1 (MQ=255)
                                                      ccGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTActc                                                                                            <  1:909222/150‑1 (MQ=255)
                                                           gCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTgcgc                                                                                       <  1:1127290/150‑1 (MQ=255)
                                                           gCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTgcgc                                                                                       <  1:405432/150‑1 (MQ=255)
                                                           gCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTgcgc                                                                                       <  1:405601/150‑1 (MQ=255)
                                                                aTCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCaaa                                                                                  >  2:1714917/1‑150 (MQ=255)
                                                                  ccGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATc                                                                                <  1:1318427/150‑1 (MQ=255)
                                                                    ggCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCAc                                                                              <  2:845723/150‑1 (MQ=255)
                                                                            gTGCACAATGCCTGTTGCGATGCTTGACGCATCTTATCAGGCCTACATTGGGTACCGGTTGGGTTGGCCGGATAAGGCGTTTACGCCGCATCCGGAAAGAATAGAGCACCTGTTAACCTAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                      <  1:520780/150‑1 (MQ=37)
                                                                            gTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                      >  2:163444/1‑150 (MQ=255)
                                                                            gTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                      >  2:163462/1‑150 (MQ=255)
                                                                            gTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                      >  2:357624/1‑150 (MQ=255)
                                                                            gTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                      >  2:357666/1‑150 (MQ=255)
                                                                             tGCACAAGGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGc                                                                     <  1:1641348/150‑1 (MQ=255)
                                                                                          atgcgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCggcggc                                                        >  2:1195484/1‑150 (MQ=255)
                                                                                          atgcgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCggcggc                                                        >  2:1195490/1‑150 (MQ=255)
                                                                                             cgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAATTCACGCCGCTGGCATATTCCGGCGgcaac                                                     >  2:611259/1‑148 (MQ=255)
                                                                                             cgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagc                                                     >  1:684052/1‑150 (MQ=255)
                                                                                              gatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagca                                                    >  2:1169802/1‑150 (MQ=255)
                                                                                                  cTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCgg                                                <  1:1468820/150‑1 (MQ=255)
                                                                                                   ttGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCggg                                               >  2:416876/1‑150 (MQ=255)
                                                                                                   ttGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCggg                                               >  2:944408/1‑150 (MQ=255)
                                                                                                        gCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTaa                                          <  1:1576417/150‑1 (MQ=255)
                                                                                                         cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGTGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                         <  2:1280006/150‑1 (MQ=255)
                                                                                                         cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                         <  2:244974/150‑1 (MQ=255)
                                                                                                         cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                         <  1:1151062/150‑1 (MQ=255)
                                                                                                              tATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTc                                    >  2:197230/1‑150 (MQ=255)
                                                                                                              tATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTc                                    <  2:1036761/150‑1 (MQ=255)
                                                                                                                tCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCa                                  <  1:74728/150‑1 (MQ=255)
                                                                                                                tCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCa                                  <  1:574667/150‑1 (MQ=255)
                                                                                                                      cTACACTGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGACTAGAGCACCAGTTAACCGAACTTACTGTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCg                            <  2:308336/150‑1 (MQ=255)
                                                                                                                      cTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCg                            <  1:659337/150‑1 (MQ=255)
                                                                                                                      cTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCg                            >  1:630485/1‑150 (MQ=255)
                                                                                                                       tACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCAACCAGCGTCGc                           >  1:1087614/1‑150 (MQ=255)
                                                                                                                            tGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCa                      <  1:1188178/150‑1 (MQ=255)
                                                                                                                             gggTCCCGGATCGGTAGGCCGGATAAGGCGTTTACGACGCATCCGGCAAGAATAGAGCACCAGGTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGAGGGATTAAGGCTTCCAGCGTCGCAGTCAg                     >  2:1511534/1‑150 (MQ=255)
                                                                                                                             gggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAg                     >  2:1511548/1‑150 (MQ=255)
                                                                                                                              ggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGa                    >  1:724773/1‑150 (MQ=255)
                                                                                                                              ggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGa                    >  1:1356025/1‑150 (MQ=255)
                                                                                                                              ggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGa                    >  1:1355586/1‑150 (MQ=255)
                                                                                                                               gTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATCTTCCGGCGGCAGCAGCGGGAGTAAGGCTTCCAGCGTCGCAGTCAGAc                   >  1:877977/1‑150 (MQ=255)
                                                                                                                               gTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGAc                   >  1:1109205/1‑150 (MQ=255)
                                                                                                                                   cGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGa               <  1:959515/150‑1 (MQ=255)
                                                                                                                                   cGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGa               <  2:390188/150‑1 (MQ=255)
                                                                                                                                     gATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAtg             >  1:1156108/1‑150 (MQ=255)
                                                                                                                                     gATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAtg             >  2:1720351/1‑150 (MQ=255)
                                                                                                                                     gATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGAAGCGGGATTAAGGCTTCAAGCGTCGCTGTCAGACGCGAtg             >  2:1161542/1‑150 (MQ=255)
                                                                                                                                       tCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAtgtg           <  1:1195490/150‑1 (MQ=255)
                                                                                                                                          ggaGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTTCTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCg        <  1:1598937/148‑1 (MQ=255)
                                                                                                                                          gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCg        >  2:202681/1‑150 (MQ=255)
                                                                                                                                          gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCg        <  1:1598953/150‑1 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGt    >  2:702931/1‑150 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGt    >  2:703072/1‑150 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTAAACGCCTTATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTACGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTc   >  1:779444/1‑150 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTc   >  2:943454/1‑150 (MQ=255)
                                                                                                                                                cGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCg  <  2:1156108/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCG  >  NC_000913/551339‑551632

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