Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,786,024 T→A 14.2% intergenic (+40/‑146) 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,0240TA14.2% 187.8 / 18.3 76intergenic (+40/‑146)rpsB/tsf30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (39/26);  new base A (6/5);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785877‑1786171
                                                                                                                                                    |                                                                                                                                                    
gAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTg                                                                                                                                                     >  6:47948/1‑149 (MQ=255)
  gCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGcc                                                                                                                                                   <  4:265984/149‑1 (MQ=255)
  gCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGcc                                                                                                                                                   >  7:308776/1‑149 (MQ=255)
      tACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCt                                                                                                                                               >  1:11874/1‑149 (MQ=255)
      tACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCt                                                                                                                                               >  3:87083/1‑149 (MQ=255)
        cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttt                                                                                                                                             >  4:104198/1‑149 (MQ=255)
         aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCtttt                                                                                                                                            >  8:276280/1‑149 (MQ=255)
         aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCtttt                                                                                                                                            >  6:24237/1‑149 (MQ=255)
          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                           <  4:346208/149‑1 (MQ=255)
          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                           <  2:115524/149‑1 (MQ=255)
          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                           <  8:116730/149‑1 (MQ=255)
           gACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTg                                                                                                                                          >  1:225178/1‑149 (MQ=255)
              ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTATACTTGGCACGCAAAAAGGGGGCTTTTCCACCTTTTTGcaa                                                                                                                                       >  7:333177/1‑147 (MQ=255)
                       tGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAAt                                                                                                                              >  4:339776/1‑149 (MQ=255)
                       tGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAAt                                                                                                                              >  6:362921/1‑149 (MQ=255)
                           gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                          <  2:58509/149‑1 (MQ=255)
                           gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                          >  3:203785/1‑149 (MQ=255)
                           gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                          >  1:96489/1‑149 (MQ=255)
                               tCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGt                                                                                                                      <  5:47948/149‑1 (MQ=255)
                               tCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGt                                                                                                                      <  5:24237/149‑1 (MQ=255)
                               tCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                       <  4:317051/132‑3 (MQ=255)
                               tCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                       >  3:317051/1‑130 (MQ=255)
                                           gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                          >  3:304602/1‑149 (MQ=255)
                                             aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGGAGAAGCGGCTGCACATGCTGCAGAGTAATAGGCGCTAGCGTCAACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                           >  5:384093/1‑114 (MQ=255)
                                                      gCCGGCGGCACCGTAGTTTACGGTGCAGATGCGTCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCGCTTTTTGCCACCTTGAAATCCTTCTGTCAGCATCGGCCCCGCATCATgg                                                                                                <  8:333177/148‑1 (MQ=255)
                                                      gCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATgg                                                                                                <  2:96489/148‑1 (MQ=255)
                                                       ccGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGTGc                                                                                              >  4:241870/1‑149 (MQ=255)
                                                         ggcggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCct                                                                                            <  4:255404/149‑1 (MQ=255)
                                                           cggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgc                                                                                          <  7:276280/149‑1 (MQ=255)
                                                            ggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                         <  2:225178/149‑1 (MQ=255)
                                                              cACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCatggtg                                                                                        >  2:182072/1‑148 (MQ=255)
                                                               aCCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtga                                                                                       >  2:93766/1‑148 (MQ=255)
                                                                ccGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaa                                                                                     <  2:11874/149‑1 (MQ=255)
                                                                   aaGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatg                                                                                   >  7:286877/1‑148 (MQ=255)
                                                                   aaGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatg                                                                                   >  4:4102/1‑148 (MQ=255)
                                                                     gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACACAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatc                                                                                                    <  4:183553/129‑1 (MQ=255)
                                                                     gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACACAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatc                                                                                                    >  3:183553/1‑129 (MQ=255)
                                                                      tttACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                               >  7:349096/1‑149 (MQ=255)
                                                                              gaagaaGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                      >  2:328967/1‑149 (MQ=255)
                                                                                   agCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                                 <  3:339776/149‑1 (MQ=255)
                                                                                   agCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCCGCTGAAATGAAAACAGCGGATTTGCa                                                                 <  4:304602/149‑1 (MQ=255)
                                                                                      ggCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAg                                                                >  1:51399/1‑148 (MQ=255)
                                                                                       gCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAt                                                             >  1:85093/1‑149 (MQ=255)
                                                                                           cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGa                                                          <  5:362921/148‑1 (MQ=255)
                                                                                           cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGa                                                          <  1:328967/148‑1 (MQ=255)
                                                                                           cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGa                                                          <  8:308776/148‑1 (MQ=255)
                                                                                              ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                        <  2:359665/84‑1 (MQ=255)
                                                                                              ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                        >  1:359665/1‑84 (MQ=255)
                                                                                               tgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGcc                                                     >  4:173715/1‑149 (MQ=255)
                                                                                                     tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                               >  7:78989/1‑149 (MQ=255)
                                                                                                     tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                               >  5:223744/1‑149 (MQ=255)
                                                                                                     tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTAAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                               >  2:178189/1‑149 (MQ=255)
                                                                                                      gAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  1:168805/1‑149 (MQ=255)
                                                                                                       aGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  7:284604/1‑148 (MQ=255)
                                                                                                          aaTATGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGAc                                           <  3:4102/148‑1 (MQ=255)
                                                                                                           aTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTATCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACgg                                         >  4:42717/1‑149 (MQ=255)
                                                                                                              tgctcTAGCGTCAACTTGGCACGCAAAAAGGGGTCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCTCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACTCCCGTGATTGACGGGt                                       <  8:349096/144‑1 (MQ=255)
                                                                                                              ggCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTATCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTg                                      <  3:42717/149‑1 (MQ=255)
                                                                                                                   tAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATg                                 >  3:185256/1‑149 (MQ=255)
                                                                                                                             cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACACGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGaaa                          <  1:93766/148‑1 (MQ=255)
                                                                                                                               tGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACt                      >  5:176785/1‑149 (MQ=255)
                                                                                                                                 gCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACtt                     <  2:316951/148‑1 (MQ=255)
                                                                                                                                  cACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc                   <  4:87083/149‑1 (MQ=255)
                                                                                                                                  cACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc                   <  3:173715/149‑1 (MQ=255)
                                                                                                                                    ggCAAAAAGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTcc                  >  6:319006/2‑148 (MQ=255)
                                                                                                                                      cAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTcc                  >  1:223146/1‑146 (MQ=255)
                                                                                                                                       aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGaa                                                         >  2:267477/1‑106 (MQ=255)
                                                                                                                                       aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  5:62167/1‑147 (MQ=255)
                                                                                                                                       aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  6:57976/1‑147 (MQ=255)
                                                                                                                                       aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga              <  2:168805/149‑1 (MQ=255)
                                                                                                                                       aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGAATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGaa                                                         <  1:267477/106‑1 (MQ=255)
                                                                                                                                       aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  5:34681/1‑147 (MQ=255)
                                                                                                                                       aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  5:343735/1‑147 (MQ=255)
                                                                                                                                       aaaaaGCGGCCTTTGCCCCCTTTTTCCCACCTTGAAAACCTGCTGTCAGCATCCGCCCCGCATCATGGGCCTGCTGAAATCAAAACAGCGGATTTGCAGAATTGAACCCCCGACATTGACGGGTGGAAAGGTTGAAAAACTTTCCaa                >  8:249880/1‑147 (MQ=255)
                                                                                                                                        aaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGt                                       >  6:55644/1‑123 (MQ=255)
                                                                                                                                        aaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGt                                       <  5:55644/123‑1 (MQ=255)
                                                                                                                                                    tGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGaaa   <  3:241870/147‑1 (MQ=255)
                                                                                                                                                    tGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  <  4:185256/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785877‑1786171

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: