Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 551,496 T→A 5.0% intergenic (+164/+31) 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,4960TA 5.0% 352.7 / 6.7 99intergenic (+164/+31)ybcF/purKputative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base T (49/45);  new base A (3/2);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 106 total aligned reads displayed.

ATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGA  >  NC_000913/551348‑551644
                                                                                                                                                    |                                                                                                                                                    
aTTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGttt                                                                                                                                                      <  3:354330/149‑1 (MQ=255)
aTTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGttt                                                                                                                                                      >  3:290470/1‑149 (MQ=255)
    gtgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                  <  1:486511/149‑1 (MQ=255)
     tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAATCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 >  1:11769/1‑149 (MQ=255)
       ccTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACAGGATCGGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                               >  4:610911/1‑149 (MQ=255)
        ctgataagACGCGTCAAGCGTCGCATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                              >  1:746011/8‑149 (MQ=25)
        cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgcctta                                                                                                                                              <  3:193926/149‑4 (MQ=37)
               gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       >  4:78588/1‑149 (MQ=255)
               gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  1:701393/149‑1 (MQ=255)
               gACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       >  2:21454/1‑149 (MQ=255)
                     tGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATa                                                                                                                                 <  1:251313/149‑1 (MQ=255)
                         gcgTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGc                                                                                                                             <  2:839194/149‑1 (MQ=255)
                          cgTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACTATGCGTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCa                                                                                                                            <  4:577812/149‑1 (MQ=255)
                            tCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                             >  1:259189/1‑126 (MQ=25)
                            tCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                             <  2:259189/130‑5 (MQ=25)
                             cTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCa                                                                                                                         <  4:378317/149‑1 (MQ=255)
                             cTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCa                                                                                                                         <  2:452977/149‑1 (MQ=255)
                                 tCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTa                                                                                                                     >  2:173027/1‑149 (MQ=255)
                                   aGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAAc                                                                                                                   >  3:107091/1‑149 (MQ=255)
                                     gCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCg                                                                                                                 <  4:107091/149‑1 (MQ=255)
                                       cTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                               >  2:315265/1‑149 (MQ=255)
                                       cTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                               >  1:923386/1‑149 (MQ=255)
                                        tACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                              <  1:36151/149‑1 (MQ=255)
                                        tACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                              <  1:342271/149‑1 (MQ=255)
                                        tACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAAc                                                                                                              >  1:777148/1‑149 (MQ=255)
                                         aCAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACt                                                                                                             >  1:725786/1‑149 (MQ=255)
                                          cAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 >  1:786075/1‑112 (MQ=255)
                                          cAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 <  2:786075/112‑1 (MQ=255)
                                             ccGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTAct                                                                                                         >  1:674796/1‑149 (MQ=255)
                                              aattGCCGGATGCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTActc                                                                                                        >  2:318393/4‑149 (MQ=255)
                                                           ggCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                               >  1:560927/1‑129 (MQ=255)
                                                           ggCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGaa                                                                                                               <  2:560927/129‑1 (MQ=255)
                                                           ggCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCa                                                                                           >  3:863227/1‑149 (MQ=255)
                                                              aaTTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                                                                                                                      >  3:212449/1‑103 (MQ=255)
                                                              aaTTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                                                                                                                      <  4:212449/103‑1 (MQ=255)
                                                              aaTTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCAcgc                                                                                        >  4:27023/1‑149 (MQ=255)
                                                                   gTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTg                                                                                   <  1:173027/149‑1 (MQ=255)
                                                                    tGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                                  >  3:150667/1‑149 (MQ=255)
                                                                    tGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTgg                                                                                  >  3:493730/1‑149 (MQ=255)
                                                                     gCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGc                                                                                 >  4:689559/1‑149 (MQ=255)
                                                                     gCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGc                                                                                 <  2:628172/149‑1 (MQ=255)
                                                                      cacaATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCa                                                                                <  2:339273/149‑1 (MQ=255)
                                                                         aaTGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCatat                                                                             >  4:895617/1‑149 (MQ=255)
                                                                         aaTGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCatat                                                                             >  3:733896/1‑149 (MQ=255)
                                                                         aaTGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCatat                                                                             >  3:246685/1‑149 (MQ=255)
                                                                              cTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcgg                                                                        <  3:369820/149‑1 (MQ=255)
                                                                              cTGATAAGATGCGTCAAGCATCGCATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcgg                                                                        >  3:237346/1‑149 (MQ=38)
                                                                               tGATGCGAAGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCCCCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggc                                                                       >  4:517828/1‑149 (MQ=255)
                                                                                 atgcgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggcgg                                                                     <  1:469053/149‑1 (MQ=255)
                                                                                 ataagatgcGTCAAGCATCGCATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggcgg                                                                     >  1:55959/5‑149 (MQ=38)
                                                                                  tgcgatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCggcggc                                                                    >  3:423716/1‑149 (MQ=255)
                                                                                  taagatgcGTCAAGCATCGCATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCggcggc                                                                    >  2:234351/4‑149 (MQ=38)
                                                                                     gatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagc                                                                 >  4:320142/1‑149 (MQ=255)
                                                                                     gatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagc                                                                 >  3:17521/1‑149 (MQ=255)
                                                                                     gatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagc                                                                 >  1:258716/1‑149 (MQ=255)
                                                                                      atgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcagca                                                                >  1:137061/1‑149 (MQ=255)
                                                                                      atgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCCGGAATATGCCagcggcgtga                                                                <  3:78588/149‑5 (MQ=38)
                                                                                          ttGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCAcgccg                                                                                      <  4:35691/123‑1 (MQ=255)
                                                                                          ttGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCAcgccg                                                                                      >  3:35691/1‑123 (MQ=255)
                                                                                          ttGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCgg                                                            <  1:199579/149‑1 (MQ=255)
                                                                                                cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTaa                                                      <  2:558697/149‑1 (MQ=255)
                                                                                                cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTaa                                                      <  1:879507/149‑1 (MQ=255)
                                                                                                 atcgcaTCAGGCATACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                                     <  3:333441/144‑1 (MQ=255)
                                                                                                 aTCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcag                                                                  <  4:694231/136‑1 (MQ=255)
                                                                                                 aTCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgcag                                                                  >  3:694231/1‑136 (MQ=255)
                                                                                                 aTCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                                     <  3:895617/149‑1 (MQ=255)
                                                                                                 aTCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                                     <  2:578550/149‑1 (MQ=255)
                                                                                                 aTCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAg                                                     <  2:202871/149‑1 (MQ=255)
                                                                                                  tCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAgg                                                    <  4:645050/149‑1 (MQ=255)
                                                                                                     tATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAAt                                                                                             >  1:229272/1‑105 (MQ=255)
                                                                                                     tATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAAt                                                                                             <  2:229272/105‑1 (MQ=255)
                                                                                                      aTCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTc                                                <  1:111457/149‑1 (MQ=255)
                                                                                                       tCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTcc                                               <  4:246685/149‑1 (MQ=255)
                                                                                                             cTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTc                                         >  1:699958/1‑149 (MQ=255)
                                                                                                              tACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCg                                        <  2:674796/149‑1 (MQ=255)
                                                                                                                cAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCa                                      <  1:662587/149‑1 (MQ=255)
                                                                                                                   tGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTc                                   <  1:200211/149‑1 (MQ=255)
                                                                                                                   tGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTc                                   >  2:321925/1‑149 (MQ=255)
                                                                                                                    gggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCa                                  <  2:777148/149‑1 (MQ=255)
                                                                                                                     ggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCa                                      >  1:818199/1‑144 (MQ=255)
                                                                                                                     ggTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCa                                      <  2:818199/144‑1 (MQ=255)
                                                                                                                      gTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGa                                >  3:606605/1‑149 (MQ=255)
                                                                                                                         ccGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGAcgc                             >  2:887067/1‑149 (MQ=255)
                                                                                                                           ggATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGa                           >  2:932737/1‑149 (MQ=255)
                                                                                                                           ggATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGa                           >  1:241305/1‑149 (MQ=255)
                                                                                                                              tCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAtgt                        <  4:863227/149‑1 (MQ=255)
                                                                                                                                ggTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATgtgt                      >  4:113896/1‑149 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGc                                                   >  1:815520/1‑119 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGc                                                   <  2:815520/119‑1 (MQ=255)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGa                           >  3:475954/3‑141 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTc               <  2:899926/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCgg             <  1:21454/149‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGt            >  1:201332/1‑149 (MQ=255)
                                                                                                                                              ggCGTTTCCGCCGCATCCGGCAAGAATAGAGCACCAGTTAACAGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCaaa        <  2:932413/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGGCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTc     >  2:511634/1‑149 (MQ=255)
                                                                                                                                                  tttACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAGTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCa    <  3:832092/149‑1 (MQ=255)
                                                                                                                                                  tttACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCt                  <  1:907985/135‑1 (MQ=255)
                                                                                                                                                  tttACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCt                  >  2:907985/1‑135 (MQ=255)
                                                                                                                                                    tACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGa  <  2:55959/149‑1 (MQ=255)
                                                                                                                                                    tACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGCCGGTCAAATTCAGa  <  4:900587/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGA  >  NC_000913/551348‑551644

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