Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,150,987 Δ1 bp intergenic (‑192/+24) pvcA ← / ← PA14_RS14420 paerucumarin biosynthesis protein PvcA/asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,150,9850C.99.8% 2336.7 / ‑2.1 473intergenic (‑190/+26)pvcA/PA14_RS14420paerucumarin biosynthesis protein PvcA/asparaginase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (277/194);  total (279/194)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 506 total aligned reads displayed.

GGCAGCGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCATTCCACCA  >  NC_008463/3150835‑3151131
                                                                                                                                                      |                                                                                                                                                  
ggCAGCGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGc                                                                                                                                                    <  2:2329728/151‑1 (MQ=255)
ggCAGCGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGc                                                                                                                                                    >  1:11318073/1‑151 (MQ=255)
 gCAGCGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGATGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGcg                                                                                                                                                   <  2:2864021/151‑2 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGTTGGGCGCCGGGCACGGGCGGGCCGCCGGCGTGGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               >  1:2826420/1‑147 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGTCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               >  2:7426358/1‑147 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               <  2:9156859/151‑5 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               <  1:6913051/151‑5 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               <  1:4179482/151‑5 (MQ=255)
     cGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTAGCGCGGGGAGCTCGGAAGGAAGCGGccgccg                                                                                                                                               >  2:940405/1‑147 (MQ=255)
       gAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCcggc                                                                                                                                            >  1:10209593/1‑151 (MQ=255)
       gAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCcggc                                                                                                                                            >  2:5932022/1‑151 (MQ=255)
        aaGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCcggcg                                                                                                                                           <  1:10950080/151‑1 (MQ=255)
            cGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCggcgggcg                                                                                                                                       >  2:620061/1‑151 (MQ=255)
              aCATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGggcggc                                                                                                                                     >  1:5619620/1‑151 (MQ=255)
                ctCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGACACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTAGCGCGGGGACCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCg                                                                                                                                   <  2:9556022/150‑1 (MQ=255)
                aTCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCCCGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCg                                                                                                                                   <  1:9499363/151‑1 (MQ=255)
                   cTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAg                                                                                                                                >  2:10636597/1‑151 (MQ=255)
                    tgtgGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAgg                                                                                                                               >  2:1311005/1‑151 (MQ=255)
                    tgtgGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAgg                                                                                                                               >  1:6241686/1‑151 (MQ=255)
                     gtgGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAggg                                                                                                                              >  1:3509407/1‑151 (MQ=255)
                          cTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTc                                                                                                                         >  2:1993053/1‑151 (MQ=255)
                             gccgccGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGt                                                                                                                      >  1:8796000/1‑151 (MQ=255)
                              ccgccgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTc                                                                                                                     >  2:4294597/1‑151 (MQ=255)
                              ccgccgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTc                                                                                                                     >  1:6786121/1‑151 (MQ=255)
                               cgccgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCg                                                                                                                    <  1:3805257/151‑1 (MQ=255)
                                gccgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCgg                                                                                                                   >  2:5121188/1‑151 (MQ=255)
                                gccgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTAGCGTATAGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCgg                                                                                                                   >  2:7159447/1‑151 (MQ=255)
                                  cgGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCg                                                                                                                 >  1:8708645/1‑151 (MQ=255)
                                      ggATAAACGGGGGACGGCCGGAACCACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCCCGGGCGGGCCGCGGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTtct                                                                                                             <  1:320965/151‑1 (MQ=255)
                                       gATAACCCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTtctc                                                                                                            <  2:2477748/151‑1 (MQ=255)
                                        aTAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTcg                                                                                                           >  2:1871031/1‑151 (MQ=255)
                                        aTAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTcg                                                                                                           >  1:10333306/1‑151 (MQ=255)
                                           aaCCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTcgccg                                                                                                        >  2:1594995/1‑151 (MQ=255)
                                                gggAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTGGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAg                                                                                                   >  2:8415830/1‑151 (MQ=255)
                                                 ggAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAgg                                                                                                  <  2:884503/151‑1 (MQ=255)
                                                   aaGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTc                                                                                                >  2:2158457/1‑151 (MQ=255)
                                                    aGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCg                                                                                               >  1:5324323/1‑151 (MQ=255)
                                                     ggCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGGTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGa                                                                                              >  1:5814205/1‑151 (MQ=255)
                                                       ccGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAgg                                                                                            >  2:6470162/1‑151 (MQ=255)
                                                       ccGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAgg                                                                                            >  1:7838880/1‑151 (MQ=255)
                                                          gAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGGCGGCGTTCTCGCCGCACAGGTCGAGGGCg                                                                                         >  2:5345851/1‑151 (MQ=255)
                                                           aaaCACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGa                                                                                        >  2:4862185/1‑151 (MQ=255)
                                                            aaCACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGaa                                                                                       >  1:1590587/1‑151 (MQ=255)
                                                             acacTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGCGGCGAGGACGTTCTAGCGGAACGCGGGCAGGTCGGAAGGTAGAGGCCGCGGTAGGGCGGCAGCAGGGCGGGCAGTCGGCGTTCACGCCGCGCAGGTCGAGGGCGagc                                                                                      >  2:4566382/1‑149 (MQ=255)
                                                              cacTGTTCGCGGAAGCTTGACGATGGGCGCCGGGCACGGGGGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAAcc                                                                                     <  1:5891428/151‑1 (MQ=255)
                                                                cTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAg                                                                                   >  2:1397553/1‑151 (MQ=255)
                                                                    tCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGtt                                                                               >  2:2996738/1‑151 (MQ=255)
                                                                      gCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTgc                                                                             <  2:10375917/151‑1 (MQ=255)
                                                                        tGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTgcgc                                                                           >  1:168217/1‑151 (MQ=255)
                                                                           aGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAGGGAAGCGGCCGACGGCGGGCGGCCGCAGGTCGGTCAGTCGGCGTTCTCGACGCACAGGTCGAGGGCGAACCAGCGTTGCGCTtc                                                                        >  1:8119169/1‑151 (MQ=255)
                                                                             cTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCTCAGGTCGAGGGCGACCCAGCGTTGCGcgtcgt                                                                      >  1:10083365/1‑151 (MQ=255)
                                                                             cTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTtcgt                                                                      <  1:2728749/151‑1 (MQ=255)
                                                                               tGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTtcgtcg                                                                    <  2:827490/151‑1 (MQ=255)
                                                                               tGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGAATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTtcgtcg                                                                    <  1:6553781/151‑1 (MQ=255)
                                                                                 aCGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGGATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGACCCAGCGTGGCGCTTCGTcgcg                                                                  >  2:3541091/1‑151 (MQ=255)
                                                                                 aCGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTAGTCGCGTATCGCGGGGAGCTCGGCAGGAAGCGGCACCCCGCGGGAGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGGACAGGTCGAGGGCGATCCTGCGTTGGGCTTCGTcgcg                                                                  >  1:4520194/1‑151 (MQ=255)
                                                                                  cGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGt                                                                 <  2:6860690/151‑1 (MQ=255)
                                                                                      gggCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGa                                                             >  1:11108980/1‑151 (MQ=255)
                                                                                      gggCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGa                                                             <  1:8856513/151‑1 (MQ=255)
                                                                                        gcgcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAgg                                                           <  1:114981/151‑1 (MQ=255)
                                                                                        gcgcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAgg                                                           >  1:4741822/1‑151 (MQ=255)
                                                                                         cgcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGc                                                          >  2:10787302/1‑151 (MQ=255)
                                                                                         cgcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGc                                                          >  2:4976462/1‑151 (MQ=255)
                                                                                          gcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTGGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGGTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGcc                                                         >  1:10622144/1‑151 (MQ=255)
                                                                                          gcCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGcc                                                         <  2:8609009/151‑1 (MQ=255)
                                                                                            cGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGGGGGGAGCTCGGAAGGAAGCGGCCGCCGGCCGGAGGACGCAGGGCGGTCAGTCGGAGTTCTCGCCGCACAGGTCGAGGGCGTACCAGCGTTGCGCTGGGTCGCGTCCGAGGCCgg                                                       >  1:4148928/1‑151 (MQ=255)
                                                                                            cGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCgg                                                       >  1:7513302/1‑151 (MQ=255)
                                                                                             gggCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGgc                                                      <  1:1358934/151‑1 (MQ=255)
                                                                                             gggCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGgc                                                      >  2:421196/1‑151 (MQ=255)
                                                                                                   gggcgggcCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAg                                                >  2:10548605/1‑151 (MQ=255)
                                                                                                   gggcgggcCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAg                                                >  2:3903701/1‑151 (MQ=255)
                                                                                                    ggcgggcCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGc                                               <  1:8980590/151‑1 (MQ=255)
                                                                                                       gggcCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCTGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAgg                                            <  1:916653/151‑1 (MQ=255)
                                                                                                         gccgccGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGcc                                          <  2:8489192/151‑1 (MQ=255)
                                                                                                           cgccgGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGGTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGa                                        >  2:7733597/1‑151 (MQ=255)
                                                                                                            gccgGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGaa                                       >  1:8224251/1‑151 (MQ=255)
                                                                                                             ccgGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAag                                      >  2:11266058/1‑151 (MQ=255)
                                                                                                                gCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGc                                   >  1:7691937/1‑151 (MQ=255)
                                                                                                                 cGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCt                                  <  2:3189118/151‑1 (MQ=255)
                                                                                                                  gTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCtt                                 >  1:4787735/1‑151 (MQ=255)
                                                                                                                  gTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCtt                                 <  2:1194967/151‑1 (MQ=255)
                                                                                                                   ttGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTg                                <  2:4134929/151‑1 (MQ=255)
                                                                                                                   ttGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTg                                <  1:1938566/151‑1 (MQ=255)
                                                                                                                     gTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGcc                              >  2:6805396/1‑151 (MQ=255)
                                                                                                                        gcgTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAg                           >  1:5465868/1‑151 (MQ=255)
                                                                                                                          gTATCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAggg                         <  2:10149433/151‑1 (MQ=255)
                                                                                                                            ctCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGgcc                       <  2:11395465/150‑1 (MQ=255)
                                                                                                                            aTCGCGGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGgcc                       >  2:6098546/1‑151 (MQ=255)
                                                                                                                              cgcgGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGgccgc                     >  1:7244860/1‑151 (MQ=255)
                                                                                                                                cgGGGAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCt                   <  1:3008856/151‑1 (MQ=255)
                                                                                                                                 ggggAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGACCCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTc                  >  1:3621972/1‑151 (MQ=255)
                                                                                                                                 ggggAGCTCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTc                  <  1:1267168/151‑1 (MQ=255)
                                                                                                                                        tCGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTc           >  2:11086087/1‑151 (MQ=255)
                                                                                                                                        gcGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTc           <  2:921803/150‑1 (MQ=255)
                                                                                                                                         cGGAAGGAAGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCa          >  2:2717802/1‑151 (MQ=255)
                                                                                                                                          ggaaggaaGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCAt         <  1:9401775/151‑1 (MQ=255)
                                                                                                                                          ggaaggaaGCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCAt         >  2:9783185/1‑151 (MQ=255)
                                                                                                                                             aggaagCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCATTcc      >  1:2092128/1‑151 (MQ=255)
                                                                                                                                                aagCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCATTccacc   >  1:2575251/1‑151 (MQ=255)
                                                                                                                                                aagCGGCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCATTccacc   <  2:7510903/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GGCAGCGGAAGGCGACATCCTGTGGACTCGCCGCCGGGGGATAAACCGGGGAAGGCCGGAAACACTGTTCGCTGAAGCTTGACGATGGGCGCCGGGCACGGGCGGGCCGCCGGCGTTGTCGCGTATCGCGGGGAGCTCGGAAGGAAGCGGCCCGCCGGCGGGCGGCCGCAGGGCGGTCAGTCGGCGTTCTCGCCGCACAGGTCGAGGGCGAACCAGCGTTGCGCTTCGTCGCGTCCGAGGCCGGCGCCGAGCAGGCCGAAGAGCTTGCCCAGGGCCGCCTCGCGGGTCATTCCACCA  >  NC_008463/3150835‑3151131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: