Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,390,498 A→C pseudogene (464/1997 nt) PA14_RS15455 → alkaline phosphatase D family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,390,4980AC99.8% 2018.6 / NA 524pseudogene (464/1997 nt)PA14_RS15455alkaline phosphatase D family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (240/283);  minor base G (1/0);  total (241/283)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 524 total aligned reads displayed.

GCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCTGTACGG  >  NC_008463/3390349‑3390644
                                                                                                                                                     |                                                                                                                                                  
gCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGa                                                                                                                                                  <  1:781060/150‑1 (MQ=255)
 cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAc                                                                                                                                                 <  1:4749301/150‑1 (MQ=255)
 cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAc                                                                                                                                                 <  1:4910159/150‑1 (MQ=255)
  ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                                <  1:3810648/150‑1 (MQ=255)
  ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                                >  1:8529814/1‑150 (MQ=255)
  ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                                <  1:1289984/150‑1 (MQ=255)
    gcgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCgg                                                                                                                                              <  2:9048256/150‑1 (MQ=255)
     cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgc                                                                                                                                             <  1:5020364/150‑1 (MQ=255)
     cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgc                                                                                                                                             <  2:1297121/150‑1 (MQ=255)
     cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgc                                                                                                                                             <  2:2482732/150‑1 (MQ=255)
         aaccGCCGAATGGGCGCCGATTCTGCTCTACCCCGCTCGCCAGCGCGCCAGCCTGGTGTTCAAGCCACTTCTGGCCCAGTGTCTCCACGGCTCATGCCGCAACCCGAATCACCCGGCGGCGGATGGCCTGCTGCAGGTCTACCGgctggt                                                                                                                                         <  2:7089467/147‑3 (MQ=255)
         aaTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                                         >  1:3993430/1‑150 (MQ=255)
         aaTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                                         <  1:757981/150‑1 (MQ=255)
          aTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGATGCGCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCAGCTGCAGGTCGACCCgctgcgg                                                                                                                                        >  1:9997486/1‑148 (MQ=255)
             gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggag                                                                                                                                     <  1:5565269/150‑1 (MQ=255)
             gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggag                                                                                                                                     <  1:8297991/150‑1 (MQ=255)
             gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggag                                                                                                                                     >  1:10646919/1‑150 (MQ=255)
              ccGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagg                                                                                                                                    <  2:10176615/150‑1 (MQ=255)
               cgcaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagga                                                                                                                                   <  2:5324205/147‑1 (MQ=255)
                  aTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaaa                                                                                                                                >  1:10154004/1‑150 (MQ=255)
                   tGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAc                                                                                                                               >  1:609142/1‑150 (MQ=255)
                     ggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACct                                                                                                                             >  2:10537297/1‑150 (MQ=255)
                      gcgcCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctg                                                                                                                            <  1:533521/150‑1 (MQ=255)
                       cgccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgc                                                                                                                           <  1:2832786/150‑1 (MQ=255)
                       cgccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCAGCCGCAATCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACcagc                                                                                                                           >  1:1144793/1‑150 (MQ=255)
                           gcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCa                                                                                                                       <  1:1477181/150‑1 (MQ=255)
                             aTCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCCGGTCGACCGGCTGCTGGAGGAAATCCTGCTGCAgg                                                                                                                     >  1:432229/1‑150 (MQ=255)
                             aTCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAgg                                                                                                                     <  2:8946212/150‑1 (MQ=255)
                             aTCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAgg                                                                                                                     <  2:7755316/150‑1 (MQ=255)
                              tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                                    <  2:926345/150‑1 (MQ=255)
                               cTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCTCCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGcc                                                                                                                   >  1:1656070/1‑150 (MQ=255)
                               cTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGcc                                                                                                                   <  1:9662525/150‑1 (MQ=255)
                                tGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTACTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCg                                                                                                                  >  2:1654738/1‑150 (MQ=255)
                                  ctctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGaa                                                                                                                <  2:2362159/150‑1 (MQ=255)
                                  ctctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGaa                                                                                                                <  2:8492893/150‑1 (MQ=255)
                                  ctctACCCCGGGCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGaa                                                                                                                <  1:7385310/150‑1 (MQ=255)
                                   tctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGaaa                                                                                                               <  2:6931319/150‑1 (MQ=255)
                                   tctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGaaa                                                                                                               >  1:4138308/1‑150 (MQ=255)
                                    ctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAg                                                                                                              >  2:3720556/1‑150 (MQ=255)
                                     tACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgc                                                                                                             <  1:6699822/150‑1 (MQ=255)
                                         ccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgcc                                                                                                         <  2:1147019/150‑1 (MQ=255)
                                         ccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgcc                                                                                                         <  1:8635161/150‑1 (MQ=255)
                                          cGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGccc                                                                                                        >  2:1098597/1‑150 (MQ=255)
                                            gTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGCCCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGACAACCTGCTGCAGGCCGAAAGCCGCCCgg                                                                                                      <  2:4342373/150‑1 (MQ=255)
                                              cGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCg                                                                                                    >  2:65950/1‑150 (MQ=255)
                                               gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGt                                                                                                   <  2:3404742/150‑1 (MQ=255)
                                               gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGACCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGt                                                                                                   >  2:3121064/1‑150 (MQ=255)
                                                ccAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGtt                                                                                                  >  1:10267155/1‑150 (MQ=255)
                                                  aGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGc                                                                                                >  2:3957058/1‑150 (MQ=255)
                                                  aGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGc                                                                                                >  2:6372972/1‑150 (MQ=255)
                                                     gcgcCAGCCTGGTGCTCACGCCACGTCTCGACCAGTTGCTCCCCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtga                                                                                             <  2:3035688/150‑1 (MQ=255)
                                                        ccAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                                          >  1:1728865/1‑150 (MQ=255)
                                                        ccAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                                          <  1:3196252/150‑1 (MQ=255)
                                                         cAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgat                                                                                         <  2:4931624/150‑1 (MQ=255)
                                                         cAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgat                                                                                         <  1:8065738/150‑1 (MQ=255)
                                                             cTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgc                                                                                     <  1:2212226/150‑1 (MQ=255)
                                                                  gCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGa                                                                                >  2:8402686/1‑150 (MQ=255)
                                                                   cTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGAc                                                                               <  1:2664859/150‑1 (MQ=255)
                                                                       aGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAgg                                                                           <  1:10895594/150‑1 (MQ=255)
                                                                         ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTc                                                                         <  2:2903472/150‑1 (MQ=255)
                                                                          cACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCt                                                                        <  2:5734762/150‑1 (MQ=255)
                                                                               ctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGcc                                                                   >  1:9529349/1‑150 (MQ=255)
                                                                               ctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGcc                                                                   <  2:10723413/150‑1 (MQ=255)
                                                                               ctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGcc                                                                   >  1:9012964/1‑150 (MQ=255)
                                                                                  gACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgac                                                                <  2:1862423/150‑1 (MQ=255)
                                                                                   aCCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgacg                                                               <  1:6444411/150‑1 (MQ=255)
                                                                                   aCCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgacg                                                               >  2:9535225/1‑150 (MQ=255)
                                                                                        ttGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTc                                                          >  1:4541830/1‑150 (MQ=255)
                                                                                          gCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGTAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGc                                                        >  1:2028154/1‑150 (MQ=255)
                                                                                            tCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCg                                                      <  2:8249209/150‑1 (MQ=255)
                                                                                             ccACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                                     >  1:3486318/1‑150 (MQ=255)
                                                                                             ccACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCGGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                                     >  1:3267414/1‑150 (MQ=255)
                                                                                                cGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGTGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGcc                                                  >  2:4843513/1‑150 (MQ=255)
                                                                                                  gCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGa                                                <  1:5523074/150‑1 (MQ=255)
                                                                                                    tCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATg                                              <  1:9234317/150‑1 (MQ=255)
                                                                                                     ccTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGc                                             <  1:6816963/150‑1 (MQ=255)
                                                                                                        gCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtgg                                          >  1:10798943/1‑150 (MQ=255)
                                                                                                            cAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGc                                      >  1:7160821/1‑150 (MQ=255)
                                                                                                            cAAGCCGCATCACCCGGCGGCCGATGGCCTGCAGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGc                                      >  1:8797935/1‑150 (MQ=255)
                                                                                                             aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcc                                     <  1:8097627/150‑1 (MQ=255)
                                                                                                             aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcc                                     >  2:6503927/1‑150 (MQ=255)
                                                                                                               gCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccat                                   <  2:8747143/150‑1 (MQ=255)
                                                                                                                 cGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatcc                                 >  1:5761878/1‑150 (MQ=255)
                                                                                                                          ccggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTg                         >  2:10349370/3‑150 (MQ=255)
                                                                                                                          ccggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTCCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTg                         >  2:5181197/3‑150 (MQ=255)
                                                                                                                             gcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATc                      >  2:6173475/1‑150 (MQ=255)
                                                                                                                                gcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGCGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGAGccgc                   >  2:2645866/1‑150 (MQ=255)
                                                                                                                                 ccGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgcc                  >  2:4624546/1‑150 (MQ=255)
                                                                                                                                   cgatggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgccg                 >  2:2129247/5‑150 (MQ=255)
                                                                                                                                    gatggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGTTGCTGGTGGCCGTCCATGGGCAGATCcgccgc                >  2:4546196/4‑150 (MQ=255)
                                                                                                                                      tggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCt              >  2:8052258/2‑150 (MQ=255)
                                                                                                                                      gggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCt              <  1:2992209/150‑1 (MQ=255)
                                                                                                                                       ggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTc             <  2:9712111/150‑1 (MQ=255)
                                                                                                                                       ggCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTc             <  2:9182753/150‑1 (MQ=255)
                                                                                                                                         ccTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCgg           <  2:10643196/150‑1 (MQ=255)
                                                                                                                                           tgctgcAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATCGGCTGATCCGCCGCCTCGGcc         <  1:5969544/150‑1 (MQ=255)
                                                                                                                                            gctgcAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCt        <  1:5229924/150‑1 (MQ=255)
                                                                                                                                            gctgcAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCt        <  2:5313057/150‑1 (MQ=255)
                                                                                                                                             ctgcAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCTg       <  2:7304768/150‑1 (MQ=255)
                                                                                                                                                  ggTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCTGTACgg  >  2:6020782/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCGCCTCGGCCTGTACGG  >  NC_008463/3390349‑3390644

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: