Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,358,513 Δ1 bp pseudogene (774/1285 nt) PA14_RS24540 → site‑specific integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,358,5100T.99.7% 661.6 / ‑3.5 301pseudogene (771/1285 nt)PA14_RS24540site‑specific integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (132/168);  total (132/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

GATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAT  >  NC_008463/5358361‑5358655
                                                                                                                                                     |                                                                                                                                                 
gATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGt                                                                                                                                                   >  1:7379651/1‑150 (MQ=32)
gATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGt                                                                                                                                                   >  2:5291182/1‑150 (MQ=32)
 aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtt                                                                                                                                                  >  1:10353381/1‑150 (MQ=32)
 aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtt                                                                                                                                                  <  2:349328/150‑1 (MQ=32)
  tAGAGCCGGATTCGCCGATGACGCCCAGTTGTTCAACGACAACCGGTTCTCACCGCACTACACGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGGGCGGATGCCCCCGCACCGttt                                                                                                                                                 <  2:548701/150‑1 (MQ=11)
  tAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                 >  2:8909152/1‑150 (MQ=32)
   agagCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttc                                                                                                                                                <  1:5159086/150‑2 (MQ=14)
    gagCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttcc                                                                                                                                               <  2:1906786/150‑3 (MQ=11)
     agCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttccg                                                                                                                                              >  2:10418429/1‑147 (MQ=1)
     agCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttccg                                                                                                                                              >  2:7453829/1‑147 (MQ=1)
                                                    cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTc                                                                                              >  1:3045538/1‑150 (MQ=1)
                                                    cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTc                                                                                              <  2:1433692/150‑1 (MQ=1)
                                                      ccGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                            >  2:3118501/1‑150 (MQ=11)
                                                      ccGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                            >  1:5782485/1‑150 (MQ=11)
                                                      ccGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                            >  1:7734225/1‑150 (MQ=11)
                                                       cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAcc                                                                                           >  2:2258180/1‑150 (MQ=11)
                                                       cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAcc                                                                                           >  2:7197423/1‑150 (MQ=11)
                                                        gCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCt                                                                                          <  1:10614878/150‑1 (MQ=11)
                                                         cACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTc                                                                                         >  1:4129681/1‑150 (MQ=11)
                                                         cACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTc                                                                                         <  1:2594123/150‑1 (MQ=11)
                                                          aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCa                                                                                        <  1:10027430/150‑1 (MQ=11)
                                                           cTACAAGAGCAAGGTGATGATCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCAc                                                                                       <  2:9049460/150‑1 (MQ=2)
                                                            tACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCCCACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACg                                                                                      >  2:8947935/1‑150 (MQ=2)
                                                              cAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCTAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGaa                                                                                    >  1:8083127/1‑150 (MQ=2)
                                                              cAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGaa                                                                                    >  2:8610234/1‑150 (MQ=11)
                                                               aaGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAg                                                                                   <  2:1080367/150‑1 (MQ=11)
                                                                agagCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGt                                                                                  <  1:8429447/150‑1 (MQ=11)
                                                                agagCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGt                                                                                  <  1:6358781/150‑1 (MQ=11)
                                                                 gagCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTg                                                                                 >  1:3258269/1‑150 (MQ=11)
                                                                   gCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGcc                                                                               >  1:4724773/1‑150 (MQ=11)
                                                                   gCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGcc                                                                               <  2:6711959/150‑1 (MQ=11)
                                                                       ggTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGTGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCGGGCCACCGACTGGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCt                                                                           <  2:416962/150‑1 (MQ=255)
                                                                         tgatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCa                                                                         >  1:5543524/1‑150 (MQ=11)
                                                                           atgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGACCATTCACCTCACGAAGTGCCTGCTCAAc                                                                       >  2:10209439/1‑150 (MQ=2)
                                                                            tgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAcc                                                                      >  2:2727543/1‑150 (MQ=11)
                                                                            tgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAcc                                                                      <  2:921238/150‑1 (MQ=11)
                                                                             gaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCa                                                                     <  2:8374459/150‑1 (MQ=11)
                                                                             gaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCa                                                                     >  1:4369100/1‑150 (MQ=11)
                                                                              aaCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCAc                                                                    <  2:4870442/150‑1 (MQ=11)
                                                                              aaCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCAc                                                                    >  2:6552608/1‑150 (MQ=11)
                                                                              aaCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCAc                                                                    <  1:3943488/150‑1 (MQ=11)
                                                                                  ccGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGa                                                                >  1:8576166/1‑150 (MQ=11)
                                                                                    gACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                              >  2:591620/1‑150 (MQ=11)
                                                                                     aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                             >  2:1478454/1‑150 (MQ=11)
                                                                                     aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                             >  1:6621854/1‑150 (MQ=11)
                                                                                     aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                             >  2:3137829/1‑150 (MQ=11)
                                                                                         ggTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCa                                                         >  1:4584173/1‑150 (MQ=11)
                                                                                         ggTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCa                                                         <  1:5312419/150‑1 (MQ=11)
                                                                                          gTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAg                                                        <  1:6912676/150‑1 (MQ=11)
                                                                                           tCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGa                                                       <  2:9623039/150‑1 (MQ=11)
                                                                                           tCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGa                                                       <  1:3572973/150‑1 (MQ=11)
                                                                                              aaGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCa                                                    <  2:3862379/150‑1 (MQ=11)
                                                                                              aaGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCa                                                    <  1:6716073/150‑1 (MQ=11)
                                                                                                gCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACg                                                  <  1:1487596/150‑1 (MQ=14)
                                                                                                  cATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGAt                                                <  1:2788438/150‑1 (MQ=17)
                                                                                                   aTGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATg                                               <  1:5239699/150‑1 (MQ=17)
                                                                                                   aTGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATg                                               <  1:5908412/150‑1 (MQ=17)
                                                                                                   aTGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATg                                               <  2:3619041/150‑1 (MQ=17)
                                                                                                     gTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAg                                             >  1:6031820/1‑150 (MQ=17)
                                                                                                      tACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGc                                            <  1:7032743/150‑1 (MQ=17)
                                                                                                      tACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGc                                            >  1:6192427/1‑150 (MQ=17)
                                                                                                       aCCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCt                                           <  2:1734700/150‑1 (MQ=17)
                                                                                                         cGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTAc                                         <  2:4456520/150‑1 (MQ=17)
                                                                                                              gTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGa                                    <  1:6240517/150‑1 (MQ=17)
                                                                                                               ttGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAg                                   <  2:8726056/150‑1 (MQ=17)
                                                                                                               ttGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAg                                   >  2:3951475/1‑150 (MQ=17)
                                                                                                               ttGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAg                                   <  1:3893669/150‑1 (MQ=17)
                                                                                                                tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAgg                                  <  2:5639490/150‑1 (MQ=17)
                                                                                                                 gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                 >  1:7475443/1‑150 (MQ=17)
                                                                                                                 gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                 >  1:10955039/1‑150 (MQ=17)
                                                                                                                 gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                 >  2:4999407/1‑150 (MQ=17)
                                                                                                                 gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                 <  1:9645328/150‑1 (MQ=17)
                                                                                                                  aCCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                >  2:3961833/1‑150 (MQ=17)
                                                                                                                   ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCg                               >  2:6986337/1‑150 (MQ=17)
                                                                                                                     gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGAc                             <  2:4107487/150‑1 (MQ=17)
                                                                                                                        aaGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCTACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTAc                          >  2:11096130/1‑150 (MQ=12)
                                                                                                                        aaGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTAc                          <  1:8837527/150‑1 (MQ=17)
                                                                                                                        aaGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTAc                          >  1:1304830/1‑150 (MQ=17)
                                                                                                                        aaGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCACCCCACCGACTTGATGCAGGCACCCGAGCGGAACATTCACCTCACCACGTGCCCGCACAACCACTCGAATCTCCAGACCCCGATGAGCTACATCGAGGCCGACTAc                          <  2:5091538/150‑1 (MQ=11)
                                                                                                                          ggTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGa                        <  1:2245750/150‑1 (MQ=17)
                                                                                                                           gTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGAt                       <  1:1886780/150‑1 (MQ=18)
                                                                                                                           gTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGAt                       >  1:595966/1‑150 (MQ=18)
                                                                                                                           gTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGAt                       >  1:7182001/1‑150 (MQ=18)
                                                                                                                            ttGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATc                      >  1:1737049/1‑150 (MQ=21)
                                                                                                                              ggggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcac                    <  2:2498061/150‑1 (MQ=25)
                                                                                                                               gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcaca                   <  2:3972725/150‑1 (MQ=25)
                                                                                                                                 gTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATg                 <  1:9529660/150‑1 (MQ=25)
                                                                                                                                  tGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGc                <  2:6229815/150‑1 (MQ=25)
                                                                                                                                   gCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcg               <  1:4878680/150‑1 (MQ=25)
                                                                                                                                    cGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgt              <  1:5053409/150‑1 (MQ=25)
                                                                                                                                    cGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgt              <  1:9322117/150‑1 (MQ=25)
                                                                                                                                       aTGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgcc           >  2:10338850/1‑150 (MQ=25)
                                                                                                                                         gACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgccgt         >  2:807390/1‑150 (MQ=25)
                                                                                                                                           ccccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGTGCTACATCGAGGCCGACTACGATCACATGcgtgccgtgc       <  1:1077626/150‑1 (MQ=19)
                                                                                                                                             ccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctg     >  2:9386968/1‑150 (MQ=17)
                                                                                                                                              cGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctgc    >  2:3060680/1‑150 (MQ=17)
                                                                                                                                               gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCa   >  1:6852195/1‑150 (MQ=17)
                                                                                                                                               gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCa   >  1:2373545/1‑150 (MQ=17)
                                                                                                                                               gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCa   <  2:3269464/150‑1 (MQ=17)
                                                                                                                                                cACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAt  <  1:6505269/150‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                                 
GATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAT  >  NC_008463/5358361‑5358655

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: