Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,288,185 C→G 7.8% R55P (CGG→CCG)  PP_3759 ← CheB‑like methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,288,1850CG 7.8% 323.0 / 15.0 102R55P (CGG→CCG) PP_3759CheB‑like methylesterase
Reads supporting (aligned to +/- strand):  ref base C (46/48);  new base G (5/3);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 108 total aligned reads displayed.

CGGGCGCGAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATCACCACT  >  NC_002947/4288037‑4288331
                                                                                                                                                    |                                                                                                                                                   
cGGGCGCGAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAgccg                                                                                                                                                    >  2:1892790/1‑150 (MQ=255)
    cgcgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCcgccgc                                                                                                                                                <  2:424534/150‑1 (MQ=255)
    cgcgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCAACAGCCAACACATCGAAGCACCCGCGGCTGATGCTCTCCTGGTCGCGTGCCTCGCATACCGGCCGATGCAGCcgccgc                                                                                                                                                <  1:1486631/150‑1 (MQ=255)
      cgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCACGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTg                                                                                                                                             <  2:229331/150‑1 (MQ=255)
       gagaAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTgc                                                                                                                                            >  1:720906/1‑150 (MQ=255)
          aaGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGACCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCTgcatc                                                                                                                                         <  2:1324498/150‑3 (MQ=255)
             tGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACg                                                                                                                                      >  2:1825935/1‑150 (MQ=255)
               cACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACgcc                                                                                                                                    <  1:575484/150‑1 (MQ=255)
               cACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACgcc                                                                                                                                    <  2:189757/150‑1 (MQ=255)
                  gggTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACgccagc                                                                                                                                 >  1:613002/1‑150 (MQ=255)
                           cTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCt                                                                                                                        <  1:1179103/150‑1 (MQ=255)
                             ggCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGAAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCtgg                                                                                                                      >  1:1190528/1‑149 (MQ=255)
                               cTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTgg                                                                                                                    >  1:268447/1‑150 (MQ=255)
                                    cGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGgtc                                                                                                               >  1:342899/1‑150 (MQ=255)
                                     gACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGgtcg                                                                                                              >  2:1742218/1‑150 (MQ=255)
                                        aGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                           >  1:1455638/1‑150 (MQ=255)
                                           gTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGc                                                                                                        >  1:708613/1‑150 (MQ=255)
                                            tCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCa                                                                                                       <  1:964607/150‑1 (MQ=255)
                                               ggTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGt                                                                                                    <  1:369536/150‑1 (MQ=255)
                                                 tCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgc                                                                                                  <  2:431075/150‑1 (MQ=255)
                                                 tCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgc                                                                                                  <  2:1632792/150‑1 (MQ=255)
                                                     gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTg                                                                                                                                             <  1:1187587/103‑1 (MQ=255)
                                                     gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTg                                                                                                                                             >  2:1187587/1‑103 (MQ=255)
                                                     gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgcagca                                                                                              >  2:1537428/1‑150 (MQ=255)
                                                      cTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCAc                                                                                             <  1:474092/150‑1 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCTgca                                                                                                                                           <  2:2354119/103‑1 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCTgca                                                                                                                                           <  2:2354095/103‑1 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCTgca                                                                                                                                           >  1:2354095/1‑103 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCTgca                                                                                                                                           >  1:2354119/1‑103 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACg                                                                                            >  2:1890638/1‑150 (MQ=255)
                                                       tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACg                                                                                            >  1:1436628/1‑150 (MQ=255)
                                                          aCCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCag                                                                                         >  2:1534807/1‑150 (MQ=255)
                                                           ccGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                        <  1:1367700/150‑1 (MQ=255)
                                                           ccGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                        <  1:1532830/150‑1 (MQ=255)
                                                            cGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGTCTGTGGc                                                                                                                   <  1:210625/122‑1 (MQ=255)
                                                            cGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGc                                                                                                                   >  2:210625/1‑122 (MQ=255)
                                                                  gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                           >  1:2435124/1‑124 (MQ=255)
                                                                  gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                           <  2:2435124/124‑1 (MQ=255)
                                                                  gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                 <  2:1361337/150‑1 (MQ=255)
                                                                  gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGACTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                 >  2:1369387/1‑150 (MQ=255)
                                                                   tGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGt                                                                                >  1:1789687/1‑150 (MQ=255)
                                                                        aCCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGgtc                                                                                                               <  2:1160637/114‑1 (MQ=255)
                                                                        aCCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGgtc                                                                                                               >  1:1160637/1‑114 (MQ=255)
                                                                          ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                         >  2:829223/1‑150 (MQ=255)
                                                                          ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                         >  2:496068/1‑150 (MQ=255)
                                                                          ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                         >  1:1141034/1‑150 (MQ=255)
                                                                           cGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGAAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGaa                                                                        <  1:2276684/150‑1 (MQ=255)
                                                                                ccGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGt                                                                                                                             >  2:1462221/1‑92 (MQ=255)
                                                                                ccGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGt                                                                                                                             <  1:1462221/92‑1 (MQ=255)
                                                                                 cGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGa                                                                  >  1:1632723/1‑150 (MQ=255)
                                                                                 cGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGa                                                                  >  1:96002/1‑150 (MQ=255)
                                                                                  gCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGAt                                                                 <  2:1141034/150‑1 (MQ=255)
                                                                                   ccACATAGATGTGCCCTCTTCTGTTGCTCTCCTTGTCGCGTGCCACGCATACCGTCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATg                                                                >  2:1910940/1‑150 (MQ=255)
                                                                                     acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGccagccag                                                                                                                              <  2:29950/86‑1 (MQ=255)
                                                                                     acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGccagccag                                                                                                                              >  1:29950/1‑86 (MQ=255)
                                                                                     acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGc                                                              <  1:1892790/150‑1 (MQ=255)
                                                                                     acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGc                                                              >  2:871457/1‑150 (MQ=255)
                                                                                        tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                           <  2:2120755/150‑1 (MQ=255)
                                                                                        tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                           <  1:1537428/150‑1 (MQ=255)
                                                                                        tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                           <  2:1789687/150‑1 (MQ=255)
                                                                                        tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                           <  2:708613/150‑1 (MQ=255)
                                                                                        tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                           <  2:342899/150‑1 (MQ=255)
                                                                                           ctCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTg                                                        <  1:636952/149‑1 (MQ=255)
                                                                                             cTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGcc                                                      <  1:1890638/150‑1 (MQ=255)
                                                                                               gCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                    <  1:926015/150‑1 (MQ=255)
                                                                                                   gctgctGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                        <  1:1120580/110‑1 (MQ=255)
                                                                                                   gctgctGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                        >  2:1120580/1‑110 (MQ=255)
                                                                                                         gATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGaa                                          <  2:837357/150‑1 (MQ=255)
                                                                                                           tGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACa                                        >  2:399714/1‑150 (MQ=255)
                                                                                                           tGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACa                                        >  2:2125674/1‑150 (MQ=255)
                                                                                                            gCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACaa                                       >  1:2313958/1‑150 (MQ=255)
                                                                                                            gCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACaa                                       >  1:2207759/1‑150 (MQ=255)
                                                                                                             ctctGCTTGTCGCGTGCCTCGCATAGCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTAGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCCGCGTGCCGAGCACCTGGAACAAt                                      <  2:1190528/150‑1 (MQ=255)
                                                                                                             ctctCCTTGTCGCGTGCCTCGCATACCGGCCGCTGAAGCCCCCGGTAGAAGCGGGCTTGCAGGTGGCTGTGGCGGTCTTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      <  1:1742218/150‑1 (MQ=255)
                                                                                                             ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      <  2:1455638/150‑1 (MQ=255)
                                                                                                             ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      <  1:256688/150‑1 (MQ=255)
                                                                                                             ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      <  2:720906/150‑1 (MQ=255)
                                                                                                                tcCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                    >  2:183309/1‑133 (MQ=255)
                                                                                                                tcCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                    <  1:183309/133‑1 (MQ=255)
                                                                                                                  cTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      <  2:2190147/145‑1 (MQ=255)
                                                                                                                  cTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                      >  1:2190147/1‑145 (MQ=255)
                                                                                                                        gcgTGCCTCGCATACCGGCCGATGTAGCCGCCGCTGCAGCACGCCATCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCa                                                 <  2:1253787/128‑1 (MQ=255)
                                                                                                                        gcgTGCCTCGCATACCGGCCGATGTAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCa                                                 >  1:1253787/1‑128 (MQ=255)
                                                                                                                          gTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                 >  2:2118373/1‑94 (MQ=255)
                                                                                                                          gTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCgccacgccac                         <  1:829223/150‑1 (MQ=255)
                                                                                                                           tGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAACACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCAcc                        >  1:2305097/1‑150 (MQ=255)
                                                                                                                              cTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGc                     >  2:878393/1‑150 (MQ=255)
                                                                                                                               tCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGcc                    >  1:1140473/1‑150 (MQ=255)
                                                                                                                                 gCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGa                  <  1:1910940/150‑1 (MQ=255)
                                                                                                                                      ccGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGgcgc             <  1:1936366/150‑1 (MQ=255)
                                                                                                                                       cGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGcc            <  2:613002/150‑1 (MQ=255)
                                                                                                                                        ggtcGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAACACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCg           >  1:391987/4‑150 (MQ=255)
                                                                                                                                           tgctGCAGCGGCGGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATc        >  1:863522/4‑150 (MQ=255)
                                                                                                                                            gctGCAGCGGCGGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATca       >  2:2425811/3‑150 (MQ=255)
                                                                                                                                               gCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                    <  2:1298494/102‑1 (MQ=255)
                                                                                                                                               gCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                    >  1:1298494/1‑102 (MQ=255)
                                                                                                                                               gCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATcacca    <  2:1632723/150‑1 (MQ=255)
                                                                                                                                                 aGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGACCCATGCCGCCACGCTATAGGCCGTGCCGCTAATCTCCACt  <  2:1969532/150‑1 (MQ=255)
                                                                                                                                                   ccgccgCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGc                                                                           <  1:768501/75‑1 (MQ=255)
                                                                                                                                                   ccgccgCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGc                                                                           >  2:768501/1‑75 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGGCGCGAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCTGCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATCACCACT  >  NC_002947/4288037‑4288331

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