Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,839 T→G 5.9% T414P (ACC→CCC)  PP_3453 ← sensor protein RstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,8390TG 5.9% 212.8 / 4.4 68T414P (ACC→CCC) PP_3453sensor protein RstB
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CTGGCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCA  >  NC_002947/3912694‑3912987
                                                                                                                                                 |                                                                                                                                                     
cTGGCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTc                                                                                                                                                      >  1:352213/1‑147 (MQ=255)
   gCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAg                                                                                                                                                   <  2:48649/147‑1 (MQ=255)
   gCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGt                                                                                                                                                 <  8:102139/149‑1 (MQ=255)
   gCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGt                                                                                                                                                 <  4:125824/149‑1 (MQ=255)
     aaGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGt                                                                                                                                                 >  2:344961/1‑147 (MQ=255)
     aaGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGt                                                                                                                                                 <  1:344961/147‑1 (MQ=255)
         aCCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCAcc                                                                                                                                           <  2:116547/149‑1 (MQ=255)
                          tCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCaa                                                                                                                                                    <  1:23187/123‑1 (MQ=255)
                          tCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCACCACCCAGCGCCGACTGCGCCAGGGTCaa                                                                                                                                                    >  2:23187/1‑123 (MQ=255)
                            cAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAg                                                                                                                                                   >  5:59225/1‑122 (MQ=255)
                            cAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAg                                                                                                                                                   <  6:59225/122‑1 (MQ=255)
                                gAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTg                                                                                                                                                >  8:350584/1‑121 (MQ=255)
                                gAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTg                                                                                                                                                <  7:350584/121‑1 (MQ=255)
                                     ggCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGCTCGCCCGCCGGCTGATg                                                                                                               >  5:27238/1‑149 (MQ=255)
                                                    gCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGcc                                                                                                 >  3:255035/1‑148 (MQ=255)
                                                       ggCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCaaaa                                                                                             >  2:287135/1‑149 (MQ=255)
                                                        gCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCAGGCTGATGGCCACCCCCAGGCCAAAac                                                                                            >  3:258756/1‑149 (MQ=255)
                                                         cccGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAac                                                                                            >  3:157767/1‑148 (MQ=255)
                                                           cGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                         <  2:352213/149‑1 (MQ=255)
                                                            gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                    >  5:54042/1‑137 (MQ=255)
                                                            gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                    <  6:54042/137‑1 (MQ=255)
                                                             tctgTTCGCGGTTGCTGCGCCCCTTTCCGCCGGCGGGCAAACGATTCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCTGGCTGATGGCCAACCCCAGGCCAAAACCACCg                                                                                      <  4:258756/149‑1 (MQ=255)
                                                                 ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                         >  2:218921/1‑143 (MQ=255)
                                                                 ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                         <  1:218921/143‑1 (MQ=255)
                                                                 ttCGCGGTTGCTGCGCCCCTATCCGCCCGCGGGCAAACGAATCCTGAAATGCGCCCCACACAGCCCCGACTTCGCCAGGGTCAACGTCCCACCCTGCCCCACGAACGCCCCCCGGCTGATGGCCCACCCCAGGCAAAAACCACCGgtgt                                                                                  >  6:236388/1‑149 (MQ=255)
                                                                    gcgGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGc                                                                               >  1:90191/1‑149 (MQ=255)
                                                                       gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                         <  1:144066/137‑1 (MQ=255)
                                                                       gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                         >  2:144066/1‑137 (MQ=255)
                                                                               gcCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGACCCACCCAGCGCGGATTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCTGCTGATGGCCAACCCCAGGCCAAAAACACCGGTGTTGCGATCGCGACt                                                                     <  5:236388/148‑1 (MQ=255)
                                                                                 cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCg                                                                  <  3:215655/149‑1 (MQ=255)
                                                                                 cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCg                                                                  >  4:215655/1‑149 (MQ=255)
                                                                                   cTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGt                                                                >  1:303204/1‑149 (MQ=255)
                                                                                      tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGTGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTc                                                               >  6:96053/1‑147 (MQ=255)
                                                                                      tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTc                                                               <  2:303204/147‑1 (MQ=255)
                                                                                        cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGa                                                           <  7:106563/149‑1 (MQ=255)
                                                                                        cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGa                                                           <  4:255035/149‑1 (MQ=255)
                                                                                              cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                     <  5:255476/149‑1 (MQ=255)
                                                                                              cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                     <  6:27238/149‑1 (MQ=255)
                                                                                                       gAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTg                                             >  8:168263/1‑148 (MQ=255)
                                                                                                          tCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGaa                                          >  6:89211/1‑148 (MQ=255)
                                                                                                            cTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCCAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                       <  1:287135/149‑1 (MQ=255)
                                                                                                                 cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACTCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGc                                   <  1:96617/149‑1 (MQ=255)
                                                                                                                 cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                      >  1:234734/1‑129 (MQ=255)
                                                                                                                 cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                      <  2:234734/129‑1 (MQ=255)
                                                                                                                    cgcCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTcgc                               >  6:279942/1‑149 (MQ=255)
                                                                                                                      ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTTCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                       <  3:70520/139‑1 (MQ=255)
                                                                                                                      ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGTAAAAc                                       >  4:70520/1‑139 (MQ=255)
                                                                                                                      ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCgcgc                             <  7:168263/149‑1 (MQ=255)
                                                                                                                       cccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGct                            >  1:320005/1‑149 (MQ=255)
                                                                                                                        ccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctc                           <  5:279942/149‑1 (MQ=255)
                                                                                                                         cacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcc                          >  4:189842/1‑149 (MQ=255)
                                                                                                                            ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                     <  6:237055/119‑1 (MQ=255)
                                                                                                                            ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                     >  5:237055/1‑119 (MQ=255)
                                                                                                                            ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGACCGGCTGGAAAACCCGCTCGCGctcctc                        >  5:189642/1‑148 (MQ=255)
                                                                                                                            ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCtcctcc                       <  2:90191/149‑1 (MQ=255)
                                                                                                                            ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCccctcc                       >  8:299590/1‑149 (MQ=255)
                                                                                                                                 gcCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAgg                  <  2:320005/149‑1 (MQ=255)
                                                                                                                                  ccGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAggg                 >  7:214756/1‑149 (MQ=255)
                                                                                                                                          tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGa                                                           >  5:173484/4‑99 (MQ=38)
                                                                                                                                          tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGa                                                           <  6:173484/96‑1 (MQ=38)
                                                                                                                                          gcCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAAt         >  4:183002/1‑149 (MQ=255)
                                                                                                                                           ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc        <  4:157767/149‑1 (MQ=255)
                                                                                                                                           ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCCCCTCCACAGCGACGCCAATc        <  7:299590/149‑1 (MQ=255)
                                                                                                                                            gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                        >  7:37015/2‑132 (MQ=255)
                                                                                                                                            gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                        <  8:37015/131‑1 (MQ=255)
                                                                                                                                            cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                   >  6:311444/1‑137 (MQ=255)
                                                                                                                                            cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                   <  5:311444/137‑1 (MQ=255)
                                                                                                                                            cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATcc       >  2:118367/1‑149 (MQ=255)
                                                                                                                                                 tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa  <  5:96053/149‑1 (MQ=255)
                                                                                                                                                 tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa  <  5:89211/149‑1 (MQ=255)
                                                                                                                                                 tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACACGGACGCCAATCCCGTCa  >  3:263468/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGGCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCA  >  NC_002947/3912694‑3912987

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: