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

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

Only 100 of 614 total aligned reads displayed.

CGACAGCGACTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGC  >  NC_008463/6069930‑6070226
                                                                                                                                                    |                                                                                                                                                    
  aCAGCGACTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGAcc                                                                                                                                                  <  1:5615195/150‑1 (MQ=255)
  aCAGCGACTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGAcc                                                                                                                                                  <  2:4532640/150‑1 (MQ=255)
  aCAGCGACTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGAcc                                                                                                                                                  <  1:3840509/150‑1 (MQ=255)
        aCTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGcc                                                                                                                                            <  2:7022011/150‑1 (MQ=255)
        aCTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGcc                                                                                                                                            <  1:8663117/150‑1 (MQ=255)
           tCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGc                                                                                                                                         <  1:6894091/150‑1 (MQ=255)
             cTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCac                                                                                                                                       <  2:8499248/150‑2 (MQ=255)
                 ccGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCacctgc                                                                                                                                   <  1:5689488/150‑6 (MQ=255)
                   gTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCacctgctc                                                                                                                                 >  1:7532956/1‑143 (MQ=255)
                       aGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGAcc                                                                                                                            >  2:126502/1‑150 (MQ=255)
                        ggCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCa                                                                                                                           <  2:1446015/150‑1 (MQ=255)
                         gccgccCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCag                                                                                                                          <  1:8128316/150‑2 (MQ=255)
                           cgccCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCaggc                                                                                                                        <  1:2619979/150‑4 (MQ=255)
                              ccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACt                                                                                                                    <  2:4359864/150‑1 (MQ=255)
                              ccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACt                                                                                                                    >  2:6732532/1‑150 (MQ=255)
                                  cAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGc                                                                                                                >  1:9090922/1‑150 (MQ=255)
                                  cAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGc                                                                                                                <  2:214713/150‑1 (MQ=255)
                                   aGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGcc                                                                                                               <  2:10443099/150‑1 (MQ=255)
                                    ggCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCt                                                                                                              <  2:10073526/150‑1 (MQ=255)
                                       cGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGtt                                                                                                           <  1:7395261/150‑1 (MQ=255)
                                        ggCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGttc                                                                                                          >  1:2966085/1‑149 (MQ=255)
                                         gcagcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGttcg                                                                                                         <  1:9612858/150‑3 (MQ=255)
                                         gcagcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGttcg                                                                                                         <  2:9726201/150‑3 (MQ=255)
                                          cagcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCgg                                                                                                        >  2:680570/1‑150 (MQ=255)
                                           agcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGc                                                                                                       >  2:396235/1‑150 (MQ=255)
                                            gcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCa                                                                                                      <  2:5566155/150‑1 (MQ=255)
                                             caTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCAt                                                                                                     <  2:1871064/150‑1 (MQ=255)
                                                  gCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGc                                                                                                >  1:100338/1‑150 (MQ=255)
                                                   cccGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCa                                                                                               <  2:3335821/150‑1 (MQ=255)
                                                    ccGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGTCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAg                                                                                              >  1:840396/1‑150 (MQ=255)
                                                     cGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGc                                                                                             <  2:9624744/150‑1 (MQ=255)
                                                     cGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGc                                                                                             <  2:8894325/150‑1 (MQ=255)
                                                      ggTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCg                                                                                            >  2:2691457/1‑150 (MQ=255)
                                                       gTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGt                                                                                           >  2:10548589/1‑150 (MQ=255)
                                                         gcgcTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTAt                                                                                         <  1:658496/150‑1 (MQ=255)
                                                           gcTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGc                                                                                       <  1:216436/150‑1 (MQ=255)
                                                            cTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCa                                                                                      <  2:4873212/150‑1 (MQ=255)
                                                            cTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCa                                                                                      <  1:5217836/150‑1 (MQ=255)
                                                             tACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCAc                                                                                     <  1:5105316/150‑1 (MQ=255)
                                                              aCTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACg                                                                                    >  2:61017/1‑150 (MQ=255)
                                                                ttCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGa                                                                                  <  1:10267544/150‑1 (MQ=255)
                                                                ttCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGa                                                                                  >  2:3012389/1‑150 (MQ=255)
                                                                   gAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGacac                                                                               >  1:9807702/1‑150 (MQ=255)
                                                                    aGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACAcc                                                                              <  1:6803241/150‑1 (MQ=255)
                                                                     gTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCt                                                                             <  1:10547496/150‑1 (MQ=255)
                                                                     gTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCt                                                                             >  1:4283807/1‑150 (MQ=255)
                                                                       cGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                           >  1:1937483/1‑150 (MQ=255)
                                                                       cGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                           <  2:9271228/150‑1 (MQ=255)
                                                                         gCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgcgc                                                                         <  1:4879190/150‑1 (MQ=255)
                                                                          ccATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGcc                                                                        >  2:8109943/1‑150 (MQ=255)
                                                                          ccATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGcc                                                                        <  2:3606781/150‑1 (MQ=255)
                                                                          ccATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGcc                                                                        <  1:4515573/150‑1 (MQ=255)
                                                                             tCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGc                                                                     <  1:3494382/150‑1 (MQ=255)
                                                                               gAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTg                                                                   >  1:6172346/1‑150 (MQ=255)
                                                                                 gcgcTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgcg                                                                 <  1:3285371/150‑1 (MQ=255)
                                                                                 gcgcTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgcg                                                                 <  1:7098542/150‑1 (MQ=255)
                                                                                  cgcTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGACCACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgcgc                                                                <  2:7916664/150‑1 (MQ=255)
                                                                                      tccgtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAgg                                                            >  2:8800124/1‑150 (MQ=255)
                                                                                       ccgtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAggg                                                           >  1:7994084/1‑150 (MQ=255)
                                                                                         gtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcg                                                         <  1:4963843/150‑1 (MQ=255)
                                                                                          tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                        >  1:4607363/1‑150 (MQ=255)
                                                                                          tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                        >  1:10673627/1‑150 (MQ=255)
                                                                                          tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                        <  2:10908778/150‑1 (MQ=255)
                                                                                          tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGGTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                        >  2:4718874/1‑150 (MQ=255)
                                                                                           ccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcg                                                       <  1:7653216/150‑1 (MQ=255)
                                                                                            cgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcgc                                                      >  1:2322351/1‑150 (MQ=255)
                                                                                             gCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCCCTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCt                                                     >  2:2015799/1‑150 (MQ=255)
                                                                                             gCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCt                                                     <  1:11056252/150‑1 (MQ=255)
                                                                                             gCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCt                                                     >  2:5046156/1‑150 (MQ=255)
                                                                                                 gTCCGCGACCGCTGGCGCGATCTCGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGAcc                                                 <  1:5512268/150‑1 (MQ=255)
                                                                                                  tCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcg                                                <  1:10771529/150‑1 (MQ=255)
                                                                                                  tCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcg                                                >  2:286006/1‑150 (MQ=255)
                                                                                                    cgcgACCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcgac                                              >  1:9383203/1‑150 (MQ=255)
                                                                                                       gACCGCTGGCGCGACCACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCCGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCgacgac                                           >  2:2531203/1‑150 (MQ=255)
                                                                                                        aCCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACt                                          >  2:7262076/1‑150 (MQ=255)
                                                                                                        aCCGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACt                                          >  1:9001678/1‑150 (MQ=255)
                                                                                                         ccGCTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGCGGGCGCGCTGACCGACGACTg                                         >  2:6295696/1‑150 (MQ=255)
                                                                                                            cTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGtt                                      <  1:9251130/150‑1 (MQ=255)
                                                                                                            cTGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGtt                                      >  2:143001/1‑150 (MQ=255)
                                                                                                             tGGCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTc                                     >  2:8283634/1‑150 (MQ=255)
                                                                                                              ggCGCGAACACGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCg                                    >  2:5881383/1‑150 (MQ=255)
                                                                                                                    cacacGTCCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCCACTGCTCCTCCAGGCACTGTGCCTGTTCGCCATAGCGAAGCGTATGCACGGACACCTGCTCCGTCTGCGCGTGGGCGCGCTGACCGACGACTGGTTCGAGGTcc                              <  2:2276259/149‑1 (MQ=255)
                                                                                                                      cacGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTg                            <  1:10371146/150‑1 (MQ=255)
                                                                                                                        cGACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCt                          <  1:5828829/150‑1 (MQ=255)
                                                                                                                         gACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTc                         <  1:1286145/150‑1 (MQ=255)
                                                                                                                         gACCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTc                         <  1:6476565/150‑1 (MQ=255)
                                                                                                                          aCCTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCg                        >  2:5375105/1‑150 (MQ=255)
                                                                                                                           ccTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCTTACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACGGGTTCGAGGTCCTGCTCGa                       >  1:1509838/1‑150 (MQ=255)
                                                                                                                           ccTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGa                       <  2:9404485/150‑1 (MQ=255)
                                                                                                                            cTGCCCGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGAc                      <  1:7254321/150‑1 (MQ=255)
                                                                                                                               cccGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTAt                   <  2:7570213/150‑1 (MQ=255)
                                                                                                                                ccGGCTCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATc                  <  2:767839/150‑1 (MQ=255)
                                                                                                                                  cggctCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATcc                 >  1:4139475/5‑150 (MQ=255)
                                                                                                                                   ggctCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCg                <  1:5963281/147‑1 (MQ=255)
                                                                                                                                   ggctCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGGGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCg                <  2:2410860/147‑1 (MQ=255)
                                                                                                                                     ctCCGGGCTGTGGTTCGACCTCGGCCCGCCCCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCa              <  1:2878600/149‑1 (MQ=255)
                                                                                                                                      tCCGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCAc             >  1:3675148/1‑150 (MQ=255)
                                                                                                                                        cGGGCTGTGGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTg           >  1:419899/1‑150 (MQ=255)
                                                                                                                                              gtgGTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCCGCGTATGCACGGACACCTGCGCCGGCTTCGCGAGGGCGCGCTGACCGACGACAGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGt     >  1:8929842/1‑150 (MQ=255)
                                                                                                                                                ggTTCGACCTCGGCCCGCACCTGCTCGACCAGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTgg   <  1:7862710/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACAGCGACTTCCTCGCCGTGGAGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGC  >  NC_008463/6069930‑6070226

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: