Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 6,070,122 T→C pseudogene (591/1086 nt) PA14_RS27735 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,070,1220TC99.8% 2538.9 / NA 643pseudogene (591/1086 nt)PA14_RS27735oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (312/330);  minor base A (1/0);  total (313/330)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 643 total aligned reads displayed.

CAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGGCGCCACGCTTC  >  NC_008463/6069972‑6070269
                                                                                                                                                      |                                                                                                                                                   
 agcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCa                                                                                                                                                 <  2:833478/151‑1 (MQ=255)
  gcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCAt                                                                                                                                                <  2:8041019/151‑1 (MQ=255)
   caTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATa                                                                                                                                               <  2:3226103/151‑1 (MQ=255)
    aTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATAc                                                                                                                                              >  2:1901927/1‑151 (MQ=255)
        gCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCa                                                                                                                                          >  1:7044716/1‑151 (MQ=255)
        gCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCa                                                                                                                                          >  2:3456395/1‑151 (MQ=255)
         cccGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAg                                                                                                                                         <  2:3993111/151‑1 (MQ=255)
         cccGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAg                                                                                                                                         <  2:9073913/151‑1 (MQ=255)
         cccGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAg                                                                                                                                         >  1:7800706/1‑151 (MQ=255)
             gTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTa                                                                                                                                     <  1:5974390/151‑1 (MQ=255)
             gTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTa                                                                                                                                     <  1:2221597/151‑1 (MQ=255)
             gTGCGCTACTTCGAGTCGGCCATCGAGCGCGTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTa                                                                                                                                     <  2:10713018/151‑1 (MQ=255)
                     cTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGa                                                                                                                             >  2:7462517/1‑151 (MQ=255)
                      ttCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGac                                                                                                                            >  1:10560562/1‑151 (MQ=255)
                      ttCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGac                                                                                                                            <  1:7455446/151‑1 (MQ=255)
                            tCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                                                                      <  2:6139770/151‑1 (MQ=255)
                            tCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                                                                      <  2:3717569/151‑1 (MQ=255)
                            tCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                                                                      <  1:1986353/151‑1 (MQ=255)
                            tCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                                                                      <  1:1777369/151‑1 (MQ=255)
                             cGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgcg                                                                                                                     <  2:1448086/151‑1 (MQ=255)
                              ggCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgcgc                                                                                                                    >  1:4041767/1‑151 (MQ=255)
                               gCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGcc                                                                                                                   <  1:4424151/151‑1 (MQ=255)
                                 cATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCgg                                                                                                                 <  2:2046084/151‑1 (MQ=255)
                                    cGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTg                                                                                                              <  2:7144281/151‑1 (MQ=255)
                                     gAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgc                                                                                                             <  2:5277795/151‑1 (MQ=255)
                                     gAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgc                                                                                                             <  2:8877355/151‑1 (MQ=255)
                                      aGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgcg                                                                                                            <  1:5329629/151‑1 (MQ=255)
                                      aGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgcg                                                                                                            >  1:1104494/1‑151 (MQ=255)
                                        cgcTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGcgcg                                                                                                          >  2:5903639/1‑151 (MQ=255)
                                        cgcTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGcgcg                                                                                                          <  2:6514302/151‑1 (MQ=255)
                                          cTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAg                                                                                                        <  1:2593726/151‑1 (MQ=255)
                                          cTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAg                                                                                                        <  2:9649267/151‑1 (MQ=255)
                                            tccgtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAggg                                                                                                      <  1:635356/151‑1 (MQ=255)
                                              cgtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcg                                                                                                    <  2:2771883/151‑1 (MQ=255)
                                               gtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                                                                   >  1:2776709/1‑151 (MQ=255)
                                                 ccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcgc                                                                                                 >  1:4925395/1‑151 (MQ=255)
                                                  cgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCt                                                                                                >  1:7320915/1‑151 (MQ=255)
                                                    cAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGa                                                                                              >  2:9813672/1‑151 (MQ=255)
                                                     aGGTCCGCGACCGCTGGCTCGAACACGACCTGCACGTC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACAGCAGCGTATGCACGTACACCTGCGCCGGCTGCGCGAGGGCGCGCAGAc                                                                                             >  2:1086358/1‑151 (MQ=255)
                                                       gTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcg                                                                                           <  1:687183/151‑1 (MQ=255)
                                                         ccGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcgac                                                                                         <  1:6053546/151‑1 (MQ=255)
                                                          cgcgACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcgacg                                                                                        >  2:1691939/1‑151 (MQ=255)
                                                          cgcgACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcgacg                                                                                        >  1:1735624/1‑151 (MQ=255)
                                                            cgACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGTCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCATGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCgacgac                                                                                      >  1:8436887/1‑151 (MQ=255)
                                                            cgACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCgacgac                                                                                      <  2:2557443/151‑1 (MQ=255)
                                                             gACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACt                                                                                     <  1:8378535/151‑1 (MQ=255)
                                                                cGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGt                                                                                  <  1:10467213/151‑1 (MQ=255)
                                                                 gCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGtt                                                                                 <  2:442858/151‑1 (MQ=255)
                                                                   tGGCGCGAACACGACCTGCCCGGC‑TCCGGGCAGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCTCTGTGCCTGTTCGGCAGACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCg                                                                               <  2:4664705/151‑1 (MQ=255)
                                                                    ggCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGa                                                                              <  1:6188196/151‑1 (MQ=255)
                                                                     gcgcGAACACGCCCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGTCC‑AGGCACTGGGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAg                                                                             <  2:6674479/151‑1 (MQ=255)
                                                                      cgcgAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCTCTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTACGAgc                                                                            <  1:306085/151‑2 (MQ=255)
                                                                      cgcgAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAgg                                                                            >  2:11163297/1‑151 (MQ=255)
                                                                       gcgAACACGACCTGCCCGGC‑TCCGGGCGGTGGGTCGTGCTCGGCCCGC‑CCCTGGTCGACC‑AGGCACTGGGCCTGGTCGGCTTACCGCGGCGGATGCACGGACACCTGCGGCGGCTGCGCGAGGGCGCGCTGGCCGTCGTGTGGGTCGTggg                                                                           >  2:2876261/1‑150 (MQ=255)
                                                                        cgAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTc                                                                          <  2:8594079/151‑1 (MQ=255)
                                                                         gAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCCGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTcc                                                                         >  1:4273326/1‑151 (MQ=255)
                                                                         gAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTcc                                                                         >  1:8469494/1‑151 (MQ=255)
                                                                           acacGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTg                                                                       <  2:5645911/151‑1 (MQ=255)
                                                                               gaacTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACCCGGGC‑CTGTGCCTGGTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACGGGTGCGAGGTCCTGcgcg                                                                   >  2:7686525/4‑151 (MQ=255)
                                                                               gACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCg                                                                   >  2:583990/1‑151 (MQ=255)
                                                                               gACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCg                                                                   >  1:10289174/1‑151 (MQ=255)
                                                                               gACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCg                                                                   <  1:1156681/151‑1 (MQ=255)
                                                                                 ccTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGAc                                                                 >  1:8231363/1‑151 (MQ=255)
                                                                                  cTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACt                                                                <  1:10016792/151‑1 (MQ=255)
                                                                                  cTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACt                                                                >  2:6259672/1‑151 (MQ=255)
                                                                                         ggctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGc                                                          >  1:9318938/4‑151 (MQ=255)
                                                                                          gctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCa                                                         <  2:10948136/149‑1 (MQ=255)
                                                                                           ctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCAc                                                        >  1:9604707/2‑151 (MQ=255)
                                                                                           ctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCAc                                                        <  1:10593529/150‑1 (MQ=255)
                                                                                              cGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGc                                                     >  1:9849443/1‑151 (MQ=255)
                                                                                               gggCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGcc                                                    >  2:9322643/1‑151 (MQ=255)
                                                                                               gggCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGcc                                                    >  1:10753709/1‑151 (MQ=255)
                                                                                                    gtgGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTg                                               >  2:6056093/1‑151 (MQ=255)
                                                                                                     tgGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTgg                                              >  2:8449040/1‑151 (MQ=255)
                                                                                                        ttcgACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGgcgc                                           >  1:3858109/4‑151 (MQ=255)
                                                                                                         tcgACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCt                                          <  1:8724600/149‑1 (MQ=255)
                                                                                                          cgACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTg                                         <  2:9449420/150‑1 (MQ=255)
                                                                                                             ccacggCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTgcgc                                      <  1:9001036/147‑1 (MQ=255)
                                                                                                                   cccGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGc                                <  2:10244162/151‑1 (MQ=255)
                                                                                                                        caCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGc                            >  1:5003038/2‑151 (MQ=255)
                                                                                                                         aCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCt                           <  2:2383346/151‑1 (MQ=255)
                                                                                                                               ctcgaccaGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTg                      <  2:10015903/144‑1 (MQ=255)
                                                                                                                                 cgaccaGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGcc                    <  1:9766770/146‑1 (MQ=255)
                                                                                                                                  gaccaGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCa                   >  1:5216408/5‑151 (MQ=255)
                                                                                                                                  gaccaGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCa                   <  2:4110810/147‑1 (MQ=255)
                                                                                                                                  gaccaGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCa                   <  1:132610/147‑1 (MQ=255)
                                                                                                                                     caGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATc                <  1:6636361/150‑1 (MQ=255)
                                                                                                                                     caGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATc                >  1:6736077/2‑151 (MQ=255)
                                                                                                                                     caGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATc                <  1:3500131/150‑1 (MQ=255)
                                                                                                                                       ggCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCg              <  2:5114641/151‑1 (MQ=255)
                                                                                                                                       ggCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCg              >  1:7153434/1‑151 (MQ=255)
                                                                                                                                        gCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCgg             <  1:2468200/151‑1 (MQ=255)
                                                                                                                                        gCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCgg             <  1:2325700/151‑1 (MQ=255)
                                                                                                                                         cACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGgc            <  2:1208751/151‑1 (MQ=255)
                                                                                                                                         cACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGgc            <  1:1561577/151‑1 (MQ=255)
                                                                                                                                         cACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGgc            >  1:4796016/1‑151 (MQ=255)
                                                                                                                                           cTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGgcgc          >  1:8104424/1‑151 (MQ=255)
                                                                                                                                           cTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGgcgc          <  1:5638334/151‑1 (MQ=255)
                                                                                                                                            tgtgCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGGCGcc         <  1:10940865/151‑1 (MQ=255)
                                                                                                                                            tgtgCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGGCGcc         <  2:4766507/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCTCGCTGCCAATCCGGCGCCACGCTTC  >  NC_008463/6069972‑6070269

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: ‑