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.3% 580.7 / 0.4 297pseudogene (771/1285 nt)PA14_RS24540site‑specific integrase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (153/141);  total (155/142)
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 318 total aligned reads displayed.

CGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATG  >  NC_008463/5358360‑5358656
                                                                                                                                                      |                                                                                                                                                  
cGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGt                                                                                                                                                    >  1:5169617/1‑151 (MQ=31)
cGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGt                                                                                                                                                    >  2:1071497/1‑151 (MQ=31)
 gATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtt                                                                                                                                                   >  2:10568850/1‑151 (MQ=31)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  >  2:6557777/1‑151 (MQ=31)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  <  1:9558789/151‑1 (MQ=31)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  >  1:7892718/1‑151 (MQ=31)
  aTAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGttt                                                                                                                                                  <  1:1784397/151‑1 (MQ=31)
   tAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttc                                                                                                                                                 >  2:9435182/1‑150 (MQ=14)
   tAGAGCCGGATTCGCCGATGACGAACAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGtttc                                                                                                                                                 <  2:822687/151‑2 (MQ=14)
                                                    tcACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTc                                                                                               >  2:3817531/1‑151 (MQ=1)
                                                     cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCa                                                                                              >  1:2521839/1‑151 (MQ=11)
                                                     cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCa                                                                                              <  2:2310772/151‑1 (MQ=11)
                                                      aCCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                             >  2:8791007/1‑151 (MQ=11)
                                                      aCCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCAc                                                                                             <  2:6374820/151‑1 (MQ=11)
                                                          cACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCa                                                                                         <  1:1210516/151‑1 (MQ=11)
                                                          cACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCa                                                                                         <  2:10988007/151‑1 (MQ=11)
                                                            cTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACg                                                                                       <  2:4804179/151‑1 (MQ=11)
                                                               cAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAg                                                                                    >  1:6246191/1‑151 (MQ=11)
                                                                aaGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGt                                                                                   <  2:1141114/151‑1 (MQ=11)
                                                                    gCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCt                                                                               <  1:4671127/151‑1 (MQ=11)
                                                                    gCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCt                                                                               <  2:1342659/151‑1 (MQ=11)
                                                                      aaGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGc                                                                             >  1:1920387/1‑151 (MQ=11)
                                                                       aGGTGATGATCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCt                                                                            <  2:10307488/151‑1 (MQ=2)
                                                                       aGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCt                                                                            >  2:4343169/1‑151 (MQ=11)
                                                                       aGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCt                                                                            >  1:1135885/1‑151 (MQ=11)
                                                                        ggTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGTCCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGATCATTCACCTCACGAAGTGCCTGCTc                                                                           >  2:819959/1‑151 (MQ=2)
                                                                           gatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAc                                                                        <  1:5702979/151‑1 (MQ=11)
                                                                           gatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAc                                                                        <  1:9953723/151‑1 (MQ=11)
                                                                            atgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAcc                                                                       >  2:7314902/1‑151 (MQ=11)
                                                                            atgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAAcc                                                                       >  1:8076500/1‑151 (MQ=11)
                                                                               aaCTCCGACCAGGTCGAAGCCATGTACCGAAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACt                                                                    <  1:1310835/151‑1 (MQ=2)
                                                                                aCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACATCACGACGTGCCTGCTCAACCACTc                                                                   <  2:2671363/151‑1 (MQ=2)
                                                                                 cTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCg                                                                  <  2:3619181/151‑1 (MQ=11)
                                                                                   ccGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGaa                                                                >  2:3702926/1‑151 (MQ=11)
                                                                                    cGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAat                                                               >  2:55367/1‑151 (MQ=11)
                                                                                     gACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                              <  2:6038566/151‑1 (MQ=11)
                                                                                      aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGACGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAatat                                                             >  2:7708346/1‑151 (MQ=2)
                                                                                      aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAatat                                                             <  2:8539835/151‑1 (MQ=11)
                                                                                        cAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATcc                                                           >  2:156122/1‑151 (MQ=11)
                                                                                        cAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATcc                                                           <  2:2119814/151‑1 (MQ=11)
                                                                                        cAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATcc                                                           <  1:8461794/151‑1 (MQ=11)
                                                                                           gTCGAAGCCATGTTCCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAAAATCCAGa                                                        <  2:7907400/151‑1 (MQ=2)
                                                                                           gTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGCCTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGa                                                        >  2:669494/1‑151 (MQ=2)
                                                                                           gTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGa                                                        <  2:6822337/151‑1 (MQ=11)
                                                                                           gTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGa                                                        >  1:7018599/1‑151 (MQ=11)
                                                                                            tCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAc                                                       <  1:4171275/151‑1 (MQ=11)
                                                                                              gAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCa                                                     >  1:9789334/1‑151 (MQ=11)
                                                                                              gAAGCCATGGACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGTCCCCCCACCGGTTACGGCACACCCGGGCCCCCGACGTGATGAAGGCCCCCGCGCGGCACCGTCACCACACGAAGTGCCTGCGCACCCCCTCGAATATCCCGAccc                                                     >  2:8844385/1‑150 (MQ=2)
                                                                                               aaGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCAc                                                    >  1:6629663/1‑151 (MQ=11)
                                                                                                aGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACg                                                   >  2:9263887/1‑151 (MQ=14)
                                                                                                 gCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCCCCGACTGGAGGAAGGCACCCGCGCGGAACCGTCACCTCCCGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACgt                                                  >  2:6683511/1‑150 (MQ=11)
                                                                                                  ccATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGAt                                                 <  1:6351700/151‑1 (MQ=17)
                                                                                                   cATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATg                                                >  1:1574532/1‑151 (MQ=17)
                                                                                                    aTGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGa                                               <  1:8447136/151‑1 (MQ=17)
                                                                                                     tGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAg                                              >  2:7057447/1‑151 (MQ=17)
                                                                                                      gTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGc                                             <  1:2848826/151‑1 (MQ=17)
                                                                                                           ggAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACAt                                        <  1:10413256/151‑1 (MQ=17)
                                                                                                             aaGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCg                                      <  1:2237278/151‑1 (MQ=17)
                                                                                                              aGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTAGAATATCCAGACCACGATGAGCTACATCGa                                     <  1:10796440/151‑1 (MQ=12)
                                                                                                               gTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAg                                    <  2:5639048/151‑1 (MQ=17)
                                                                                                                ttGACCGAGATGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAgg                                   <  1:7594962/151‑1 (MQ=12)
                                                                                                                ttGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAgg                                   <  1:9127364/151‑1 (MQ=17)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  >  1:2004081/1‑151 (MQ=17)
                                                                                                                 tGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGc                                  <  2:4463728/151‑1 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 <  1:10487653/151‑1 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 <  2:5787972/151‑1 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGcc                                 <  1:10096895/151‑1 (MQ=17)
                                                                                                                  gACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACAGCGAGGcc                                 >  2:10195998/1‑151 (MQ=12)
                                                                                                                     cGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGAc                              <  1:5494515/151‑1 (MQ=17)
                                                                                                                      gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             >  1:2779313/1‑151 (MQ=17)
                                                                                                                      gagaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACt                             <  2:9109715/151‑1 (MQ=17)
                                                                                                                       agaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCTACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTa                            <  2:7829984/151‑1 (MQ=12)
                                                                                                                       agaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTa                            >  2:8960921/1‑151 (MQ=17)
                                                                                                                        gaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTAc                           <  2:3169742/151‑1 (MQ=17)
                                                                                                                        gaAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTAc                           <  2:4191634/151‑1 (MQ=17)
                                                                                                                         aaGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACg                          <  1:7500631/151‑1 (MQ=17)
                                                                                                                          aGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCCCTCGAAGATCCAGACCACGATGAGCTACATCGAGGCCGACTACGa                         >  2:1418275/1‑151 (MQ=12)
                                                                                                                              tGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcac                     >  1:393328/1‑151 (MQ=21)
                                                                                                                               ggggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATcaca                    >  1:4280766/1‑151 (MQ=21)
                                                                                                                                gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACAt                   <  1:2152274/151‑1 (MQ=21)
                                                                                                                                gggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACAt                   >  1:3620957/1‑151 (MQ=21)
                                                                                                                                 ggTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATg                  >  2:10633969/1‑151 (MQ=21)
                                                                                                                                  gTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGc                 >  1:5129588/1‑151 (MQ=21)
                                                                                                                                    gCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgt               >  1:3039508/1‑151 (MQ=21)
                                                                                                                                    gCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgt               <  2:3201457/151‑1 (MQ=21)
                                                                                                                                         tGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgccgt          <  1:693576/151‑1 (MQ=21)
                                                                                                                                          gACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgccgtg         <  2:2884837/151‑1 (MQ=21)
                                                                                                                                           cccccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGACGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGACTATCCAGACCACCATGAGCTACATCGAGGCCGACTACGATCACCTGcgtgccgtgc        <  2:8320956/150‑1 (MQ=1)
                                                                                                                                           aCCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGcgtgccgtgc        <  1:3760837/151‑1 (MQ=21)
                                                                                                                                             cccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctg      <  1:262518/151‑1 (MQ=18)
                                                                                                                                             cccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctg      >  2:1928378/1‑151 (MQ=18)
                                                                                                                                              ccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctgc     >  2:5214237/1‑151 (MQ=17)
                                                                                                                                              ccGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGtgctgc     <  1:8658300/151‑1 (MQ=17)
                                                                                                                                               cGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCa    <  1:420959/151‑1 (MQ=17)
                                                                                                                                               cGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCACACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCa    >  1:8782325/1‑151 (MQ=12)
                                                                                                                                                gCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCAt   >  1:830817/1‑151 (MQ=17)
                                                                                                                                                 cACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATg  >  2:9525262/1‑151 (MQ=17)
                                                                                                                                                 cACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATg  >  2:739863/1‑151 (MQ=17)
                                                                                                                                                 cACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATg  >  2:3954338/1‑151 (MQ=17)
                                                                                                                                                 cACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATTCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACGATGAGCTACATCGAGGCCGACTACGATCACATGCGTGCCGTGCTGCATg  >  2:492276/1‑151 (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: