Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,229 A→T 7.1% M122K (ATG→AAG) ‡ trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,2290AT 7.1% 591.5 / 23.6 183M122K (ATG→AAG) ‡trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (94/76);  new base T (6/7);  total (100/83)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 185 total aligned reads displayed.

CGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTC  >  NC_002947/93080‑93378
                                                                                                                                                     |                                                                                                                                                     
cGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCa                                                                                                                                                       >  1:805464/1‑150 (MQ=255)
 ggTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCAt                                                                                                                                                      <  1:5318542/150‑1 (MQ=255)
       ggcggTAACTGGCACGATCTCGGCGCCAACCAGCTTCATGCGGAACACGTTCGCCAGCTGGCGCTCGATGTCGGTGGCACCCATGAAGATCACGCAAGGCAGGACGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCTGGCATGCCGgc                                                                                                                                                >  2:1778483/1‑150 (MQ=255)
          ggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGcc                                                                                                                                             <  2:4082154/150‑1 (MQ=255)
              aCTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGtt                                                                                                                                         >  2:1747238/1‑150 (MQ=255)
                  gCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCgg                                                                                                                                     >  1:113248/1‑150 (MQ=255)
                    tttacttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCcgcg                                                                                                                                   >  2:946012/7‑150 (MQ=255)
                    aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGTGCCGGTTTCGGCg                                                                                                                                   <  2:1834480/150‑1 (MQ=255)
                    aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                                                   >  2:4012083/1‑150 (MQ=255)
                    aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                                                   >  1:1399276/1‑150 (MQ=255)
                      gATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCACGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCAGTTTCGGCGAt                                                                                                                                 >  2:2558994/1‑150 (MQ=255)
                      gATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCAAGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGAt                                                                                                                                 >  1:960318/1‑150 (MQ=255)
                       aTCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGGCGGTGGCACCCATGTAGATCACGCAAGGCCGGCCGACGCGCGCGGCGACCGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATc                                                                                                                                >  2:298255/1‑150 (MQ=255)
                        tctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCa                                                                                                                               <  2:5380206/150‑1 (MQ=255)
                         ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                                                                              >  1:4534651/2‑147 (MQ=255)
                         ctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAg                                                                                                                              >  2:2614657/1‑150 (MQ=255)
                               gcCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGcc                                                                                                                                             >  1:2252265/1‑129 (MQ=255)
                               gcCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGcc                                                                                                                                             <  2:2252265/129‑1 (MQ=255)
                                cccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttt                                                                                                                       <  2:3178376/150‑1 (MQ=255)
                                 ccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                                                                      >  1:4068880/1‑150 (MQ=255)
                                  cagcagCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAGCGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGgtttt                                                                                                                     >  1:5388746/1‑150 (MQ=255)
                                  cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                                     <  2:127624/150‑1 (MQ=255)
                                  cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                                     <  2:4409882/150‑1 (MQ=255)
                                  cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                                     >  1:2648231/1‑150 (MQ=255)
                                  cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                                     <  1:2386149/150‑1 (MQ=255)
                                   agcagcTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTg                                                                                                                    <  1:2143037/150‑1 (MQ=255)
                                    gcagcTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGc                                                                                                                   <  2:2439389/150‑1 (MQ=255)
                                        cTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCAt                                                                                                               <  1:1094864/150‑1 (MQ=255)
                                            aTGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgcgc                                                                                                           >  1:3825565/1‑150 (MQ=255)
                                             tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCa                                                                                                                               <  1:688068/129‑1 (MQ=255)
                                             tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCt                                                                                                          >  2:5396128/1‑150 (MQ=255)
                                              gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                                                         <  2:3538905/150‑1 (MQ=255)
                                              gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                                                         <  2:624357/150‑1 (MQ=255)
                                                    cacGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcag                                                                                                   >  2:2337066/1‑150 (MQ=255)
                                                      cGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagca                                                                                                 <  1:194442/150‑1 (MQ=255)
                                                       gTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTg                                                                                                                    >  2:137408/1‑130 (MQ=255)
                                                       gTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagcag                                                                                                >  2:2912915/1‑150 (MQ=255)
                                                       gTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagcag                                                                                                >  2:1672327/1‑150 (MQ=255)
                                                          cGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCAc                                                                                             >  2:4014297/1‑150 (MQ=255)
                                                               gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                                                                        <  1:3986700/150‑1 (MQ=255)
                                                                ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                                                                       <  1:3383445/150‑1 (MQ=255)
                                                                   gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                                    >  2:1079777/1‑150 (MQ=255)
                                                                   gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                                    >  2:4247756/1‑150 (MQ=255)
                                                                      cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                 >  2:760183/1‑150 (MQ=255)
                                                                      cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                 <  1:2427726/150‑1 (MQ=255)
                                                                      cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                 >  2:5218795/1‑150 (MQ=255)
                                                                       tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                                                                                <  1:4645392/150‑1 (MQ=255)
                                                                             tcagTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgtta                                                                          <  2:3024314/147‑2 (MQ=255)
                                                                                gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                                                                       <  2:1686963/150‑1 (MQ=255)
                                                                                gTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgttaatc                                                                       >  1:3459595/1‑146 (MQ=255)
                                                                                 tGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCt                                                                      >  1:1312630/1‑150 (MQ=255)
                                                                                     aCCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtgtg                                                                  >  2:4890271/1‑150 (MQ=255)
                                                                                      cccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgc                                                                 <  1:1712426/150‑1 (MQ=255)
                                                                                       ccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTAGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcg                                                                >  2:671428/1‑150 (MQ=255)
                                                                                        cATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgc                                                               >  1:2767397/1‑150 (MQ=255)
                                                                                         aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                                              >  2:3321481/1‑150 (MQ=255)
                                                                                         aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                                              >  2:4769517/1‑150 (MQ=255)
                                                                                            tAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCg                                                           <  2:4448358/150‑1 (MQ=255)
                                                                                              gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                                                         <  2:3577549/150‑1 (MQ=255)
                                                                                              gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                                                         <  2:3651322/150‑1 (MQ=255)
                                                                                               aTCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTa                                                        >  1:1307909/1‑150 (MQ=255)
                                                                                                tCACGCAAGGCAGGCCGAAGCGCGCGGTGACAGTGGAGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCGTGTGCGCGCCGGTAt                                                       >  2:4875624/1‑150 (MQ=255)
                                                                                                 cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                             >  1:3342989/1‑127 (MQ=255)
                                                                                                 cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                             <  2:3342989/127‑1 (MQ=255)
                                                                                                 cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATg                                                      >  1:3760350/1‑150 (MQ=255)
                                                                                                 cACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATg                                                      <  1:4405528/150‑1 (MQ=255)
                                                                                                    gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                   >  1:4276336/1‑150 (MQ=255)
                                                                                                    gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                   <  1:3311439/150‑1 (MQ=255)
                                                                                                    gCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGtt                                                   >  1:2948899/1‑150 (MQ=255)
                                                                                                      aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTCGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa                                                 <  1:507703/150‑1 (MQ=255)
                                                                                                          caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                                             >  2:1522086/1‑150 (MQ=255)
                                                                                                          caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                                             <  2:3480785/150‑1 (MQ=255)
                                                                                                           aggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGctct                                            >  1:1372309/1‑150 (MQ=255)
                                                                                                            ggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCtt                                           <  1:1202865/150‑1 (MQ=255)
                                                                                                             gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCgg                                                          >  1:1262762/1‑134 (MQ=255)
                                                                                                             gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCgg                                                          <  2:1262762/134‑1 (MQ=255)
                                                                                                             gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTc                                          >  1:111192/1‑150 (MQ=255)
                                                                                                              ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                                                        >  2:880294/1‑87 (MQ=255)
                                                                                                                 aaGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGc                                      <  1:2707144/150‑1 (MQ=255)
                                                                                                                   gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                    <  2:5226135/150‑1 (MQ=255)
                                                                                                                   gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                    >  2:1655942/1‑150 (MQ=255)
                                                                                                                     gcgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTga                                  <  1:3313713/150‑1 (MQ=255)
                                                                                                                      cgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaa                                 <  2:3906596/150‑1 (MQ=255)
                                                                                                                      cgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaa                                 <  2:1020967/150‑1 (MQ=255)
                                                                                                                       gcggcgACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaag                                >  2:2189276/1‑150 (MQ=255)
                                                                                                                               atgcacggtgtcgCCACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                    >  1:5188199/6‑138 (MQ=255)
                                                                                                                               aGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCtt                        >  1:4525732/1‑150 (MQ=255)
                                                                                                                               aGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCtt                        >  1:4525777/1‑150 (MQ=255)
                                                                                                                                gtggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTc                       <  1:3957457/150‑1 (MQ=255)
                                                                                                                                 tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                                                        >  2:2161446/1‑68 (MQ=255)
                                                                                                                                 tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCg                      >  2:5346256/1‑150 (MQ=255)
                                                                                                                                   gcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCt                    >  1:594735/1‑150 (MQ=255)
                                                                                                                                   gcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCt                    >  2:4221908/1‑150 (MQ=255)
                                                                                                                                     ggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTcc                  >  2:198458/1‑150 (MQ=255)
                                                                                                                                        ggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTccgcc               >  2:3855112/1‑150 (MQ=255)
                                                                                                                                            acacCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCaa           <  2:113248/150‑1 (MQ=255)
                                                                                                                                            acacCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCaa           >  2:2615311/1‑150 (MQ=255)
                                                                                                                                                cGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCt       >  1:4220504/1‑150 (MQ=255)
                                                                                                                                                  tGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCg     <  1:489149/150‑1 (MQ=255)
                                                                                                                                                     aTGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTc  >  1:2607933/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTC  >  NC_002947/93080‑93378

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: