Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,518,715 T→G 5.2% K269Q (AAG→CAG)  clpA ← ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,518,7150TG 5.2% 244.9 / 4.0 77K269Q (AAG→CAG) clpAATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base G (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

CCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACGCCTGC  >  NC_002947/4518568‑4518862
                                                                                                                                                   |                                                                                                                                                    
ccGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGACCAGCGCCTTGAAGCGCTtc                                                                                                                                                     >  5:32375/1‑149 (MQ=255)
   acgacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCg                                                                                                                                                  <  3:254207/149‑1 (MQ=255)
   acgacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGAAGCGCTTCAAGGCGCTGCTCg                                                                                                                                                  <  4:317344/149‑1 (MQ=255)
   acgacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGAAGCGCTTCAAGGCGCTGCTCg                                                                                                                                                  <  3:376517/149‑1 (MQ=255)
    cgacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGa                                                                                                                                                 <  6:339389/149‑1 (MQ=255)
    cgacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGa                                                                                                                                                 <  5:352446/149‑1 (MQ=255)
     gacAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGaa                                                                                                                                                <  4:177943/149‑1 (MQ=255)
        agcagcGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGt                                                                                                                                              >  2:271899/1‑148 (MQ=255)
         gaagcGTCGTGAGCAGTTTGGAGGCATCCATAACCCCGCACGACGCCGAACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGACTTGAAGCGCTTCTCGAAGTca                                                                                                                                            <  4:307850/147‑1 (MQ=255)
          cagcGACTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcac                                                                                                                                           >  7:276874/1‑149 (MQ=255)
            gcGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacca                                                                                                                                         <  5:265921/149‑1 (MQ=255)
             cGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcaccac                                                                                                                                        >  3:311077/1‑149 (MQ=255)
                cTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAACGCCTTGAAGCGCTTCTCGAAGTCACCACGGt                                                                                                                                     <  1:271899/149‑1 (MQ=255)
                   gAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTa                                                                                                                                    <  5:184005/147‑1 (MQ=255)
                   gAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTAtt                                                                                                                                  >  5:284143/1‑149 (MQ=255)
                         gTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGCCGCTTGCGCAGCTCACCGAGCAGCGCCGTGACGCGGTTCTCGAAGTCACCGCGGGAGTTGGTAc                                                                                                                            >  8:292244/1‑149 (MQ=255)
                            ggAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCg                                                                                                                          <  4:284234/148‑1 (MQ=255)
                             gAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCga                                                                                                                         <  6:409783/148‑2 (MQ=255)
                               ggCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCa                                                                                                                      <  8:148410/149‑1 (MQ=255)
                                      atcACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTAc                                                                                                                            <  7:385199/136‑1 (MQ=255)
                                      atcACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTAc                                                                                                                            >  8:385199/1‑136 (MQ=255)
                                      atcACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGgc                                                                                                              >  4:213505/1‑149 (MQ=255)
                                         accatcaccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGcc                                                                                                           >  5:249981/8‑149 (MQ=255)
                                           cccgcccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGcc                                                                                                           >  2:376990/1‑147 (MQ=255)
                                           cccgcccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGcc                                                                                                           >  2:371290/1‑147 (MQ=255)
                                           cccgcccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGcc                                                                                                           >  8:332403/1‑147 (MQ=255)
                                           cccgcccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCa                                                                                                         <  1:236966/149‑1 (MQ=255)
                                            ccgcccgACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGcc                                                                                                           >  7:165207/1‑146 (MQ=255)
                                                  gACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAg                                                                                                  >  5:322480/1‑149 (MQ=255)
                                                  gACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAg                                                                                                  >  2:287426/1‑149 (MQ=255)
                                                   aCGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAgg                                                                                                 <  6:284143/149‑1 (MQ=255)
                                                     gccgcACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGACCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGa                                                                                               <  6:32375/149‑1 (MQ=255)
                                                      ccgcACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAg                                                                                              <  8:276874/149‑1 (MQ=255)
                                                      ccgcACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAg                                                                                              <  5:207531/149‑1 (MQ=255)
                                                          aCCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTaca                                                                                          >  5:244525/1‑149 (MQ=255)
                                                            cGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTAcacc                                                                                        <  5:273042/149‑1 (MQ=255)
                                                                     tgatgGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTgc                                                                               >  6:385796/1‑149 (MQ=255)
                                                                          gtgtGAATTTCGTCGATGAACAGGATGGTCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCcag                                                                          <  3:321728/149‑1 (MQ=255)
                                                                             tgAATTTCGTCGATGAACAGGATGGCCTGCGTACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCAGCAAGGCGCCCAg                                                                                                        >  1:265568/1‑115 (MQ=255)
                                                                             tgAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCAGCAAGGCGCCCAg                                                                                                        <  2:265568/115‑1 (MQ=255)
                                                                              gAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAgg                                                                      <  7:332403/149‑1 (MQ=255)
                                                                                aTTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTc                                                                    <  3:178893/149‑1 (MQ=255)
                                                                                    cgtcgATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAgg                                                                 <  5:378979/148‑1 (MQ=255)
                                                                                     gtcgATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGTCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCa                                                               >  5:17508/1‑149 (MQ=255)
                                                                                      tcgATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCAc                                                              >  4:116827/1‑149 (MQ=255)
                                                                                      tcgATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCAc                                                              >  4:240128/1‑149 (MQ=255)
                                                                                        gATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCt                                                            >  5:422519/1‑149 (MQ=255)
                                                                                          tGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGa                                                          >  3:421083/1‑149 (MQ=255)
                                                                                          tGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGa                                                          <  4:311077/149‑1 (MQ=255)
                                                                                            aaCAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGAcc                                                        <  3:213505/149‑1 (MQ=255)
                                                                                                    ggCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATcgacga                                                >  3:2387/1‑149 (MQ=255)
                                                                                                     gCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGAt                                               >  6:77307/1‑149 (MQ=255)
                                                                                                     gCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGAt                                               <  2:201538/149‑1 (MQ=255)
                                                                                                       cTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGc                                             <  6:249981/149‑1 (MQ=255)
                                                                                                          cGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGtt                                          >  7:96923/1‑149 (MQ=255)
                                                                                                          cGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGt                                           <  1:165754/148‑1 (MQ=255)
                                                                                                          cGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGt                                           <  3:240128/148‑1 (MQ=255)
                                                                                                            gACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTg                                        <  3:116827/149‑1 (MQ=255)
                                                                                                                 ttGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAgg                                                                                                 <  5:92276/87‑1 (MQ=255)
                                                                                                                 ttGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAgg                                                                                                 >  6:92276/1‑87 (MQ=255)
                                                                                                                   gcgcAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAgg                                  <  8:165207/148‑1 (MQ=255)
                                                                                                                       aGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTc                             >  5:24803/1‑149 (MQ=255)
                                                                                                                       aGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTc                             <  5:252680/149‑1 (MQ=255)
                                                                                                                          tCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGc                          >  3:296449/1‑149 (MQ=255)
                                                                                                                          tCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGc                          >  1:100902/1‑149 (MQ=255)
                                                                                                                            ttcGAGAAGCGCTTCAAGGCGCTGCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGa                        >  3:66599/3‑149 (MQ=255)
                                                                                                                             ccGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGAt                       >  8:388462/1‑149 (MQ=255)
                                                                                                                              cGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAGAGGCGCCCAGATCAAGGGTGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCCGCCCTTCGGCGAt                       >  1:336930/1‑149 (MQ=255)
                                                                                                                              cGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATg                      <  7:292244/149‑1 (MQ=255)
                                                                                                                              cGAGAAGCGCTTCAAGGCGCTGCTCGAAGACACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACAACACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATg                      >  1:33437/1‑149 (MQ=255)
                                                                                                                                  cagcGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCg                  >  2:930/1‑149 (MQ=255)
                                                                                                                                      gcCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCtt              <  8:96923/149‑1 (MQ=255)
                                                                                                                                       ccTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTg             >  8:222226/1‑149 (MQ=255)
                                                                                                                                       ccTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTg             >  5:304682/1‑149 (MQ=255)
                                                                                                                                        cTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGc            <  1:376990/149‑1 (MQ=255)
                                                                                                                                        cTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGc            <  1:371290/149‑1 (MQ=255)
                                                                                                                                        cTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGc            <  7:388462/149‑1 (MQ=255)
                                                                                                                                           gAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCa         >  6:436162/1‑149 (MQ=255)
                                                                                                                                                  ttCTCGAAGTCACCACCGTTTTCGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACgcctgc  <  6:302247/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACGCCTGC  >  NC_002947/4518568‑4518862

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: