Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,024 T→G 10.3% E190D (GAA→GAC trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,0240TG10.3% 615.7 / 49.8 203E190D (GAA→GACtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (87/95);  new base G (7/14);  total (94/109)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 205 total aligned reads displayed.

CGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATG  >  NC_002947/92876‑93171
                                                                                                                                                    |                                                                                                                                                   
cgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCtt                                                                                                                                                    >  1:2722407/1‑150 (MQ=255)
 gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTc                                                                                                                                                   >  2:5005577/1‑150 (MQ=255)
    gCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGAc                                                                                                                                                >  1:3818923/1‑150 (MQ=255)
      aGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGt                                                                                                                                              >  1:900795/1‑150 (MQ=255)
       ggCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGtt                                                                                                                                             >  2:3598824/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)
              aGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAcc                                                                                                                                      <  1:2768716/150‑1 (MQ=255)
                gCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCa                                                                                                                                    <  2:2618623/150‑1 (MQ=255)
                gCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCa                                                                                                                                    <  1:5264550/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)
                      tGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACg                                                                                                                              >  2:5259270/1‑150 (MQ=255)
                       ggCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGTAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGc                                                                                                                             >  2:4917071/1‑150 (MQ=255)
                                 cttccttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGt                                                                                                                   <  1:5065111/150‑1 (MQ=255)
                                   tccttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTc                                                                                                                 <  1:3938222/150‑1 (MQ=255)
                                     cttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGTGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                                                                               <  1:5376182/150‑1 (MQ=255)
                                     cttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                                                                               <  1:4087102/150‑1 (MQ=255)
                                      ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCtt                                                                                                                      <  2:2941778/142‑1 (MQ=255)
                                      ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCtt                                                                                                                      >  1:2941778/1‑142 (MQ=255)
                                      ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATg                                                                                                              >  2:2338592/1‑150 (MQ=255)
                                      ttcTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATg                                                                                                              >  1:4815417/1‑150 (MQ=255)
                                       tctcTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATgg                                                                                                             <  2:5023317/150‑1 (MQ=255)
                                           ttGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTc                                                                                                         >  2:1601233/1‑150 (MQ=255)
                                           ttGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTc                                                                                                         >  1:3257484/1‑150 (MQ=255)
                                           ttGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTc                                                                                                         >  2:2247505/1‑150 (MQ=255)
                                             gCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCtt                                                                                                       >  1:1674948/1‑150 (MQ=255)
                                             gCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCtt                                                                                                       >  1:1834480/1‑150 (MQ=255)
                                             gCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCtt                                                                                                       >  2:4480881/1‑150 (MQ=255)
                                               aaCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCAc                                                                                                                               <  1:4002166/124‑1 (MQ=255)
                                               aaCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTc                                                                                                     >  1:1048863/1‑150 (MQ=255)
                                                aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                                                                    <  1:2329811/150‑1 (MQ=255)
                                                aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                                                                    <  1:696908/150‑1 (MQ=255)
                                                 cTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                                                                                                                                            <  1:2764044/109‑1 (MQ=255)
                                                         gggTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCgg                                                                                           >  1:3178376/1‑150 (MQ=255)
                                                             ttcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcc                                                                                       >  1:3273329/6‑150 (MQ=255)
                                                              tcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcca                                                                                      <  1:4266510/146‑2 (MQ=255)
                                                              tcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcca                                                                                      <  2:146561/146‑2 (MQ=255)
                                                              tcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcca                                                                                      <  2:4666654/146‑2 (MQ=255)
                                                               ctTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcgg                                                                                     >  1:127624/1‑150 (MQ=255)
                                                                tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                                                                               <  1:4047890/123‑1 (MQ=255)
                                                                tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAt                                                                                                               >  2:4047890/1‑123 (MQ=255)
                                                                tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                                                                    >  1:1570/1‑150 (MQ=255)
                                                                tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                                                                    >  1:624357/1‑150 (MQ=255)
                                                                 ttGCCGAGGACCGCCTTGACGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGCTCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACCCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggct                                                                                   <  2:3297721/150‑2 (MQ=255)
                                                                 gtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCcagcc                                                                                   <  1:2535909/149‑2 (MQ=255)
                                                                  tGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcgg                                                                                  <  2:2666445/150‑1 (MQ=255)
                                                                    ccGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGtt                                                                                >  1:2211317/1‑149 (MQ=255)
                                                                      gatcaccGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTAc                                                                              <  2:3767896/144‑1 (MQ=255)
                                                                           tCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCa                                                                         <  1:1090417/150‑1 (MQ=255)
                                                                           tCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGACCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCa                                                                         <  2:4363513/150‑1 (MQ=255)
                                                                            cGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGctt                                                                        >  2:1341355/1‑148 (MQ=255)
                                                                              aCTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTACCTGGCACGa                                                                      >  2:2886617/1‑150 (MQ=255)
                                                                              aCTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGa                                                                      >  1:1357728/1‑150 (MQ=255)
                                                                               cTGGAAGTCGCGCACCATGGCGGGGTACTGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCAAGGCGTCTTTCAGGCTACCGGTGCCGGCGGTAACTGGCACGAt                                                                     >  2:2143037/1‑150 (MQ=255)
                                                                               cTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAt                                                                     >  2:3849657/1‑150 (MQ=255)
                                                                                tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAgg                                                                                                  >  1:2643880/1‑120 (MQ=255)
                                                                                   aaGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCg                                                                 <  2:5233575/150‑1 (MQ=255)
                                                                                   aaGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTATTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCg                                                                 <  1:768565/150‑1 (MQ=255)
                                                                                    aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg                                                                <  2:2969802/150‑1 (MQ=255)
                                                                                    aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg                                                                >  2:17241/1‑150 (MQ=255)
                                                                                     gTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGgc                                                               <  2:3027670/150‑1 (MQ=255)
                                                                                        gcgcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGcc                                                            >  1:4409882/1‑150 (MQ=255)
                                                                                         cgcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGccc                                                           >  1:4315975/1‑150 (MQ=255)
                                                                                                  gCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTc                                                  >  1:2065135/1‑150 (MQ=255)
                                                                                                            ggTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                   >  2:2738747/1‑107 (MQ=255)
                                                                                                            ggTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAaca                                        <  2:2545376/150‑1 (MQ=255)
                                                                                                             gtgtGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAacac                                       >  2:2427726/1‑150 (MQ=255)
                                                                                                               gtgGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCCTCATGCGGAACACGt                                     <  1:764022/150‑1 (MQ=255)
                                                                                                                 gggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTc                                   <  1:349754/150‑1 (MQ=255)
                                                                                                                  ggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCg                                  >  2:2386149/1‑150 (MQ=255)
                                                                                                                    ccagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGcc                                >  2:3999072/1‑150 (MQ=255)
                                                                                                                     cggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGtttggct                               >  2:66133/3‑145 (MQ=255)
                                                                                                                     cagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCct                               >  1:3577549/1‑150 (MQ=255)
                                                                                                                       gccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCctgc                             <  1:4917071/150‑1 (MQ=255)
                                                                                                                         caCTGTGCCGATCAGGTAGTAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCctgctg                           <  1:2379088/150‑1 (MQ=255)
                                                                                                                          aCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCa                                                                         <  1:2514133/103‑1 (MQ=255)
                                                                                                                            ggTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGgcg                        <  2:1075183/150‑1 (MQ=255)
                                                                                                                             gTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGGTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGacgt                       <  2:4575657/150‑2 (MQ=255)
                                                                                                                               gCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTc                     <  2:1959006/150‑1 (MQ=255)
                                                                                                                                  gATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGAt                  >  2:194442/1‑150 (MQ=255)
                                                                                                                                  gATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGAt                  >  1:3097639/1‑150 (MQ=255)
                                                                                                                                   aTCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATg                 >  2:694563/1‑150 (MQ=255)
                                                                                                                                   aTCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATg                 >  1:5152676/1‑150 (MQ=255)
                                                                                                                                    tCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGt                >  1:4448358/1‑150 (MQ=255)
                                                                                                                                    tCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGt                <  1:5005589/150‑1 (MQ=255)
                                                                                                                                      aGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCg              <  2:4777603/150‑1 (MQ=255)
                                                                                                                                      aGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCg              <  1:669017/150‑1 (MQ=255)
                                                                                                                                        gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt            >  2:3313713/1‑150 (MQ=255)
                                                                                                                                        gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt            <  2:1372439/150‑1 (MQ=255)
                                                                                                                                        gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt            >  1:1126703/1‑150 (MQ=255)
                                                                                                                                        gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt            <  2:2722407/150‑1 (MQ=255)
                                                                                                                                        gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGGACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt            <  2:2922013/150‑1 (MQ=255)
                                                                                                                                          aGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTgg          <  1:3430464/150‑1 (MQ=255)
                                                                                                                                           gAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGc         >  2:5318542/1‑150 (MQ=255)
                                                                                                                                           ataGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTc                     <  1:5382881/136‑1 (MQ=255)
                                                                                                                                            aaGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCa        >  1:4074267/1‑150 (MQ=255)
                                                                                                                                                tgtCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCa    <  1:291837/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATG  >  NC_002947/92876‑93171

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: