Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,824,715 C→A 5.9% intergenic (+139/+35) PP_5102 → / ← trmB hypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,824,7150CA 5.9% 318.5 / 7.9 101intergenic (+139/+35)PP_5102/trmBhypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base C (50/45);  new base A (4/2);  total (54/47)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-02

Only 100 of 102 total aligned reads displayed.

CGCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTG  >  NC_002947/5824568‑5824862
                                                                                                                                                    |                                                                                                                                                     
cgcgGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCa                                                                                                                                                      <  3:119859/149‑1 (MQ=255)
 gcgGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCac                                                                                                                                                     <  4:323738/149‑2 (MQ=255)
   ggCTGAGCGAATGACTCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAggg                                                                                                                                                   <  2:303374/149‑1 (MQ=255)
    gCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGc                                                                                                                                                 >  2:158870/1‑149 (MQ=255)
       gAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGt                                                                                                                                              <  4:97443/149‑1 (MQ=255)
           gAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                         >  8:173432/1‑149 (MQ=255)
           gAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGttttt                                                                                                                                          >  1:30007/1‑148 (MQ=255)
            aaTGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                         >  1:81556/1‑149 (MQ=255)
            aaTGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAc                                                                                                                                        <  3:258463/149‑1 (MQ=255)
            aaTGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAc                                                                                                                                        <  2:240962/149‑1 (MQ=255)
            aaTGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGtttttt                                                                                                                                         <  2:30007/149‑3 (MQ=255)
                aCCCCTTGTGGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGttttttctt                                                                                                                                      <  8:42779/148‑6 (MQ=255)
                       gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                              <  7:173432/148‑1 (MQ=255)
                       gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGc                                                                                                                               <  2:81556/147‑1 (MQ=255)
                           gAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                              <  8:127000/144‑1 (MQ=255)
                           gAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                              >  7:127000/1‑144 (MQ=255)
                           gAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                          >  1:194860/1‑148 (MQ=255)
                            aGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                          >  8:351724/1‑147 (MQ=255)
                            aGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                          >  2:248607/1‑147 (MQ=255)
                                gATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCTGCCg                                                                                                                    >  7:80964/1‑149 (MQ=255)
                                gATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCg                                                                                                                    <  2:194860/149‑1 (MQ=255)
                                     aaCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAACGGCCGTATCAGCCGCTCaa                                                                                                               >  4:173138/1‑149 (MQ=255)
                                     aaCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCaa                                                                                                               <  1:158870/149‑1 (MQ=255)
                                     aaCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCa                                                                                                                >  7:40204/1‑148 (MQ=255)
                                      aCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAAt                                                                                                              >  3:272994/1‑149 (MQ=255)
                                         gcgAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCACGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCa                                                                                                           >  1:110658/1‑149 (MQ=255)
                                           gAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCAcc                                                                                                         >  7:307307/1‑149 (MQ=255)
                                              aGGGCCGGTTTTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTtc                                                                                                      <  1:316986/149‑1 (MQ=255)
                                                ggccggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTtctc                                                                                                    <  6:363127/149‑1 (MQ=255)
                                                   cggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGaa                                                                                                 <  5:368701/149‑1 (MQ=255)
                                                       ccTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACtt                                                                                              >  2:85126/1‑148 (MQ=255)
                                                           gctgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAg                                                                                          >  2:89309/1‑149 (MQ=255)
                                                            ctgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCaa                                                                                           <  1:98003/148‑1 (MQ=255)
                                                             tgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTc                                                                                        >  3:285219/1‑149 (MQ=255)
                                                               gCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTcc                                                                                       <  5:209372/148‑1 (MQ=255)
                                                                 cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGt                                                                                         >  1:95507/1‑144 (MQ=255)
                                                                 cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGt                                                                                         <  2:95507/144‑1 (MQ=255)
                                                                 cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCca                                                                                     >  8:114410/1‑148 (MQ=255)
                                                                 cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCca                                                                                     >  5:171841/1‑148 (MQ=255)
                                                                  aaTCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGttttt                                                                                                                                          <  6:112096/95‑2 (MQ=18)
                                                                  aaTCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGttttt                                                                                                                                          >  5:112096/1‑94 (MQ=18)
                                                                   aTCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCacac                                                                                  >  4:205177/1‑149 (MQ=255)
                                                                    tctcCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACAcc                                                                                 >  5:351604/1‑149 (MQ=255)
                                                                      tcCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACAcccc                                                                               >  4:259292/1‑149 (MQ=255)
                                                                      tcCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACAcccc                                                                               >  5:345268/1‑149 (MQ=255)
                                                                        cACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCATTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGACGCTCAATCCACCATCTCGAACTTCAAGTCCCACCCCCCGt                                                                             >  6:364946/1‑149 (MQ=255)
                                                                         aCCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTg                                                                            >  8:211748/1‑149 (MQ=255)
                                                                           ccGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGCGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGcc                                                                          >  8:199447/1‑149 (MQ=255)
                                                                                cTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGc                                                                     >  7:114710/1‑149 (MQ=255)
                                                                                   cACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCTATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCg                                                                   <  1:265069/148‑1 (MQ=255)
                                                                                    aCCGGGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGc                                                                  <  8:3295/148‑1 (MQ=255)
                                                                                      cGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCa                                                               <  7:114410/149‑1 (MQ=255)
                                                                                        tGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAAACGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCAcc                                                             >  6:290847/1‑149 (MQ=255)
                                                                                          ggCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAGCGAAAAGGCCGTATCAGACGCTCAATCCACAATCTCGAACTTCAAGTCACACACCCCGTGCACCAGCCGCGCAAcgg                                                           >  4:12124/1‑148 (MQ=38)
                                                                                                gcTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGtt                                                     >  8:57822/1‑149 (MQ=255)
                                                                                                gcTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGtt                                                     <  5:364946/149‑1 (MQ=255)
                                                                                                 cTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTc                                                    <  3:259292/149‑1 (MQ=255)
                                                                                                 cTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTc                                                    <  5:290847/149‑1 (MQ=255)
                                                                                                    ccGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGaa                                                 >  2:245132/1‑149 (MQ=255)
                                                                                                    ccGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGaa                                                 <  4:285219/149‑1 (MQ=255)
                                                                                                     cGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAAc                                                >  8:164816/1‑149 (MQ=255)
                                                                                                     cGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGACCCAGCCGCTCACCGCGGCGTTCGAAc                                                >  4:264068/1‑149 (MQ=255)
                                                                                                       tATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACtt                                              >  8:248093/1‑149 (MQ=255)
                                                                                                         ttCGCGTGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTCTCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                             <  6:138048/148‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTgg                                            >  7:163924/1‑149 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTgg                                            >  7:159035/1‑149 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTgg                                            <  6:345268/149‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                             <  2:110658/148‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                             <  7:164816/148‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                             <  3:176302/148‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                             <  4:272994/148‑1 (MQ=255)
                                                                                                            gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGt                                           >  1:97409/1‑148 (MQ=255)
                                                                                                              gggTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATc                                       <  7:351724/149‑1 (MQ=255)
                                                                                                                gTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCgg                                     <  8:40204/149‑1 (MQ=255)
                                                                                                                   aaCCCGCTCCTACAAAGAAAAAACCTCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGgc                                   <  8:307307/148‑1 (MQ=255)
                                                                                                                    aCCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGgcgc                                 <  6:171841/149‑1 (MQ=255)
                                                                                                                     cccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGct                                >  1:241539/1‑149 (MQ=255)
                                                                                                                      ccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGctc                               >  6:24006/1‑149 (MQ=255)
                                                                                                                          tCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGACCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCg                           >  5:289481/1‑149 (MQ=255)
                                                                                                                            cTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCg                           <  8:114710/147‑1 (MQ=255)
                                                                                                                            cTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCg                           <  1:245132/147‑1 (MQ=255)
                                                                                                                                aaaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAACGGCCGTATCAGCCGCTCACTCCACGTTCTCGAACTTCAAGTCCCCCACCGCGTGGCCCAGCCGCTCACCGCGGCGTTCGAACAGGGTGCGCGGGCGCTCATCGGGgcgc                     >  4:280297/1‑149 (MQ=255)
                                                                                                                                aaaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCCCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgcgc                     >  8:218525/1‑149 (MQ=255)
                                                                                                                                 aaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGcgcg                    >  5:60953/1‑149 (MQ=255)
                                                                                                                                  aGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCgg                   >  6:54836/1‑149 (MQ=255)
                                                                                                                                  aGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCgg                   <  3:205177/149‑1 (MQ=255)
                                                                                                                                  aGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCgg                   >  5:363070/1‑149 (MQ=255)
                                                                                                                                   ggtaaaaCCGCCCTGAAAAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTATTCGGCGCGCGGc                  >  6:308002/4‑149 (MQ=255)
                                                                                                                                   gAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGc                  <  7:49192/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAACCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGca                 >  4:336692/1‑149 (MQ=255)
                                                                                                                                      aaaaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGcaca               <  3:12124/149‑1 (MQ=255)
                                                                                                                                      aaaaCCGCCCTGAAAAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGcaca               >  1:268428/1‑149 (MQ=255)
                                                                                                                                         aCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgcgc                     <  3:54322/140‑1 (MQ=255)
                                                                                                                                         aCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAg            <  6:351604/149‑1 (MQ=255)
                                                                                                                                         aCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACAACCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgcgc                     >  4:54322/1‑140 (MQ=255)
                                                                                                                                           cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGt          <  7:97033/149‑1 (MQ=255)
                                                                                                                                           cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGt          <  7:199447/149‑1 (MQ=255)
                                                                                                                                            gCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTg         >  2:13330/1‑149 (MQ=255)
                                                                                                                                               cTTTTCACGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCg      <  3:148181/149‑1 (MQ=255)
                                                                                                                                                    cAGGGCGGTTTTACCGAAAAGGCCGTATCAACCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTg  <  3:336692/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTG  >  NC_002947/5824568‑5824862

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