Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,585,549 A→C 12.5% intergenic (‑156/+15) treSB ← / ← glgE bifunctional trehalose synthase B/maltokinase/alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,585,5490AC12.5% 334.4 / 5.4 145intergenic (‑156/+15)treSB/glgEbifunctional trehalose synthase B/maltokinase/alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (98/25);  new base C (18/4);  total (116/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

TCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTC  >  NC_002947/4585401‑4585698
                                                                                                                                                      |                                                                                                                                                     
tCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGcac                                                                                                                                                      >  2:5888032/1‑150 (MQ=255)
  cTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAg                                                                                                                                                    >  1:8065890/1‑150 (MQ=255)
  cTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGGTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGcaccg                                                                                                                                                    >  1:4179753/1‑148 (MQ=255)
  cTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACAGCGGTGCCGAGGTGCAGGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCCCTCCTACCGGGGGGGGGGcaccg                                                                                                                                                    >  1:8065995/1‑148 (MQ=255)
  cTTTCAATTGAATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTAACCGCGAAGATGACACAGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGAGGATAAACCCGCTCCTACGCGGGGGGGGGcccaa                                                                                                                                                    >  2:7546488/1‑149 (MQ=38)
   tttCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACAg                                                                                                                                                    >  1:3139341/1‑150 (MQ=255)
   tttCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAgg                                                                                                                                                   >  1:985312/1‑149 (MQ=255)
       aaTTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGa                                                                                                                                               >  2:2234234/1‑150 (MQ=255)
         ttGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAAc                                                                                                                                             >  2:204346/1‑150 (MQ=255)
           gcatgcaCAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCGTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCAGTCTGAACCt                                                                                                                                           >  1:3815806/1‑150 (MQ=255)
                 acaGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAGGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCGTGTTaa                                                                                                                                     >  1:8380940/1‑150 (MQ=255)
                 acaGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACCGTCTGAACCGTGTTaa                                                                                                                                     >  2:3103082/1‑150 (MQ=255)
                 acaGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAAc                                                                                                                                    >  2:4509171/1‑150 (MQ=255)
                   aGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGCTAAAACCGCTCCTACGGGGGGGGGGGGCACAGTCTGAACCTTGTTAAc                                                                                                                                    >  1:5531584/1‑150 (MQ=255)
                   aGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTCCGGGGGGGGGGGGCCCAGGCGTAACCTTGTTAAc                                                                                                                                    >  1:962130/1‑150 (MQ=38)
                       aTCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCCGTCTGAACCGTGTTAACGGGCt                                                                                                                               >  2:2706576/1‑150 (MQ=255)
                       aTCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCt                                                                                                                               >  2:6209587/1‑150 (MQ=255)
                        tCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCtt                                                                                                                              >  2:2611104/1‑150 (MQ=255)
                        tCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCATGTTAACGGGCtt                                                                                                                              >  1:8505760/1‑150 (MQ=255)
                           gAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACCGTCGGAACCGTGTTAACGGGCtttt                                                                                                                            >  1:6713378/1‑150 (MQ=38)
                            aGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCAGGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTc                                                                                                                          >  1:4535232/1‑150 (MQ=255)
                               ttGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTAGAt                                                                                                                       >  1:5943490/1‑150 (MQ=255)
                                 gCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCAGGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCGACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATgc                                                                                                                     >  1:4280335/1‑150 (MQ=255)
                                       tCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGa                                                                                                               >  2:256517/1‑150 (MQ=255)
                                          CCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATg                                                                                                             <  1:5386558‑M1/150‑1 (MQ=255)
                                          CCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGGCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGc                                                                                                            <  1:2898391‑M1/150‑1 (MQ=255)
                                           CCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGcc                                                                                                           <  2:858440‑M1/150‑1 (MQ=255)
                                           CCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGcc                                                                                                           <  2:3391087‑M1/150‑1 (MQ=255)
                                              TAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGTGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGaa                                                                                                        <  1:2985999‑M1/150‑1 (MQ=255)
                                                GGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACg                                                                                                      <  1:7153607‑M1/150‑1 (MQ=255)
                                                 GAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACgg                                                                                                     <  2:4069966‑M1/150‑1 (MQ=255)
                                                 GAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACgg                                                                                                     <  2:4069976‑M1/150‑1 (MQ=255)
                                                 GAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCCCCGTCTGGACCGTGTTAACGGGCTTTTTCGTTGCGCCAGAAGCCGAACGGc                                                                                                    <  1:1166975‑M1/150‑1 (MQ=255)
                                                     GGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAg                                                                                                  <  1:1298054‑M1/150‑1 (MQ=255)
                                                       GTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCTCGGATAAACCCGCTCCTACGGGGGGGGGGGCAAAGTCAGAAACATGTTAAAGGGCTTTTTCGATGTGCCATATGCCGAACGGCAGGTg                                                                                               <  2:6320718‑M1/150‑1 (MQ=255)
                                                       GTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTg                                                                                               <  1:1393324‑M1/150‑1 (MQ=255)
                                                        gtttatccgcgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTg                                                                                               >  1:8522594‑M2/87‑150 (MQ=255)
                                                             ccGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAggg                                                                                         >  1:9007906/1‑150 (MQ=255)
                                                              cgcgAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAggg                                                                                         >  1:1708168/1‑150 (MQ=255)
                                                              cgcgAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAggg                                                                                         >  1:1707889/1‑150 (MQ=255)
                                                               gcgAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGGCGGAACCGGGGTAACGGGCGTTTGCGATGCGCCCCATGCCCAACGGCCGGGGCCCGGGcg                                                                                       >  2:1556261/1‑149 (MQ=25)
                                                                 cgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGctc                                                                                      >  1:267719‑M2/78‑150 (MQ=255)
                                                                    ccatgcacCTCGGCACCGCGGTGTATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATg                                                                                  >  2:2026271/6‑150 (MQ=38)
                                                                     tatgcacCTCGGCACCGCGGTGCATAGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATgcg                                                                                >  2:2878032/5‑150 (MQ=38)
                                                                     aTGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATg                                                                                  >  1:7146278/1‑150 (MQ=255)
                                                                       gACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATgcgc                                                                               >  1:1903001/1‑150 (MQ=255)
                                                                       gACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTACCAGGGCTCTATgcgc                                                                               >  1:5955262/1‑150 (MQ=255)
                                                                          acCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTAGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATc                                                                            >  2:7554570/1‑150 (MQ=255)
                                                                           ccGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATcc                                                                           >  1:1858419/1‑150 (MQ=255)
                                                                           ccGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATcc                                                                           >  1:1858349/1‑150 (MQ=255)
                                                                            ctcggcaccGCGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCAt                                                                         >  1:1002935/8‑150 (MQ=38)
                                                                            ctcggcaccGCGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCAAAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCAt                                                                         >  1:1003056/8‑150 (MQ=38)
                                                                            cgcgGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCa                                                                          >  1:379147/1‑150 (MQ=255)
                                                                             gcgGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCAtt                                                                        >  2:7504943/1‑150 (MQ=255)
                                                                                  gCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGt                                                                    >  2:7942894/1‑150 (MQ=255)
                                                                                     gAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCtt                                                                 >  2:4443454/1‑150 (MQ=255)
                                                                                       ggTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTg                                                                >  1:4533787/1‑150 (MQ=255)
                                                                                        gTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGCACCGTCGGCACCTTGGTAACGGGGGTTTTCGATGCGCCCGTTGCCGAACGGCAGGTGCCCGGGCGCTATGCGCATCCCTTTGGTATTGc                                                               >  2:3415938/1‑150 (MQ=25)
                                                                                          gCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCGGAACCTTGTTAACGGGCTTTTTCGATTCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCGATGCGCATCCATTGGGTCTTGCCGt                                                            >  2:1498679/1‑150 (MQ=255)
                                                                                           acctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacggggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGt                                                            >  2:6855072‑M2/52‑150 (MQ=255)
                                                                                            aTGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTAc                                                          >  1:9124968/1‑150 (MQ=255)
                                                                                             tGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTAcc                                                         >  1:210119/1‑150 (MQ=255)
                                                                                                     gCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATTCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGc                                                >  1:314880/1‑150 (MQ=255)
                                                                                                            cGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGc                                          >  1:8746065/1‑150 (MQ=255)
                                                                                                            cGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGGCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCGATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGc                                          >  1:8746207/1‑150 (MQ=255)
                                                                                                              gtgtTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGt                                       >  2:8927426/1‑150 (MQ=255)
                                                                                                               tgtTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTCGAAGCGCCAGATGCCGAACGGCAGGTGCCAGGGCGCGATGCGCCTCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGt                                       >  2:505903/1‑150 (MQ=255)
                                                                                                                gtTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGtt                                      <  1:5207416/150‑1 (MQ=255)
                                                                                                                    gcgGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGGTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTtcatc                                  >  2:2745005/1‑150 (MQ=255)
                                                                                                                    gcgGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATTCGCCAGAAGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTtcatc                                  >  1:3783351/1‑150 (MQ=255)
                                                                                                                       gATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAgg                               >  1:60660/1‑150 (MQ=255)
                                                                                                                       gATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAgg                               >  1:60681/1‑150 (MQ=255)
                                                                                                                        aTAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGGTCATCAGGt                              >  1:7625802/1‑150 (MQ=255)
                                                                                                                         tAAACCCGCTCCTACGGGGGGGGGGGCCCAGTCTGAACCGTGTTAACGGGCTTTTTCGATGCCCCAGATGCCGAACGGCACGTGCCAGGGCACGATGCGCCTCCAATGGGTGTTGCCGTACCACGCACAGCCGTGGCCGTTCATCAGGTc                             >  2:7956703/1‑150 (MQ=255)
                                                                                                                          taaacccgctcctacgggggGGGGGGCCCAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCAGTGGGTCTTGCCGTACCACGACCACCGGGGGCCGTTCCTCAGGTc                             >  2:5474453‑M2/21‑150 (MQ=255)
                                                                                                                              ccgctcctacggggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGAAGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCa                         >  2:6721409‑M2/17‑150 (MQ=255)
                                                                                                                              cccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCa                         >  2:844537‑M2/17‑150 (MQ=255)
                                                                                                                              cccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCa                         >  2:6397216‑M2/17‑150 (MQ=255)
                                                                                                                              cccgctcctacgggggGGGGGGCACAGTCTGAACCTTGGTAACGGGGTTTTTCGATGCGCCAGATGCCGAAAGGCAGGTGCCAGGGCTCGATGCGCATCCAATGGGTCTTGCCCTACCACGACCAGCGGTGGCCGTTCATCAGGtcttaa                         >  1:4064282‑M2/17‑148 (MQ=255)
                                                                                                                               ccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCAc                        >  1:5096636‑M2/16‑150 (MQ=255)
                                                                                                                               ccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCAc                        >  2:6151900‑M2/16‑150 (MQ=255)
                                                                                                                               ccgctcctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTAGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGGTCATCAGGTCTTCAc                        >  1:5770260‑M2/16‑150 (MQ=255)
                                                                                                                                ccgctcctacggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCAcc                       >  2:939436‑M2/15‑150 (MQ=255)
                                                                                                                                   ctacgggggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTg                    >  1:3985730‑M2/12‑150 (MQ=255)
                                                                                                                                   ccccccgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTg                    <  2:5943490‑M2/139‑1 (MQ=255)
                                                                                                                                    ccccggggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTgg                   <  2:8065995‑M2/140‑1 (MQ=255)
                                                                                                                                    cccccgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTgg                   <  2:8065890‑M2/140‑1 (MQ=255)
                                                                                                                                     ctacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTggg                  >  1:5327922‑M2/10‑150 (MQ=255)
                                                                                                                                      tacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGt                 >  2:8823325‑M2/9‑150 (MQ=255)
                                                                                                                                      tacgggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTACTATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCAACCATTGGGTCTTGCCGTACCACCACCAGCGGTGGCCGTTCAACAGGTCTTCACCCTGGGt                 >  1:7645379‑M2/9‑150 (MQ=255)
                                                                                                                                        cggcggGGGGGGCCCCGCCTGCACCGTGGGATCGGGCTCGGTCGAGGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTCTGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGttt               <  2:8704711‑M2/144‑1 (MQ=255)
                                                                                                                                           gggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTc                          >  2:5303111‑M2/4‑136 (MQ=255)
                                                                                                                                           gggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCgg            <  1:6666953‑M2/147‑1 (MQ=255)
                                                                                                                                             gGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCg          <  1:4442490‑M2/149‑1 (MQ=255)
                                                                                                                                                 gggCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATCCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgt      >  1:6355532/1‑150 (MQ=255)
                                                                                                                                                  ggCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgtc     <  1:4507692/150‑1 (MQ=255)
                                                                                                                                                   gCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgtcg    >  1:7354501/1‑150 (MQ=255)
                                                                                                                                                   gCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgtcg    >  1:7354596/1‑150 (MQ=255)
                                                                                                                                                      acaGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgtcgtc  <  2:8505760/150‑1 (MQ=255)
                                                                                                                                                      acaGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTgtcgtc  >  1:1152511/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTC  >  NC_002947/4585401‑4585698

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