Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,970 C→G 5.7% M208I (ATG→ATC) ‡ trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,9700CG 5.7% 709.0 / 23.1 212M208I (ATG→ATC) ‡trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (98/102);  new base G (6/6);  total (104/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 215 total aligned reads displayed.

ACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAG  >  NC_002947/92821‑93119
                                                                                                                                                     |                                                                                                                                                     
aCGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCAcc                                                                                                                                                       <  2:3872231/150‑1 (MQ=255)
  gCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCCCGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCAt                                                                                                                                                     >  2:5065111/1‑150 (MQ=255)
          ctTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGt                                                                                                                                             >  1:5448359/1‑150 (MQ=255)
           ttCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTa                                                                                                                                            >  1:4436463/1‑150 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGTAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         >  2:1299403/1‑150 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         <  2:2267999/150‑1 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         <  1:899895/150‑1 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         <  1:1970405/150‑1 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         <  1:5202606/150‑1 (MQ=255)
              gAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCGTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                                                                                         >  2:3205085/1‑150 (MQ=255)
                 gAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgtg                                                                                                                                      <  1:1548853/150‑1 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   <  2:4885907/150‑1 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   <  2:294022/150‑1 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   >  2:5335518/1‑150 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   <  2:1020984/150‑1 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   <  2:4067923/150‑1 (MQ=255)
                    cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                                                                                                   <  1:2180705/150‑1 (MQ=255)
                     tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTggg                                                                                                                                  >  1:2922013/1‑150 (MQ=255)
                     tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTggg                                                                                                                                  >  2:686944/1‑150 (MQ=255)
                     tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTggg                                                                                                                                  >  2:349788/1‑150 (MQ=255)
                      cATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgc                                                                                                                                 >  1:3603580/1‑150 (MQ=255)
                       aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgcc                                                                                                                                >  1:2666445/1‑150 (MQ=255)
                       aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgcc                                                                                                                                >  1:4777603/1‑150 (MQ=255)
                       aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgcc                                                                                                                                >  2:3441316/1‑150 (MQ=255)
                       aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgcc                                                                                                                                >  2:2091085/1‑150 (MQ=255)
                           aaCAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgccagcc                                                                                                                            >  1:4731707/1‑150 (MQ=255)
                              aGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACg                                                                                                                         >  1:5245835/1‑150 (MQ=255)
                                  ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                                                                                                     >  2:5352614/1‑150 (MQ=255)
                                   cATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGTCCAGCCACGGTGca                                                                                                                    <  2:1894921/150‑2 (MQ=255)
                                   cATGGCATTGGAGCCACCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGcc                                                                                                                    <  2:4635414/150‑1 (MQ=255)
                                     tGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCg                                                                                                                                                 <  2:4441423/119‑1 (MQ=255)
                                          ttGGAGCCGCCGCCGACGCAGGCAACCTGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGTAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAgg                                                                                                             <  2:2894118/150‑1 (MQ=255)
                                           tGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGt                                                                                                            <  2:1442586/150‑1 (MQ=255)
                                                  gccgccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGgt                                                                                                     <  2:4743778/150‑1 (MQ=255)
                                                      ccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgc                                                                                                                                 <  1:1198418/118‑1 (MQ=255)
                                                        gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTc                                                                                               >  2:5005589/1‑150 (MQ=255)
                                                          cGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGa                                                                                             <  2:2413971/150‑1 (MQ=255)
                                                           gCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGAc                                                                                            >  1:3818923/1‑150 (MQ=255)
                                                               gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAgg                                                                                                             <  2:3626319/129‑1 (MQ=255)
                                                               gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                        <  1:5070033/150‑1 (MQ=255)
                                                               gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                        >  1:1589770/1‑150 (MQ=255)
                                                               gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                        <  2:425352/150‑1 (MQ=255)
                                                               gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                        <  1:3901813/150‑1 (MQ=255)
                                                                 aaCCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGt                                                                                      >  1:212931/1‑150 (MQ=255)
                                                                    cAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAc                                                                                   >  2:2735886/1‑150 (MQ=255)
                                                                     aGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAcc                                                                                  <  1:2768716/150‑1 (MQ=255)
                                                                       gCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCa                                                                                <  2:5448359/150‑1 (MQ=255)
                                                                         tGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGt                                                                              >  2:597370/1‑150 (MQ=255)
                                                                         tGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGt                                                                              >  2:1120538/1‑150 (MQ=255)
                                                                                      gCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTc                                                                 <  2:455528/150‑1 (MQ=255)
                                                                                            cttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGTGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                           <  1:5376182/150‑1 (MQ=255)
                                                                                             ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCtt                                                                  >  1:2941778/1‑142 (MQ=255)
                                                                                             ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATg                                                          >  2:2338592/1‑150 (MQ=255)
                                                                                              tctcTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCggg                                                                                                                                               <  2:3302273/64‑1 (MQ=255)
                                                                                              tctcTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCggg                                                                                                                                               >  1:3302273/1‑64 (MQ=255)
                                                                                              tctcTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATgg                                                         <  2:5070785/150‑1 (MQ=255)
                                                                                                  ttGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTc                                                     >  2:1601233/1‑150 (MQ=255)
                                                                                                  ttGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTc                                                     >  1:3257484/1‑150 (MQ=255)
                                                                                                    gCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCtt                                                   >  1:1674948/1‑150 (MQ=255)
                                                                                                    gCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCtt                                                   >  1:1834480/1‑150 (MQ=255)
                                                                                                     cAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTTCGGGTGTGGGCCAGCCACg                                                                                                                         <  2:4998592/79‑1 (MQ=255)
                                                                                                     cAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACg                                                                                                                         >  1:4998592/1‑79 (MQ=255)
                                                                                                      aaCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCAc                                                                           >  2:4002166/1‑124 (MQ=255)
                                                                                                      aaCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTc                                                 >  1:1048863/1‑150 (MQ=255)
                                                                                                       aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                <  1:696908/150‑1 (MQ=255)
                                                                                                       aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                <  2:19514/150‑1 (MQ=255)
                                                                                                       aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                <  1:2329811/150‑1 (MQ=255)
                                                                                                        cTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                        <  1:2764044/109‑1 (MQ=255)
                                                                                                             cACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTAc                                          <  2:261241/150‑1 (MQ=255)
                                                                                                              aCGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                           >  2:4826907/1‑132 (MQ=255)
                                                                                                              aCGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                           <  1:4826907/132‑1 (MQ=255)
                                                                                                                gggTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCgg                                       >  1:3178376/1‑150 (MQ=255)
                                                                                                                     tcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcca                                  <  1:4266510/146‑2 (MQ=255)
                                                                                                                      ctTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcgg                                 >  1:127624/1‑150 (MQ=255)
                                                                                                                       tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                           <  1:4047890/123‑1 (MQ=255)
                                                                                                                       tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                >  1:624357/1‑150 (MQ=255)
                                                                                                                        ttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                               >  1:4452232/1‑150 (MQ=255)
                                                                                                                        gtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCcagcc                               <  1:2535909/149‑2 (MQ=255)
                                                                                                                           ccGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGtt                            >  1:2211317/1‑149 (MQ=255)
                                                                                                                            cGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTaa                           <  1:5335518/150‑1 (MQ=255)
                                                                                                                             gatcaccGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTAc                          <  2:3767896/144‑1 (MQ=255)
                                                                                                                                  tCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCa                     <  1:1090417/150‑1 (MQ=255)
                                                                                                                                  tCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGACCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCa                     <  2:4363513/150‑1 (MQ=255)
                                                                                                                                   tgACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCAc                    <  1:393932/149‑1 (MQ=255)
                                                                                                                                   cGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGctt                    >  2:1341355/1‑148 (MQ=255)
                                                                                                                                     aCTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGa                  >  1:1357728/1‑150 (MQ=255)
                                                                                                                                      cTGGAAGTCGCGCACCATGGCGGGGTACTGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAAGGCGTCTTTCAGGCTACCGGTGCCGGCGGTAACTGGCACGAt                 >  2:2143037/1‑150 (MQ=255)
                                                                                                                                      cTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCCTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAt                 >  2:4461264/1‑150 (MQ=255)
                                                                                                                                       tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAgg                                              <  2:2643880/120‑1 (MQ=255)
                                                                                                                                       tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAtc                >  2:1206719/1‑150 (MQ=255)
                                                                                                                                          aaGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTATTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCg             <  1:768565/150‑1 (MQ=255)
                                                                                                                                           aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg            >  2:17241/1‑150 (MQ=255)
                                                                                                                                           aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg            >  2:774676/1‑150 (MQ=255)
                                                                                                                                            ttCGCGCACCATGGCGGGGTACGTGTGTGGGCCAGCCACGTTGCCGATCAGGTAGAAGGTGTCTTCGACTTTGGTAACCCAGTCACGCAGGTCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATATCTgc           <  2:3027330/149‑1 (MQ=255)
                                                                                                                                             tCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                >  2:2955602/1‑112 (MQ=255)
                                                                                                                                             tCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGgcg          >  1:4082154/1‑150 (MQ=255)
                                                                                                                                               gcgcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGcc        >  1:4409882/1‑150 (MQ=255)
                                                                                                                                                 gcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGaa                                                                                                        >  2:3045568/1‑52 (MQ=255)
                                                                                                                                                 gcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGaa                                                                                                        <  1:3045568/52‑1 (MQ=255)
                                                                                                                                                   aCCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCcagc    <  2:2373750/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAG  >  NC_002947/92821‑93119

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: