Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,786,022 T→A 6.2% intergenic (+38/‑148) rpsB → / → tsf 30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,786,0220TA 6.2% 233.3 / 2.5 78intergenic (+38/‑148)rpsB/tsf30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base A (3/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGA  >  NC_002947/1785875‑1786168
                                                                                                                                                    |                                                                                                                                                  
tCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCt                                                                                                                                                    <  7:130718/148‑1 (MQ=255)
 cGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCttt                                                                                                                                                  >  7:206107/1‑149 (MQ=255)
   aGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                <  4:355357/149‑1 (MQ=255)
     cTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccc                                                                                                                                              >  1:317515/1‑149 (MQ=255)
       gTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccccc                                                                                                                                            >  1:232400/1‑149 (MQ=255)
          cATGACTTCGATGGCTGACTCATTCATCGGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGGACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACTCAAAAAGGGGGCTTTGCCCCCttt                                                                                                                                         <  7:266390/149‑1 (MQ=255)
           aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCtttt                                                                                                                                        <  7:58191/149‑1 (MQ=255)
            tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAGCGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttttt                                                                                                                                       <  4:399126/149‑1 (MQ=255)
            tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttttt                                                                                                                                       <  4:337779/149‑1 (MQ=255)
            tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttttt                                                                                                                                       >  7:216359/1‑149 (MQ=255)
            tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                       <  1:234730/149‑1 (MQ=255)
                  cGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAcc                                                                                                                                 >  2:75169/1‑149 (MQ=255)
                       gCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGaa                                                                                                                            <  3:253016/149‑1 (MQ=255)
                            cGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCct                                                                                                                       <  4:115328/149‑1 (MQ=255)
                               aGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAAGGGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgc                                                                                                                     >  1:337553/1‑149 (MQ=255)
                               aGTCATCCGCGGCCGAAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGAGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                    <  4:17793/149‑1 (MQ=255)
                                    tCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAg                                                                                                               >  2:155404/1‑149 (MQ=255)
                                      cgcgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAg                                                                                                               >  3:242659/1‑147 (MQ=255)
                                      cgcgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGc                                                                                                              <  1:182826/148‑1 (MQ=255)
                                        cgGCCGCAAGATCGTTGCCGGCGGCACCGATGTTTTCGCTGAAG‑AAGCGGCTGCACCTGCTGCTGTGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATc                                                                                                           <  6:309273/149‑1 (MQ=255)
                                         ggCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCg                                                                                                          <  1:155404/149‑1 (MQ=255)
                                             gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                      >  5:39326/1‑149 (MQ=255)
                                             gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                      >  7:138392/1‑149 (MQ=255)
                                                acaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGc                                                                                                   <  8:206107/149‑1 (MQ=255)
                                                 caacGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGca                                                                                                  >  6:392069/1‑149 (MQ=255)
                                                     gTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatcat                                                                                              >  6:206499/1‑149 (MQ=255)
                                                      ttGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATg                                                                                             <  1:261762/149‑1 (MQ=255)
                                                        gCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATgg                                                                                            <  6:70095/148‑1 (MQ=255)
                                                        gCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATgg                                                                                            <  5:206499/148‑1 (MQ=255)
                                                          cggcggCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGcc                                                                                         >  8:66600/1‑149 (MQ=255)
                                                          aggcggcACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGcc                                                                                         <  7:365572/148‑1 (MQ=255)
                                                            gcggcACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctg                                                                                       <  5:392069/149‑1 (MQ=255)
                                                                 aCCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTc                                                                                                                 <  7:83618/118‑1 (MQ=255)
                                                                 aCCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTc                                                                                                                 >  8:83618/1‑118 (MQ=255)
                                                                        tttACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                           >  2:132135/1‑149 (MQ=255)
                                                                        tttACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                           <  8:42949/149‑1 (MQ=255)
                                                                         ttACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                           <  2:195641/148‑1 (MQ=255)
                                                                             gCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCACCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCg                                                                      >  5:155374/1‑149 (MQ=255)
                                                                             gCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCAAGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCg                                                                      >  7:82843/1‑149 (MQ=255)
                                                                                 aag‑aagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGa                                                                    <  8:216359/147‑1 (MQ=255)
                                                                                 aag‑aagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAtt                                                                  <  5:45713/149‑1 (MQ=255)
                                                                                   g‑aagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                <  1:132135/149‑1 (MQ=255)
                                                                                      agCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGc                                                               <  4:344135/148‑1 (MQ=255)
                                                                                      agCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                              <  2:337553/149‑1 (MQ=255)
                                                                                       gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatgaaa                                                                            <  4:399314/134‑1 (MQ=255)
                                                                                       gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatgaaa                                                                            <  4:405692/134‑1 (MQ=255)
                                                                                       gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatgaaa                                                                            >  3:405692/1‑134 (MQ=255)
                                                                                       gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatgaaa                                                                            >  3:399314/1‑134 (MQ=255)
                                                                                       gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                >  2:345082/1‑146 (MQ=255)
                                                                                           cTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                         <  2:232400/149‑1 (MQ=255)
                                                                                             gCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAATGGGGCTTTGCCCCCTTTTTGACACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTCATGAAATGAAAACAGCGGATTTGAAGAATTg                                                        >  6:177699/1‑148 (MQ=255)
                                                                                               aCCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAc                                                     >  7:30728/1‑149 (MQ=255)
                                                                                                      gctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGt                                               >  5:79207/1‑148 (MQ=255)
                                                                                                       ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGa                                             >  1:277792/1‑149 (MQ=255)
                                                                                                       ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGAGGATTTGCAGAATTGAACGCCCGTGa                                             >  7:372346/1‑149 (MQ=255)
                                                                                                        tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                            >  3:89090/1‑149 (MQ=255)
                                                                                                         gAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                           >  6:259732/1‑149 (MQ=255)
                                                                                                              aTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACgg                                      <  4:89090/149‑1 (MQ=255)
                                                                                                                      tAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATg                              >  3:15383/1‑149 (MQ=255)
                                                                                                                       aGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCGTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGATTgg                             >  1:286634/1‑149 (MQ=255)
                                                                                                                       aGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATgg                             <  5:1254/149‑1 (MQ=255)
                                                                                                                          gTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCATAATTGAACGCCCGTTATTGACGGGTGGAATGGTTg                          <  8:138392/149‑1 (MQ=255)
                                                                                                                          gTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                           >  1:2276/1‑148 (MQ=255)
                                                                                                                           tCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGaa                                                           <  1:396356/115‑1 (MQ=255)
                                                                                                                           tCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGaa                                                           >  2:396356/1‑115 (MQ=255)
                                                                                                                           tCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTg                          >  7:179104/1‑148 (MQ=255)
                                                                                                                           tCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTg                          >  5:241856/1‑148 (MQ=255)
                                                                                                                           tCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGa                         <  1:345082/149‑1 (MQ=255)
                                                                                                                               cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGaaaa                      <  6:39326/148‑1 (MQ=255)
                                                                                                                               cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGaa                        <  2:277792/146‑1 (MQ=255)
                                                                                                                                     aCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAAGCCCGTGATTGCCGGGTGGAATGGTTGAAAAACTTTc                >  7:399593/1‑148 (MQ=255)
                                                                                                                                         aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa             >  6:212228/1‑147 (MQ=255)
                                                                                                                                         aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCCGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga           <  1:75169/149‑1 (MQ=255)
                                                                                                                                         aaaaaGGGGGCAAAGCCCCCTTTTTGGCACCTTGAAATCCTGCTGTAAGCATCGGCCCCGCATCATGGGCCGGCTGAAATGAAAAAAGCGGGTTTGCAGAATTTAACGCCCGGGGTTGACGGGTGGAATGGTTGTAAAACTTTACaa             >  8:181251/1‑147 (MQ=255)
                                                                                                                                         aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa             >  4:365174/1‑147 (MQ=255)
                                                                                                                                          aaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag          >  5:357693/1‑149 (MQ=255)
                                                                                                                                          aaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAACACAGCGGATTTGCAGAATTGAAAGCCCGTGATGGACGGGTGGAATGGTGGAAAAACGTGCCAagag          >  6:89079/1‑149 (MQ=255)
                                                                                                                                           aaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag          >  6:288258/1‑148 (MQ=255)
                                                                                                                                            aaGGGGGCCTTGCCCCCTTTTTTCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGTGCCTGCTGAAATTAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGa        <  2:313569/149‑1 (MQ=255)
                                                                                                                                                  gCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGa  <  6:79207/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAG‑AAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGA  >  NC_002947/1785875‑1786168

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