Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,842 A→G 5.9% L413L (TTG→CTG)  PP_3453 ← sensor protein RstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,8420AG 5.9% 203.9 / 2.3 68L413L (TTG→CTG) PP_3453sensor protein RstB
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base G (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATC  >  NC_002947/3912697‑3912989
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                                     <  6:59225/122‑1 (MQ=255)
                         cAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAg                                                                                                                                                     >  5:59225/1‑122 (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)
                                                          tctgTTCGCGGTTGCTGCGCCCCTTTCCGCCGGCGGGCAAACGATTCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCTGGCTGATGGCCAACCCCAGGCCAAAACCAC‑Cg                                                                                        <  4:258756/149‑1 (MQ=255)
                                                              ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           <  1:218921/143‑1 (MQ=255)
                                                              ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           >  2:218921/1‑143 (MQ=255)
                                                              ttCGCGGTTGCTGCGCCCCTATCCGCCCGCGGGCAAACGAATCCTGAAATGCGCCCCACACAGCCCCGACTTCGCCAGGGTCAACGTCCCACCCTGCCCCACGAACGCCCCCCGGCTGATGGCCCACCCCAGGCAAAAACCAC‑CGgtgt                                                                                    >  6:236388/1‑149 (MQ=255)
                                                                 gcgGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGc                                                                                 >  1:90191/1‑149 (MQ=255)
                                                                    gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           <  1:144066/137‑1 (MQ=255)
                                                                    gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           >  2:144066/1‑137 (MQ=255)
                                                                            gcCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGACCCACCCAGCGCGGATTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCTGCTGATGGCCAACCCCAGGCCAAAAACAC‑CGGTGTTGCGATCGCGACt                                                                       <  5:236388/148‑1 (MQ=255)
                                                                              cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                    <  3:215655/149‑1 (MQ=255)
                                                                              cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                    >  4:215655/1‑149 (MQ=255)
                                                                                cTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGt                                                                  >  1:303204/1‑149 (MQ=255)
                                                                                   tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGTGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                 >  6:96053/1‑147 (MQ=255)
                                                                                   tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                 <  2:303204/147‑1 (MQ=255)
                                                                                     cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  7:106563/149‑1 (MQ=255)
                                                                                     cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  4:255035/149‑1 (MQ=255)
                                                                                           cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  5:255476/149‑1 (MQ=255)
                                                                                           cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  6:27238/149‑1 (MQ=255)
                                                                                                    gAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTg                                               >  8:168263/1‑148 (MQ=255)
                                                                                                       tCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGaa                                            >  6:89211/1‑148 (MQ=255)
                                                                                                         cTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCCAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                         <  1:287135/149‑1 (MQ=255)
                                                                                                              cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACTCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGc                                     <  1:96617/149‑1 (MQ=255)
                                                                                                              cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                        >  1:234734/1‑129 (MQ=255)
                                                                                                              cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                        <  2:234734/129‑1 (MQ=255)
                                                                                                                 cgcCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTcgc                                 >  6:279942/1‑149 (MQ=255)
                                                                                                                   ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTTCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                         <  3:70520/139‑1 (MQ=255)
                                                                                                                   ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGTAAAAc                                         >  4:70520/1‑139 (MQ=255)
                                                                                                                   ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCgcgc                               <  7:168263/149‑1 (MQ=255)
                                                                                                                    cccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGct                              >  1:320005/1‑149 (MQ=255)
                                                                                                                     ccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctc                             <  5:279942/149‑1 (MQ=255)
                                                                                                                      cacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcc                            >  4:189842/1‑149 (MQ=255)
                                                                                                                         ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       >  5:237055/1‑119 (MQ=255)
                                                                                                                         ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  6:237055/119‑1 (MQ=255)
                                                                                                                         ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGACCGGCTGGAAAACCCGCTCGCGctcctc                          >  5:189642/1‑148 (MQ=255)
                                                                                                                         ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCtcctcc                         <  2:90191/149‑1 (MQ=255)
                                                                                                                         ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCccctcc                         >  8:299590/1‑149 (MQ=255)
                                                                                                                              gcCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAgg                    <  2:320005/149‑1 (MQ=255)
                                                                                                                               ccGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAggg                   >  7:214756/1‑149 (MQ=255)
                                                                                                                                       tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  6:173484/96‑1 (MQ=38)
                                                                                                                                       tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             >  5:173484/4‑99 (MQ=38)
                                                                                                                                       gcCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAAt           >  4:183002/1‑149 (MQ=255)
                                                                                                                                        ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc          <  4:157767/149‑1 (MQ=255)
                                                                                                                                        ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCCCCTCCACAGCGACGCCAATc          <  7:299590/149‑1 (MQ=255)
                                                                                                                                         gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                          >  7:37015/2‑132 (MQ=255)
                                                                                                                                         gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                          <  8:37015/131‑1 (MQ=255)
                                                                                                                                         cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                     >  6:311444/1‑137 (MQ=255)
                                                                                                                                         cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                     <  5:311444/137‑1 (MQ=255)
                                                                                                                                         cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATcc         >  2:118367/1‑149 (MQ=255)
                                                                                                                                              tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa    <  5:96053/149‑1 (MQ=255)
                                                                                                                                              tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa    <  5:89211/149‑1 (MQ=255)
                                                                                                                                              tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACACGGACGCCAATCCCGTCa    >  3:263468/1‑149 (MQ=255)
                                                                                                                                                aaGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCAtc  <  8:214756/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAAGGACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATC  >  NC_002947/3912697‑3912989

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: ‑