Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,851 A→T 5.1% intergenic (‑137/+43) aroQ‑I ← / ← dsbD‑I type II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,8510AT 5.1% 1139.5 / 30.5 333intergenic (‑137/+43)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base A (190/126);  new base T (9/8);  total (199/134)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 336 total aligned reads displayed.

GTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGCAA  >  NC_002947/650702‑651000
                                                                                                                                                     |                                                                                                                                                     
     gTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCg                                                                                                                                                  <  2:4592140/150‑1 (MQ=255)
       tGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGt                                                                                                                                                <  2:5523239/150‑1 (MQ=255)
        gCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTc                                                                                                                                               <  2:2261399/150‑1 (MQ=255)
                  gATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCa                                                                                                                                     <  1:3401625/150‑1 (MQ=255)
                       cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtt                                                                                                                                >  2:1686541/1‑150 (MQ=255)
                        ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAttt                                                                                                                               >  2:6202272/1‑150 (MQ=255)
                        ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAttt                                                                                                                               >  2:2627090/1‑150 (MQ=255)
                        ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAttt                                                                                                                               >  1:2103285/1‑150 (MQ=255)
                        ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAAACCGGTCACTCATTTCAGGAttt                                                                                                                               >  2:5093614/1‑150 (MQ=255)
                              ttGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGa                                                                                                                         <  1:6175246/150‑1 (MQ=255)
                                    ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcgc                                                                                                                   >  1:2570400/1‑150 (MQ=255)
                                    ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGAGATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcgc                                                                                                                   >  1:4006955/1‑150 (MQ=255)
                                     gCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcgcg                                                                                                                  >  1:4346394/1‑150 (MQ=255)
                                       aaTTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCggg                                                                                                                <  1:5005530/150‑1 (MQ=255)
                                         ttCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                              <  2:4160373/150‑1 (MQ=255)
                                         ttCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                              >  1:780718/1‑150 (MQ=255)
                                         ttCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                              <  1:1498242/150‑1 (MQ=255)
                                         ttCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                              >  1:966990/1‑150 (MQ=255)
                                            gACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                           >  1:5635941/1‑150 (MQ=255)
                                            gACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                           >  1:5922728/1‑150 (MQ=255)
                                               acaGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGAATGTGACCGGATATATAAATCCGGTCACTCATTACAGGATTTTTCCGAGCTCGCGGGAAGTtata                                                                                                        >  1:5094816/1‑150 (MQ=255)
                                                 aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                      >  2:5268267/1‑150 (MQ=255)
                                                 aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                      >  1:620745/1‑150 (MQ=255)
                                                 aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                      >  2:3049563/1‑150 (MQ=255)
                                                 aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                      >  1:5520194/1‑150 (MQ=255)
                                                    gcggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCg                                                                                                   >  1:4860150/1‑150 (MQ=255)
                                                      ggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGAc                                                                                                 >  2:657567/1‑150 (MQ=255)
                                                         gAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTGTAGACCGACttt                                                                                              >  2:3857327/1‑150 (MQ=255)
                                                         gAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACttt                                                                                              <  1:1058397/150‑1 (MQ=255)
                                                           gtgtGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGc                                                                                            >  1:1947214/1‑150 (MQ=255)
                                                            tgtgCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGcc                                                                                           >  1:493470/1‑150 (MQ=255)
                                                               gCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcgcg                                                                                        >  1:3579100/1‑150 (MQ=255)
                                                                ccGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCgcgc                                                                                       >  2:607800/1‑150 (MQ=255)
                                                                 cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCa                                                                                      >  2:1430515/1‑150 (MQ=255)
                                                                 cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCa                                                                                      >  2:5391156/1‑150 (MQ=255)
                                                                 cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCa                                                                                      >  2:5757592/1‑150 (MQ=255)
                                                                    aaaGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgc                                                                                   <  2:945438/150‑1 (MQ=255)
                                                                      aGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgc                                                                                 <  2:1526544/150‑1 (MQ=255)
                                                                      aGCAACGGCAAGTGTCCAGGTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgc                                                                                 <  1:3894016/150‑1 (MQ=255)
                                                                       gCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgcg                                                                                >  1:2292837/1‑150 (MQ=255)
                                                                        cAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGt                                                                               <  1:6674886/150‑1 (MQ=255)
                                                                         aaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTc                                                                              >  1:1688059/1‑150 (MQ=255)
                                                                          aCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCa                                                                             <  1:2936210/150‑1 (MQ=255)
                                                                           cGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAg                                                                            <  1:53268/150‑1 (MQ=255)
                                                                               aaGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATc                                                                        >  2:3757834/1‑150 (MQ=255)
                                                                               aaGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATc                                                                        >  2:254569/1‑150 (MQ=255)
                                                                                 gtgtCCAGTTCCCGGCAATAACCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCGCTCCTTTCAGGATTTTTCCGAGCTCGCGGGAAGGTATAGACCGACATTGCCGCGCACGCGCGGCAGAATCTg                                                                      >  2:2155271/1‑150 (MQ=255)
                                                                                  tgtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                     <  2:3469366/150‑1 (MQ=255)
                                                                                  tgtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                     <  1:3459193/150‑1 (MQ=255)
                                                                                  tgtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                     <  2:4283337/150‑1 (MQ=255)
                                                                                    tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCCCTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                                   >  1:5767501/1‑150 (MQ=255)
                                                                                    tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                                   >  2:2315286/1‑150 (MQ=255)
                                                                                    tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                                   >  2:3107515/1‑150 (MQ=255)
                                                                                    tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                                   >  2:4095915/1‑150 (MQ=255)
                                                                                      cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGa                                                                 >  1:4743045/1‑150 (MQ=255)
                                                                                      cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGa                                                                 >  2:2854650/1‑150 (MQ=255)
                                                                                      cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGa                                                                 >  2:5127802/1‑150 (MQ=255)
                                                                                       aGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGaa                                                                >  1:72419/1‑150 (MQ=255)
                                                                                        gTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAt                                                               >  1:4494890/1‑150 (MQ=255)
                                                                                         ttCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                              >  1:2768908/1‑150 (MQ=255)
                                                                                         ttCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                              >  2:4672412/1‑150 (MQ=255)
                                                                                            ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                           <  1:5606277/150‑1 (MQ=255)
                                                                                            ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                           >  2:197891/1‑150 (MQ=255)
                                                                                            ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                           >  1:3148505/1‑150 (MQ=255)
                                                                                            ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                           >  2:2784995/1‑150 (MQ=255)
                                                                                             cGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAcc                                                          >  1:2804389/1‑150 (MQ=255)
                                                                                               gCAAGAGCCAGCACGAGGTCCGAGGGTCGCGCGAAGTCTGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTGCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCgg                                                        <  2:5902253/150‑1 (MQ=255)
                                                                                                cAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGc                                                       <  2:1077035/150‑1 (MQ=255)
                                                                                                 aaTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCa                                                      <  1:5665338/150‑1 (MQ=255)
                                                                                                 aaTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCa                                                      <  1:3638909/150‑1 (MQ=255)
                                                                                                  aTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAt                                                     >  1:380511/1‑150 (MQ=255)
                                                                                                      ccAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAt                                                 >  1:3050203/1‑150 (MQ=255)
                                                                                                      ccAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAt                                                 >  1:5375922/1‑150 (MQ=255)
                                                                                                       cAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtc                                                >  2:834517/1‑150 (MQ=255)
                                                                                                         gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                              <  1:5274917/150‑1 (MQ=255)
                                                                                                         gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                              >  1:5830692/1‑150 (MQ=255)
                                                                                                         gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                              >  1:2207434/1‑150 (MQ=255)
                                                                                                           aCGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGc                                            <  1:1382631/150‑1 (MQ=255)
                                                                                                           aCGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGc                                            >  1:6758510/1‑150 (MQ=255)
                                                                                                                    gAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAccc                                   >  2:6348782/1‑150 (MQ=255)
                                                                                                                        ttCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATc                               >  1:414823/1‑150 (MQ=255)
                                                                                                                                aaGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                       >  2:1254223/1‑150 (MQ=255)
                                                                                                                                aaGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                       >  1:4270153/1‑150 (MQ=255)
                                                                                                                                  atgagtgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAAtt                     >  1:4643446/5‑150 (MQ=255)
                                                                                                                                  atgagtgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAAtt                     >  1:2435862/5‑150 (MQ=255)
                                                                                                                                   tATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAAttt                    <  2:4743262/150‑1 (MQ=255)
                                                                                                                                      gtgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc                 >  1:456540/1‑150 (MQ=255)
                                                                                                                                       tgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCt                >  2:1842516/1‑150 (MQ=255)
                                                                                                                                       tgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCt                >  1:1313355/1‑150 (MQ=255)
                                                                                                                                        gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTg               >  1:2004481/1‑150 (MQ=255)
                                                                                                                                         aCCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGc              >  2:2688949/1‑150 (MQ=255)
                                                                                                                                           cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGACTTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCg            >  2:6443074/1‑150 (MQ=255)
                                                                                                                                           cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCg            >  1:1597870/1‑150 (MQ=255)
                                                                                                                                            ggATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGt           >  2:3526150/1‑150 (MQ=255)
                                                                                                                                             gATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGtt          >  1:3327108/1‑150 (MQ=255)
                                                                                                                                                atTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgc       >  2:2379985/1‑150 (MQ=255)
                                                                                                                                                 tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca      <  1:2654546/150‑1 (MQ=255)
                                                                                                                                                 tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca      <  2:5612233/150‑1 (MQ=255)
                                                                                                                                                 tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca      <  1:6508543/150‑1 (MQ=255)
                                                                                                                                                 tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca      <  1:1173418/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGCAA  >  NC_002947/650702‑651000

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: