Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,155,407 C→G 8.4% G300A (GGC→GCC)  PP_2771 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,155,4070CG 8.4% 366.2 / 18.8 117G300A (GGC→GCC) PP_2771sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (47/59);  new base G (5/5);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 120 total aligned reads displayed.

ATTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGCA  >  NC_002947/3155259‑3155552
                                                                                                                                                    |                                                                                                                                                 
aTTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTTAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGgccgcc                                                                                                                                                  >  1:2035383/1‑150 (MQ=255)
 ttGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGccgccg                                                                                                                                                 <  2:2391264/150‑1 (MQ=255)
   gACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGc                                                                                                                                               >  2:568062/1‑150 (MQ=255)
    accaccAGCTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTAGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGcc                                                                                                                                              <  2:1784710/150‑1 (MQ=255)
      caccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTg                                                                                                                                            >  1:1904771/1‑150 (MQ=255)
      caccaGTTCGCCGATGGCGATCAGTCGGGCGCGAGTAACGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTg                                                                                                                                            <  1:738191/150‑1 (MQ=255)
       accaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGc                                                                                                                                           <  1:1889205/150‑1 (MQ=255)
       accaGTTCGCCGATGGCGATCAGTAGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGACAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGc                                                                                                                                           >  1:932082/1‑150 (MQ=255)
        ccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGAACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCCCGTGCCGGCCGCCGGCCTGcc                                                                                                                                          >  1:1030616/1‑149 (MQ=255)
          aGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCCGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTTGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGctcc                                                                                                                                        <  1:2168042/150‑1 (MQ=255)
               gCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCCGGTTGCAGTGGCAGAACCGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGaa                                                                                                                                   <  1:1328888/150‑1 (MQ=255)
                ccGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAAc                                                                                                                                  <  1:1008554/150‑1 (MQ=255)
                ccGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAAc                                                                                                                                  >  1:585334/1‑150 (MQ=255)
                 cGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTTGCAGAAACGATAGGCACTGTTTATCAGACGGTCCTCGACTTTCTCTTACAGGCTTCCTGCCTGCTCCCCGcaca                                                                                                                                 >  2:508110/1‑150 (MQ=255)
                   gaccggATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGcc                                                                                                                                              >  2:688261/6‑135 (MQ=255)
                    tGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGgc                                                                                                                              <  2:1904771/150‑1 (MQ=255)
                      gCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGgcgc                                                                                                                            <  1:1280685/150‑1 (MQ=255)
                      gCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGgcgc                                                                                                                            <  2:1920514/150‑1 (MQ=255)
                        gATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTg                                                                                                                          <  1:287052/150‑1 (MQ=255)
                        gATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGAGCAGGCCGGCGGCCTGcac                                                                                                                                         <  2:49175/135‑3 (MQ=255)
                        gATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGAGCAGGCCGGCGGCCTGcac                                                                                                                                         >  1:49175/1‑133 (MQ=255)
                            aGTCGGGCGCGGGGGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAg                                                                                                                      <  2:911767/150‑1 (MQ=255)
                             gTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAgg                                                                                                                     >  2:797659/1‑150 (MQ=255)
                                gggCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTc                                                                                                                  <  2:1771175/150‑1 (MQ=255)
                                gggCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTc                                                                                                                  <  1:1382287/150‑1 (MQ=255)
                                 ggCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCg                                                                                                                 <  1:289779/150‑1 (MQ=255)
                                      gggTGATGTTGCGGCTCGGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCgcagca                                                                                                            >  1:1927117/1‑150 (MQ=255)
                                       ggTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCaa                                                                                                           <  1:568062/150‑1 (MQ=255)
                                       ggTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCaa                                                                                                           >  2:1971059/1‑150 (MQ=255)
                                         tGATGTTGCGGCTCAGGCCCAGCACGCTCGCCGTGATCACCTTGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACt                                                                                                         <  1:1362143/150‑1 (MQ=255)
                                          gATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTg                                                                                                        <  2:496502/150‑1 (MQ=255)
                                             gTTGCGGCTTAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttc                                                                                                     <  2:2035383/150‑1 (MQ=255)
                                                 cGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGCTAGGAACAAACCAGAGGGCTGTCCTCGACCTTCTCGGGCAGGCCGCCGGCCTGATCCTCGCACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCg                                                                                                 <  1:172700/150‑1 (MQ=255)
                                                     tCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttcttc                                                                                                  <  1:1083497/145‑1 (MQ=255)
                                                     tCAGGCCCAGCAGGCTGGCAGTGTGAACCGGGTTGCAGTGGCAGAAACGATATGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCTCTGAAAGGCTCGCAGCAACTGTTCTTCGgcgc                                                                                             >  2:1001784/1‑150 (MQ=255)
                                                         gCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTc                                                                                         >  1:931957/1‑150 (MQ=255)
                                                         gCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTc                                                                                         >  1:471451/1‑150 (MQ=255)
                                                           ccAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTcgg                                                                                       >  1:1902353/1‑150 (MQ=255)
                                                             agcagGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTcggcg                                                                                     >  1:1294294/1‑150 (MQ=255)
                                                              gcagGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACGGTTCTTCGGCGCTTTcggcgg                                                                                    >  2:216726/1‑150 (MQ=255)
                                                                  gCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCACCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCt                                                                                >  1:2318902/1‑150 (MQ=255)
                                                                    tGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAg                                                                                                                      >  1:1508473/1‑110 (MQ=255)
                                                                    tGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCa                                                                              <  2:585334/150‑1 (MQ=255)
                                                                    tGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCa                                                                              >  1:1756984/1‑150 (MQ=255)
                                                                    tGGCCATGTGCACATGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCa                                                                              >  2:1956297/1‑150 (MQ=255)
                                                                        cGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCgg                                                                          >  2:631548/1‑150 (MQ=255)
                                                                           gtgCACCTGGTTGCAGTGGCAGAAACGATAGGCACTCAGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGa                                                                       <  1:1109345/150‑1 (MQ=255)
                                                                            tgCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGACGGCCTGCTCATCGAACAGGCGGTGAAAGGCTCGCAGCAACTGTTCTTCCGCGCTTTCCGCGATGCTCACCGGCGAc                                                                      <  2:1428129/150‑1 (MQ=255)
                                                                              cACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCGTCGGCGCTTTCGGCGGTGCTCACCGGCGAttt                                                                    <  2:699588/150‑1 (MQ=255)
                                                                                ccTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGCCTg                                                                                                                                            <  2:156567/76‑1 (MQ=255)
                                                                                ccTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGCCTg                                                                                                                                            >  1:156567/1‑76 (MQ=255)
                                                                                  tGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTcc                                                                >  2:391525/1‑150 (MQ=255)
                                                                                  tGGTTGCAGTGGCAGAAACGACAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTcc                                                                <  2:932082/150‑1 (MQ=255)
                                                                                   ggTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCt                                                               <  2:1902353/150‑1 (MQ=255)
                                                                                    gTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCAGGTCAACTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTg                                                              <  1:762091/150‑1 (MQ=255)
                                                                                        cAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCg                                                                                                 <  1:2408390/111‑1 (MQ=255)
                                                                                               aGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATgc                                                   >  1:785196/1‑150 (MQ=255)
                                                                                                  aaCGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCa                                                <  2:444129/150‑1 (MQ=255)
                                                                                                   aCGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTg                                                                                                                                            <  2:1362111/57‑1 (MQ=255)
                                                                                                     gATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGa                                                                       <  1:2238397/124‑1 (MQ=255)
                                                                                                     gATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGa                                                                       >  2:2238397/1‑124 (MQ=255)
                                                                                                       gaGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTNTNGNCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGc                                           <  1:1338797/149‑1 (MQ=255)
                                                                                                        aGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAggacaa                                          <  2:2204358/150‑5 (MQ=255)
                                                                                                          gCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGAc                                        <  1:391525/150‑1 (MQ=255)
                                                                                                            aCTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAg                                      <  2:1756984/150‑1 (MQ=255)
                                                                                                            aCTACGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAg                                      >  1:552733/1‑150 (MQ=255)
                                                                                                              tGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGAc                                                                      >  2:1814730/1‑116 (MQ=255)
                                                                                                              tGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGAc                                                                      <  1:1814730/116‑1 (MQ=255)
                                                                                                              tGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGgc                                    >  1:734726/1‑150 (MQ=255)
                                                                                                                cGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAACGGATTACAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGgcgc                                  <  1:1956297/150‑1 (MQ=255)
                                                                                                                  tAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAg                                <  2:1864043/150‑1 (MQ=255)
                                                                                                                  tAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAg                                >  2:2221219/1‑150 (MQ=255)
                                                                                                                     cAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAg                                               >  1:1086986/1‑132 (MQ=255)
                                                                                                                     cAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAg                                               <  2:1086986/132‑1 (MQ=255)
                                                                                                                      aGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGgtcg                            <  2:1137921/150‑1 (MQ=255)
                                                                                                                        gCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGCCTGca                                                                                                                                          >  1:889020/1‑37 (MQ=25)
                                                                                                                        gCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGCCTGca                                                                                                                                          <  2:889020/38‑2 (MQ=37)
                                                                                                                        gCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCg                                                                                                 <  1:566873/79‑1 (MQ=255)
                                                                                                                        gCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGTTGCGACAGGCCCAGgtcgtc                          >  2:503275/1‑150 (MQ=255)
                                                                                                                            tCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTcc                      <  2:1375492/150‑1 (MQ=255)
                                                                                                                            tCCGCGACCTTCGCGTGCAGGCCGCCGGCCGGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCCGCGGTGCTCCCCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTcc                      <  2:1030616/150‑1 (MQ=255)
                                                                                                                             ccTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCt                     <  2:952403/150‑1 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCgcag                                                                                                              >  2:1998937/1‑60 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCgcag                                                                                                              <  1:1998937/60‑1 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGgtcgtcgtc                       <  2:437783/147‑1 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGgtcgtcgtc                       >  1:437783/1‑147 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTg                    >  1:2135822/1‑150 (MQ=255)
                                                                                                                               tCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTTGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGc                   >  2:1582783/1‑150 (MQ=255)
                                                                                                                               tCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGc                   >  2:1960440/1‑150 (MQ=255)
                                                                                                                               tCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGc                   >  1:1598134/1‑150 (MQ=255)
                                                                                                                                cGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCa                  <  1:797659/150‑1 (MQ=255)
                                                                                                                                    cTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGgcgc                                                          >  1:1814683/1‑106 (MQ=255)
                                                                                                                                    cTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTg                    >  1:340310/1‑144 (MQ=255)
                                                                                                                                    cTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTg                    <  2:340310/144‑1 (MQ=255)
                                                                                                                                           tGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTg       >  1:1857987/1‑150 (MQ=255)
                                                                                                                                           agCAGGCCGGCGGCCTGCTCCTCGAACAGGCGCTGAAAGGATCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTg       >  1:2297845/2‑150 (MQ=255)
                                                                                                                                            gCAGGCCGGCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGc                                                                                 <  1:451128/75‑1 (MQ=255)
                                                                                                                                            gCAGGCCGGCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGc                                                                                 >  2:451128/1‑75 (MQ=255)
                                                                                                                                              aGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCgg    <  2:552733/150‑1 (MQ=255)
                                                                                                                                                gccgccGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGAAGGGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGca  >  1:1355618/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGCA  >  NC_002947/3155259‑3155552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: