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% 594.7 / ‑3.4 290pseudogene (771/1285 nt)PA14_RS24540site‑specific integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (144/145);  total (144/146)
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 312 total aligned reads displayed.

CGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATG  >  NC_008463/5358360‑5358656
                                                                                                                                                      |                                                                                                                                                  
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGACCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  >  2:3822693/1‑151 (MQ=14)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  <  1:10005399/151‑1 (MQ=31)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  <  2:9235964/151‑1 (MQ=31)
     gagCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttccg                                                                                                                                               <  1:4184285/151‑4 (MQ=1)
                                                    tcACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTc                                                                                               >  2:10764559/1‑151 (MQ=1)
                                                     cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCa                                                                                              >  1:3671930/1‑151 (MQ=11)
                                                      aCCGCACTACAAGATCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                             <  1:3099404/151‑1 (MQ=2)
                                                      aCCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGGTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTACCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGACCATTCAc                                                                                             >  2:6773652/1‑151 (MQ=2)
                                                        cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCt                                                                                           >  1:8116170/1‑151 (MQ=11)
                                                        cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCt                                                                                           <  1:839265/151‑1 (MQ=11)
                                                        cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCt                                                                                           <  1:9450786/151‑1 (MQ=11)
                                                        cGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCt                                                                                           <  1:4979891/151‑1 (MQ=11)
                                                          cACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCGATGTACCGGATGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCa                                                                                         <  2:6738258/151‑1 (MQ=2)
                                                                acgagCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGt                                                                                   <  2:5557263/149‑1 (MQ=11)
                                                                aaGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGt                                                                                   >  1:601579/1‑151 (MQ=11)
                                                                  gagCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGTGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGc                                                                                 <  1:1243095/151‑1 (MQ=2)
                                                                   agCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGcc                                                                                <  1:5215077/151‑1 (MQ=11)
                                                                     cAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTg                                                                              >  1:1394593/1‑151 (MQ=11)
                                                                     cAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTg                                                                              <  2:5764755/151‑1 (MQ=11)
                                                                      aaGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGc                                                                             >  1:8962923/1‑151 (MQ=11)
                                                                         gTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCa                                                                          >  2:1523859/1‑151 (MQ=11)
                                                                         gTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCa                                                                          >  1:4439153/1‑151 (MQ=11)
                                                                           gatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAc                                                                        <  2:8465928/151‑1 (MQ=11)
                                                                           gatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAc                                                                        >  1:9815577/1‑151 (MQ=11)
                                                                               aaCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCCCt                                                                    >  2:7996324/1‑151 (MQ=2)
                                                                                aCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTc                                                                   <  1:3389892/151‑1 (MQ=11)
                                                                                aCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTc                                                                   <  2:271396/151‑1 (MQ=11)
                                                                                 cTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCg                                                                  <  2:4325856/151‑1 (MQ=11)
                                                                                    cGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                               >  2:1373134/1‑151 (MQ=11)
                                                                                    cGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                               <  1:9014183/151‑1 (MQ=11)
                                                                                    cGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                               <  1:25659/151‑1 (MQ=11)
                                                                                    cGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGAGGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTGGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                               >  2:1368801/1‑151 (MQ=2)
                                                                                     gACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                              >  1:3634388/1‑151 (MQ=11)
                                                                                     gACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                              >  1:8292085/1‑151 (MQ=11)
                                                                                        cAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATcc                                                           <  1:5842041/151‑1 (MQ=11)
                                                                                          ggTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAg                                                         <  1:10105207/151‑1 (MQ=11)
                                                                                            tCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAc                                                       <  2:8566381/151‑1 (MQ=11)
                                                                                             cGAAGCCATGTACCGGAAGTTGACCGAGCAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGTACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAcc                                                      >  2:5076473/1‑151 (MQ=2)
                                                                                             cGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAcc                                                      <  2:181874/151‑1 (MQ=11)
                                                                                             cGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAcc                                                      <  1:5470747/151‑1 (MQ=11)
                                                                                              gAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAccc                                                     >  2:5841965/1‑150 (MQ=11)
                                                                                              gAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCa                                                     >  1:1745693/1‑151 (MQ=11)
                                                                                               aaGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCAc                                                    <  1:3151171/151‑1 (MQ=11)
                                                                                                 gCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGACCATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGa                                                  >  1:4400133/1‑151 (MQ=9)
                                                                                                  ccATGTACCGGAAGTTGACCGAGAAGGGTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGTAGTGCCTGCTCAACCACTCGAATATCCAGACCACGAt                                                 >  2:8085669/1‑151 (MQ=12)
                                                                                                    aTGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGa                                               >  2:2867594/1‑151 (MQ=17)
                                                                                                       tACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCt                                            <  1:5663085/151‑1 (MQ=17)
                                                                                                         ccGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTAc                                          <  1:6427880/151‑1 (MQ=17)
                                                                                                           ggAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACAt                                        <  1:3803887/151‑1 (MQ=17)
                                                                                                            gAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATc                                       <  2:7527154/151‑1 (MQ=17)
                                                                                                            gAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACACGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATc                                       <  2:2585544/151‑1 (MQ=12)
                                                                                                             aaGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCg                                      >  1:8512399/1‑151 (MQ=17)
                                                                                                             aaGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCg                                      >  1:7099268/1‑151 (MQ=17)
                                                                                                              cgTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGCAGTGCCAGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGa                                     <  1:3370688/150‑1 (MQ=255)
                                                                                                              aGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACCTCGa                                     <  2:2359624/151‑1 (MQ=12)
                                                                                                              aGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGa                                     <  1:9916065/151‑1 (MQ=17)
                                                                                                              aGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGa                                     <  2:3912585/151‑1 (MQ=17)
                                                                                                              aGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGa                                     >  2:4353392/1‑151 (MQ=17)
                                                                                                               gTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAg                                    <  2:9696593/151‑1 (MQ=17)
                                                                                                                ttGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAgg                                   <  2:2939687/151‑1 (MQ=17)
                                                                                                                 tGACCGCGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGTTGAAGGCACCCGAGCGGACCATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  >  2:9097124/1‑151 (MQ=255)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  >  1:7461991/1‑151 (MQ=17)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  >  2:3446953/1‑151 (MQ=17)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  <  1:4676400/151‑1 (MQ=17)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  >  1:9081887/1‑151 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 >  2:8380326/1‑151 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 >  2:7715216/1‑151 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTCCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 >  2:4956889/1‑151 (MQ=12)
                                                                                                                   cccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCg                                <  2:10108554/150‑1 (MQ=17)
                                                                                                                    ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGa                               <  2:702651/151‑1 (MQ=17)
                                                                                                                    ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGa                               <  2:5880864/151‑1 (MQ=17)
                                                                                                                    ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGa                               >  1:6592264/1‑151 (MQ=17)
                                                                                                                    ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGa                               <  1:3242275/151‑1 (MQ=17)
                                                                                                                    ccGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGa                               <  1:10749103/151‑1 (MQ=17)
                                                                                                                     cGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCCCTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGAc                              <  2:6421972/151‑1 (MQ=12)
                                                                                                                      gagaCGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             <  2:6098212/151‑1 (MQ=12)
                                                                                                                      gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             >  1:2259778/1‑151 (MQ=17)
                                                                                                                      gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             >  1:7679045/1‑151 (MQ=17)
                                                                                                                      gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             >  1:10245919/1‑151 (MQ=17)
                                                                                                                       agaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTa                            >  2:3152972/1‑151 (MQ=17)
                                                                                                                           ggTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGAt                        >  1:6195744/1‑151 (MQ=18)
                                                                                                                           ggTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACAGTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGAt                        >  2:222/1‑151 (MQ=16)
                                                                                                                             ttGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATca                      <  1:7611913/151‑1 (MQ=21)
                                                                                                                              tGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCTCTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcac                     >  2:8829462/1‑151 (MQ=19)
                                                                                                                               ggggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcaca                    <  1:4394169/151‑1 (MQ=21)
                                                                                                                                gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACAt                   <  1:7511200/151‑1 (MQ=21)
                                                                                                                                gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACAt                   >  1:9538758/1‑151 (MQ=21)
                                                                                                                                gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACAt                   <  2:4617010/151‑1 (MQ=21)
                                                                                                                                 ggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATg                  <  2:6907515/151‑1 (MQ=21)
                                                                                                                                 ggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATg                  >  2:6532728/1‑151 (MQ=21)
                                                                                                                                 ggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATg                  >  1:3178410/1‑151 (MQ=21)
                                                                                                                                   tGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcg                <  2:4232508/151‑1 (MQ=21)
                                                                                                                                     cGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtg              >  1:6751698/1‑151 (MQ=21)
                                                                                                                                     cGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtg              >  1:7913606/1‑151 (MQ=21)
                                                                                                                                           aCCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgccgtgc        >  2:9583011/1‑151 (MQ=21)
                                                                                                                                            ccccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGACTATCCAGACCACGATGAGCTCCATCGAGGCCGACTACGATCACATGCGTGCCGtgct       >  2:10218294/1‑151 (MQ=255)
                                                                                                                                            ccccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgct       >  1:5096842/1‑151 (MQ=18)
                                                                                                                                             cccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctg      <  1:555766/151‑1 (MQ=18)
                                                                                                                                                gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAt   >  2:10504986/1‑151 (MQ=17)
                                                                                                                                                gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAt   <  2:2831999/151‑1 (MQ=17)
                                                                                                                                                      |                                                                                                                                                  
CGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATG  >  NC_008463/5358360‑5358656

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: