Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,886,982 T→C 5.9% H292R (CAC→CGC)  PP_3430 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,886,9820TC 5.9% 199.3 / 2.1 67H292R (CAC→CGC) PP_3430two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CCGCCAACTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATCGCCTGCACCTGCTCTTGTG  >  NC_002947/3886836‑3887129
                                                                                                                                                    |                                                                                                                                                      
ccgccAACTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTgg                                                                                                                                                      >  6:101024/1‑149 (MQ=255)
ccgccAACTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTgg                                                                                                                                                      <  1:310274/149‑1 (MQ=255)
 cgccAACTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGt                                                                                                                                                     <  4:306523/149‑1 (MQ=255)
     aaCTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGt                                                                                                                                                  >  4:315277/1‑148 (MQ=255)
     aaCTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTg                                                                                                                                                 <  2:84209/149‑1 (MQ=255)
           cccACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCGTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCGGCCGGCAAACCTCGCACTCCCACGAATCGAGTGGTAGTGACCATc                                                                                                                                           >  4:251646/1‑149 (MQ=255)
           cccACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCACGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATc                                                                                                                                           >  3:364501/1‑149 (MQ=255)
            ccaccaCCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCCCGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATcg                                                                                                                                          >  5:261613/1‑149 (MQ=255)
            ccaccaCCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGAACTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCAAACGAAACCAGTGGTAGTGACAATc                                                                                                                                           >  1:111136/1‑149 (MQ=255)
            ccaccaCCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATcg                                                                                                                                          >  5:165200/1‑149 (MQ=255)
                             tCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAt                                                                                                                         >  6:188015/1‑149 (MQ=255)
                             tCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTCATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAt                                                                                                                         >  3:312454/1‑149 (MQ=255)
                               tcTGTTGTTCATCATCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGTCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAt                                                                                                                         <  8:335762/146‑1 (MQ=255)
                                cTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAtctc                                                                                                                      <  2:111136/149‑1 (MQ=255)
                                      tgCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAgctg                                                                                                                >  5:6046/1‑149 (MQ=255)
                                         aTCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCg                                                                                                                            <  1:218385/134‑1 (MQ=255)
                                         aTCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCg                                                                                                                            >  2:218385/1‑134 (MQ=255)
                                         aTCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTc                                                                                                             >  6:207788/1‑149 (MQ=255)
                                          tCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTcc                                                                                                            >  7:105447/1‑149 (MQ=255)
                                             gcgcCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCcgcg                                                                                                         >  4:118103/1‑149 (MQ=255)
                                               gcCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTg                                                                                                       <  1:149588/149‑1 (MQ=255)
                                                ccTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTgg                                                                                                      >  7:263645/1‑149 (MQ=255)
                                                        tGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGCGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                                               >  2:349299/1‑148 (MQ=255)
                                                        tGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGCGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                                               >  2:343600/1‑148 (MQ=255)
                                                          aCGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTgg                                                                                             <  3:28267/148‑1 (MQ=255)
                                                          aCGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGCGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                                               <  1:343600/146‑1 (MQ=255)
                                                          aCGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGCGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                                               <  1:349299/146‑1 (MQ=255)
                                                          aCGCGCAGCCGCTTCACGGTCTGCGGTGCTGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGGCAAACCTCGCACTCCAACGAACCCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTg                                                                                              <  8:121924/147‑1 (MQ=255)
                                                                 gCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCACGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGccc                                                                                     >  2:274669/1‑149 (MQ=255)
                                                                   ctctTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCgg                                                                                   <  3:315277/149‑1 (MQ=255)
                                                                      ttCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGa                                                                                  <  5:101024/147‑1 (MQ=255)
                                                                      ttCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGa                                                                                  <  6:6046/147‑1 (MQ=255)
                                                                          cGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑Gc                                                                         >  5:138759/1‑149 (MQ=255)
                                                                           ggTATGCGTTGATGTCCCGGCCAACCGCCGTCAGGAGGCTTCCGTCAAACCTCGCACTCCAACGAAACCCGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑Gcc                                                                        <  4:296241/149‑1 (MQ=255)
                                                                                   ggTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                                                        >  1:42583/1‑113 (MQ=255)
                                                                                   ggTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                                                        <  2:42583/113‑1 (MQ=255)
                                                                                   ggTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTa                                                                 <  7:120480/149‑1 (MQ=255)
                                                                                   ggTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTa                                                                 <  5:188015/149‑1 (MQ=255)
                                                                                     tGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACt                                                               <  5:207788/149‑1 (MQ=255)
                                                                                            ccGGCCAACCGCCGTCAGGAGGCTGCCGTACAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAt                                                                                                                         >  1:218905/1‑88 (MQ=255)
                                                                                            ccGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAt                                                                                                                         <  2:218905/87‑1 (MQ=255)
                                                                                            ccGGCCAACCGCCGTCAGGAGGCTGCAGTCAAACCTCGCACTCCAACGAAACCAGTCGTAGTGACCATCGCGGTGGCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTa                                                         <  2:286983/148‑1 (MQ=255)
                                                                                             cGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACg                                                       >  6:340987/1‑149 (MQ=255)
                                                                                                  aaCCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTg                                                  >  1:54688/1‑149 (MQ=255)
                                                                                                     agccgTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGCAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAATGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGCTACGCTCTGTCg                                               >  6:293242/2‑149 (MQ=255)
                                                                                                       ccgTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCGCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCggg                                             <  4:312454/149‑1 (MQ=255)
                                                                                                          tCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGc                                          <  5:29004/149‑1 (MQ=255)
                                                                                                            aggtggCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCa                                         <  8:263645/148‑1 (MQ=255)
                                                                                                                   tGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTc                                  >  2:8751/1‑148 (MQ=255)
                                                                                                                    gCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGa                                                                                  >  1:129503/1‑102 (MQ=255)
                                                                                                                    gCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGa                                                                                  <  2:129503/102‑1 (MQ=255)
                                                                                                                    gCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAAt                                                                     <  2:70680/112‑1 (MQ=255)
                                                                                                                    gCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAAt                                                                     >  1:70680/1‑112 (MQ=255)
                                                                                                                       gTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTACTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCa                             >  7:700/1‑149 (MQ=255)
                                                                                                                          cAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCAcc                           >  6:381120/1‑149 (MQ=255)
                                                                                                                               cTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACAt                      >  5:135176/1‑149 (MQ=255)
                                                                                                                               cTCGCACTCCAACGAAACCAGTGGTAGTGACCATCACGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACAt                      <  1:274669/149‑1 (MQ=255)
                                                                                                                                tCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCAc                                        >  8:17607/1‑130 (MQ=255)
                                                                                                                                tCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCAc                                        <  7:17607/130‑1 (MQ=255)
                                                                                                                                tCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATc                     >  8:92180/1‑149 (MQ=255)
                                                                                                                                    aCTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAATGGGGGGGGAACTAGGAGGAATAGAAATTCGATACGCTCTGTCGGGTGCACCAGCACAATCCCCGAAATcc                    >  1:121060/1‑148 (MQ=255)
                                                                                                                                        cAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACTAGCTCCATCACCGACATCGCCTGCAc             <  1:8751/149‑1 (MQ=255)
                                                                                                                                              aaCCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTGCAGCTGCTCCGCGTGGCTGCATTGTACCAGCCCGGACACT‑‑CCCAATCGTACTTCTGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATCGCCTGCACCTGctc       <  7:163040/149‑1 (MQ=255)
                                                                                                                                              aaCCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGCCACT‑‑GCCAATCGTACTTCGGTCCGCTCTGTCGGGTGCACCAGCTCCATCCCCGCCAACGCCTGCCCCTGctc       >  8:190456/1‑149 (MQ=255)
                                                                                                                                               acccttGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCCCGTGGCTGCATTGGAACAGCCCTGACACT‑‑GCCAATCGTACTTCGGTACGGTCTGTCGGGGGCACCAGGTCCATCACGTACATCTCCTGCACCTGctct      <  5:293242/144‑1 (MQ=255)
                                                                                                                                                  aGTGGTAGTTACCATCGCGGTGTCGCAGGCGGATGTCCTGCTTCTCCGCGTGGCTGCATTGGAACAGCCCTGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATCGCCTGCACCTGCTCTtgt   <  6:165200/149‑1 (MQ=255)
                                                                                                                                                  aGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATCGCCTGCACCTGCTCTtgt   <  3:118103/149‑1 (MQ=255)
                                                                                                                                                   gTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAAGAGCCCTGACACT‑‑GCCAATCGTACTTCTGTACGCTCTGTCGGGTGAACCAGCTCCATCACCGACATCGCCTGCACCTGCTCTtgtg  <  4:364501/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCGCCAACTGGCCCACCACCTCCATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACT‑‑GCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGACATCGCCTGCACCTGCTCTTGTG  >  NC_002947/3886836‑3887129

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