Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,665,182 G→T 10.0% intergenic (+35/‑194) metK → / → ligB methionine adenosyltransferase/DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,665,1820GT10.0% 1017.3 / 79.9 319intergenic (+35/‑194)metK/ligBmethionine adenosyltransferase/DNA ligase B
Reads supporting (aligned to +/- strand):  ref base G (137/149);  new base T (12/20);  total (149/170)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 323 total aligned reads displayed.

GGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCAC  >  NC_002947/5665034‑5665331
                                                                                                                                                    |                                                                                                                                                     
    aCTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATc                                                                                                                                                  >  1:739456/1‑150 (MQ=255)
    aCTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATc                                                                                                                                                  <  2:293245/150‑1 (MQ=255)
         ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                                             <  1:3116312/150‑1 (MQ=255)
         ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                                             <  1:5570507/150‑1 (MQ=255)
         ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                                             <  1:4760541/150‑1 (MQ=255)
         ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGAGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAg                                                                                                                                             <  2:1926201/150‑1 (MQ=255)
          gCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCGTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                                            <  1:3224526/150‑1 (MQ=255)
            cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCGTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgggg                                                                                                                                          <  1:2865354/150‑1 (MQ=255)
              tACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCt                                                                                                                                        <  2:680677/150‑1 (MQ=255)
                 tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttt                                                                                                                                     >  1:6453830/1‑150 (MQ=255)
                   cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCTCCTGGGAGCGCACTGAACGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCCGCAGGGGCgtgttt                                                                                                                                   >  1:6850564/1‑150 (MQ=255)
                   cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                                                   <  2:2427168/150‑1 (MQ=255)
                    gcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtg                                                                                                                                  <  1:5236797/150‑1 (MQ=255)
                       gcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                               <  1:1159445/150‑1 (MQ=255)
                       gcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                               <  1:463805/150‑1 (MQ=255)
                       gcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                               <  1:5825864/150‑1 (MQ=255)
                       gcagAAGACTGTCGGCGACGACACCTTCACCACCGTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                               <  1:868844/150‑1 (MQ=255)
                         agaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                                             <  1:512224/150‑1 (MQ=255)
                         agaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                                             <  2:3634846/150‑1 (MQ=255)
                         agaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                                             >  2:1511754/1‑150 (MQ=255)
                         agaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                                             <  2:5631582/150‑1 (MQ=255)
                          gaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGGCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCCGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTGTTGTGCCTg                                                                                                                            >  2:2518231/1‑150 (MQ=255)
                                tGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACgg                                                                                                                      >  2:6581173/1‑150 (MQ=255)
                                   cGGCGACGACACCTTCACCACCTTCCCCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAt                                                                                                                   <  2:1058503/150‑1 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCa                                                                                                                                              <  2:5869809/121‑1 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt                                                                                                                 >  1:3501440/1‑150 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt                                                                                                                 >  2:7419065/1‑150 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt                                                                                                                 >  1:1003900/1‑150 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt                                                                                                                 <  1:2411098/150‑1 (MQ=255)
                                       gacgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTc                                                                                                               <  2:6209826/150‑1 (MQ=255)
                                         cgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                                             <  1:3027986/150‑1 (MQ=255)
                                         cgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                                             <  1:3986452/150‑1 (MQ=255)
                                             accttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGATCGCTCGGCAgggg                                                                                                                                          >  1:1274435/1‑117 (MQ=255)
                                              ccttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGccc                                                                                                        <  2:4371302/150‑1 (MQ=255)
                                              ccttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGccc                                                                                                        <  1:1905681/150‑1 (MQ=255)
                                                ttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                                      >  2:5258/1‑150 (MQ=255)
                                                 tcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTc                                                                                                     <  2:2082146/150‑1 (MQ=255)
                                                     caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                 <  1:4896479/150‑1 (MQ=255)
                                                     caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                 >  1:2354491/1‑150 (MQ=255)
                                                     caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                 >  1:1964825/1‑150 (MQ=255)
                                                      accttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATc                                                                                                >  1:2818639/1‑150 (MQ=255)
                                                          tcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATc                                                                                                >  1:3413653/1‑146 (MQ=255)
                                                          tcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATc                                                                                                <  2:3413653/146‑1 (MQ=255)
                                                                gggAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCcgcg                                                                                      >  1:5899967/1‑150 (MQ=255)
                                                                    gcgcACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCg                                                                                  >  1:2254891/1‑150 (MQ=255)
                                                                      gcACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                                                <  2:1342317/150‑1 (MQ=255)
                                                                       cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGttt                                                                               <  2:321889/150‑1 (MQ=255)
                                                                        accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTc                                                                              >  1:5316226/1‑150 (MQ=255)
                                                                        accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTc                                                                              >  2:5426461/1‑150 (MQ=255)
                                                                         ccgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCt                                                                             <  2:5983377/150‑1 (MQ=255)
                                                                          cgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                            <  2:6955682/150‑1 (MQ=255)
                                                                              cgcgCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTc                                                                        >  2:3353137/1‑150 (MQ=255)
                                                                                   ccAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTcc                                                                       <  1:5663461/146‑1 (MQ=255)
                                                                                   ccAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGtct                                                                   >  2:1620388/1‑150 (MQ=255)
                                                                                    cAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGtctc                                                                  >  2:7987/1‑150 (MQ=255)
                                                                                      gTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCgg                                                                <  1:1683896/150‑1 (MQ=255)
                                                                                        cGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGc                                                              <  2:5226169/150‑1 (MQ=255)
                                                                                          ttGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGCCGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCtt                                                            <  2:5514041/150‑1 (MQ=255)
                                                                                           tGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                                           <  1:7220097/150‑1 (MQ=255)
                                                                                               tGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                                       <  1:1305743/150‑1 (MQ=255)
                                                                                               tGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                                       >  1:1482385/1‑150 (MQ=255)
                                                                                                    cTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGc                                                  <  2:4571264/150‑1 (MQ=255)
                                                                                                     tGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCa                                                 >  2:713128/1‑150 (MQ=255)
                                                                                                     tGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCa                                                 <  1:6966915/150‑1 (MQ=255)
                                                                                                        cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa                                              <  2:2018748/150‑1 (MQ=255)
                                                                                                        cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCAGCCGAGCGATCGGCAGGGGCTTATTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa                                              <  1:2466032/150‑1 (MQ=255)
                                                                                                        cGGCCTGTAAGCGTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTGTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa                                              <  2:7350104/150‑1 (MQ=255)
                                                                                                          gCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGc                                            <  1:409762/150‑1 (MQ=255)
                                                                                                            cTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGa                                          <  2:2371592/150‑1 (MQ=255)
                                                                                                            cTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGa                                          <  1:6333740/150‑1 (MQ=255)
                                                                                                             tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaa                                         <  2:7009076/150‑1 (MQ=255)
                                                                                                             tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaa                                         >  2:2816950/1‑150 (MQ=255)
                                                                                                             tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaa                                         >  2:3471378/1‑150 (MQ=255)
                                                                                                              gTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaaa                                        <  1:4227511/150‑1 (MQ=255)
                                                                                                               tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAat                                       >  1:540366/1‑150 (MQ=255)
                                                                                                                   cTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                                   <  2:1687727/150‑1 (MQ=255)
                                                                                                                    ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc                                  >  1:1855578/1‑150 (MQ=255)
                                                                                                                    ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc                                  <  1:6133242/150‑1 (MQ=255)
                                                                                                                    ttCCTGCACCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc                                  <  2:2215387/150‑1 (MQ=255)
                                                                                                                     tCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                            <  2:3113263/107‑1 (MQ=255)
                                                                                                                     tCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAcc                                 <  1:252433/150‑1 (MQ=255)
                                                                                                                        tGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa                              >  1:3815713/1‑150 (MQ=255)
                                                                                                                         gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGAc                             >  2:3486491/1‑150 (MQ=255)
                                                                                                                          cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCtt                                                            >  2:4024489/1‑118 (MQ=255)
                                                                                                                             cgcCAGTCGAAATCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTTACTTTTACGGCATTTTCATCGCCCGTCATATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCTGTCTTGGAGACGAGCTCCAGCGAAATATCACCTGACATCt                         <  2:2789885/150‑1 (MQ=255)
                                                                                                                                cAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgca                      <  2:3868676/150‑1 (MQ=255)
                                                                                                                                     gAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgcac                 >  2:6103104/1‑150 (MQ=255)
                                                                                                                                     gAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgcac                 >  2:1191011/1‑150 (MQ=255)
                                                                                                                                     aaaaGCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgcac                 >  1:980329/2‑150 (MQ=255)
                                                                                                                                       aaGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccg               <  2:4665344/150‑1 (MQ=255)
                                                                                                                                          ccccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgccc            >  1:6423268/1‑150 (MQ=255)
                                                                                                                                           cccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccg           >  2:1746216/1‑150 (MQ=255)
                                                                                                                                            ccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTTGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccgc          >  2:4998163/1‑150 (MQ=255)
                                                                                                                                            ccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccgc          >  1:3896977/1‑150 (MQ=255)
                                                                                                                                            ccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccgc          >  1:6014211/1‑150 (MQ=255)
                                                                                                                                               gCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGc                                                  <  2:2849356/107‑1 (MQ=255)
                                                                                                                                               gCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGcccgcccc       <  1:7008012/150‑1 (MQ=255)
                                                                                                                                               gCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGcccgcccc       <  1:445985/150‑1 (MQ=255)
                                                                                                                                                   aGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCa   >  1:3716327/1‑150 (MQ=255)
                                                                                                                                                    gCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCAc  <  2:2559685/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCAC  >  NC_002947/5665034‑5665331

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