Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 72,440 T→C K132E (AAG→GAG)  PA14_RS00315 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00846372,4400TC99.8% 1743.3 / NA 451K132E (AAG→GAG) PA14_RS00315hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (170/280);  minor base A (0/1);  total (170/281)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 454 total aligned reads displayed.

GGCAGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTTGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATCGCGCCAGAGCGGGC  >  NC_008463/72293‑72586
                                                                                                                                                   |                                                                                                                                                  
 gCAGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTcgcg                                                                                                                                                 <  1:5433167/150‑1 (MQ=255)
  cAGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCgg                                                                                                                                                <  2:1216202/150‑1 (MQ=255)
   aGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGc                                                                                                                                               <  1:7583142/150‑1 (MQ=255)
   aGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGc                                                                                                                                               <  1:1600230/150‑1 (MQ=255)
     cGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCcg                                                                                                                                             <  2:5290563/150‑1 (MQ=255)
        tCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCcgcgc                                                                                                                                          >  1:2884999/1‑150 (MQ=255)
           ggCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGcc                                                                                                                                       <  1:10454512/150‑1 (MQ=255)
           ggCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGcc                                                                                                                                       <  1:5754614/150‑1 (MQ=255)
             cgcAGCCGGCGAGGAGCCGCAGCGCGAGCGGGACGCGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCCTAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCaa                                                                                                                                     <  2:9487303/150‑1 (MQ=255)
               cAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAAtt                                                                                                                                   <  1:6713265/150‑1 (MQ=255)
               cAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAAtt                                                                                                                                   <  2:4330832/150‑1 (MQ=255)
                  ccGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTccc                                                                                                                                <  2:7544448/150‑1 (MQ=255)
                   cGGCGCGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGCGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCt                                                                                                                               <  2:7170174/150‑1 (MQ=255)
                   cGGCGAGGAGCAGGAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCTTATCGCCCTCGCCTTGCGCCGTTCAGTCCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCt                                                                                                                               <  1:10171544/150‑1 (MQ=255)
                    ggCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCtg                                                                                                                              >  2:6565391/1‑150 (MQ=255)
                    ggCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCtg                                                                                                                              >  1:10844068/1‑150 (MQ=255)
                    ggCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCAGGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCACTCGCGGCCGCGCGCCAATTCCCtg                                                                                                                              <  2:6798213/150‑1 (MQ=255)
                        aggagCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCtgttgt                                                                                                                          <  1:8722679/150‑1 (MQ=255)
                         ggagCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTc                                                                                                                         <  2:5682056/150‑1 (MQ=255)
                          gagCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTcc                                                                                                                        <  2:4770260/150‑1 (MQ=255)
                           agcagcagcGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCt                                                                                                                       <  1:6792235/150‑1 (MQ=255)
                              agcagcGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGt                                                                                                                    >  1:4368636/1‑150 (MQ=255)
                               gcagcGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTc                                                                                                                   >  1:9055285/1‑150 (MQ=255)
                                cagcGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCg                                                                                                                  <  1:958119/150‑1 (MQ=255)
                                 agcGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGt                                                                                                                 <  1:2374537/150‑1 (MQ=255)
                                    gcgAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTaga                                                                                                              <  2:4651397/150‑1 (MQ=255)
                                        gCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCtt                                                                                                          <  1:7512225/150‑1 (MQ=255)
                                         cGGGACGAGACGGGTGAGTGGACGCCTGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCAGTTGTCCTGGTCGTAGAGCGTc                                                                                                         <  2:8242431/150‑1 (MQ=255)
                                           ggACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTCGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCACTTCCCTGTTGTCCTGGTCGTAGAGCTTCGc                                                                                                       >  2:9549993/1‑150 (MQ=255)
                                            gACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGcc                                                                                                      >  1:1440196/1‑150 (MQ=255)
                                              cGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCa                                                                                                    >  1:8123439/1‑150 (MQ=255)
                                                agaCGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAgg                                                                                                  >  2:4536817/1‑150 (MQ=255)
                                                agaCGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTACAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAgg                                                                                                  <  1:10360491/150‑1 (MQ=255)
                                                 gaCGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAggg                                                                                                 >  1:6977579/1‑150 (MQ=255)
                                                  aCGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGa                                                                                                <  1:4965554/150‑1 (MQ=255)
                                                      gTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGa                                                                                            <  1:2162935/150‑1 (MQ=255)
                                                       tGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCGAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAAGTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGac                                                                                           <  1:6492093/150‑2 (MQ=255)
                                                        gAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAg                                                                                          <  1:3164698/150‑1 (MQ=255)
                                                         aGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCTCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGc                                                                                         <  1:205346/150‑1 (MQ=255)
                                                          gTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGcc                                                                                        >  2:1512346/1‑150 (MQ=255)
                                                              aCGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGAtg                                                                                    <  1:2754764/150‑1 (MQ=255)
                                                              aCGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGAtg                                                                                    <  1:2949461/150‑1 (MQ=255)
                                                                gCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGAtgtg                                                                                  >  2:10503943/1‑150 (MQ=255)
                                                                  aTGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGcc                                                                                >  1:4570938/1‑150 (MQ=255)
                                                                  aTGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCATGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGcc                                                                                >  2:396285/1‑150 (MQ=255)
                                                                  aTGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTAGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGAAGTGcc                                                                                <  2:1409094/150‑1 (MQ=255)
                                                                      ggcgggcTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCTCAGCGCCCCTCGCGGCCGCGCGCCTATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGc                                                                            <  1:8338878/150‑1 (MQ=255)
                                                                        cgggcTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCt                                                                          <  2:4961712/150‑1 (MQ=255)
                                                                         gggcTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTc                                                                         <  1:2577268/150‑1 (MQ=255)
                                                                           gcTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCa                                                                       <  2:10090583/150‑1 (MQ=255)
                                                                              ccGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGc                                                                    >  2:7341285/1‑150 (MQ=255)
                                                                              ccGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGc                                                                    <  2:6325915/150‑1 (MQ=255)
                                                                                ggAGGGCGTTCAGCTCCTATCATAGCGCCCTAGCCTTTCGCCGTTCATCCCTCCGCGGATCTCCCCACGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCgg                                                                  <  1:5174688/150‑1 (MQ=255)
                                                                                 gAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCGTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGt                                                                 >  1:10945432/1‑150 (MQ=255)
                                                                                  aGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTa                                                                <  2:1673567/150‑1 (MQ=255)
                                                                                   gggCGTTCAGCTCCGAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAg                                                               <  2:225167/150‑1 (MQ=255)
                                                                                    ggCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGGCGTAGAGCTTCGCCCCGGGACGGACGCCGATGTGCCCGGCCTCCAGGCGGTAgc                                                              >  2:5077690/1‑150 (MQ=255)
                                                                                       gTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCt                                                           <  1:9486529/150‑1 (MQ=255)
                                                                                        gtCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAAGTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTg                                                          <  2:24976/149‑1 (MQ=255)
                                                                                          cAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGc                                                        >  1:7312632/1‑150 (MQ=255)
                                                                                            gCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCGTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCg                                                      >  2:2838743/1‑150 (MQ=255)
                                                                                             cTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCgg                                                     <  1:478792/150‑1 (MQ=255)
                                                                                              tCCTAGCCTCGCGCCCTCGCCGTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGgc                                                    <  2:10260018/150‑1 (MQ=255)
                                                                                              tCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTAGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGgc                                                    >  2:9076177/1‑150 (MQ=255)
                                                                                                  aGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTg                                                >  2:7801937/1‑150 (MQ=255)
                                                                                                   gCATAGCGCCCTCGCCTTGCGCCGGTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTCGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGa                                               <  1:9214666/150‑1 (MQ=255)
                                                                                                    cATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGaa                                              <  1:9897195/150‑1 (MQ=255)
                                                                                                    cATACCGCCCGCGCCTTGCGCCGATCAGCCCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCACTCCCCGGCTGTCCTGGTCTTAGTGCTTCTCCCAGGGACGGACGCCGATGTGCCGGCCCTCCAGGCGGTAGCGCTGGCCGGCGCTGaa                                              <  1:4253442/150‑1 (MQ=255)
                                                                                                     aTAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAg                                             <  1:9676200/150‑1 (MQ=255)
                                                                                                      tAGCGCCCTCGCCTTGCGCCGGTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGGGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGt                                            <  2:8911648/150‑1 (MQ=255)
                                                                                                       aGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTc                                           <  1:7815380/150‑1 (MQ=255)
                                                                                                       aGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAAGTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTc                                           <  2:4095464/150‑1 (MQ=255)
                                                                                                         cgcCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCt                                         <  1:7969345/150‑1 (MQ=255)
                                                                                                            ccTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGa                                      >  2:10782627/1‑150 (MQ=255)
                                                                                                             cTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCGGGGGTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGaa                                     <  2:2644917/150‑1 (MQ=255)
                                                                                                               cGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACt                                   <  1:5285369/150‑1 (MQ=255)
                                                                                                                 ccTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTg                                 >  1:10092716/1‑150 (MQ=255)
                                                                                                                  cTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGa                                >  2:3027186/1‑150 (MQ=255)
                                                                                                                     gcgcCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGt                             <  1:1156631/150‑1 (MQ=255)
                                                                                                                     gcgcCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGt                             <  1:5611988/150‑1 (MQ=255)
                                                                                                                      cgcCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGtt                            <  1:414097/150‑1 (MQ=255)
                                                                                                                      cgcCGTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCGTGAACGTGAGGtt                            >  1:3470323/1‑150 (MQ=255)
                                                                                                                          gTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGc                        <  1:6216086/150‑1 (MQ=255)
                                                                                                                          gTTCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGc                        <  2:6389969/150‑1 (MQ=255)
                                                                                                                           ttCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCt                       <  2:1843372/150‑1 (MQ=255)
                                                                                                                            tCAGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCGTTGAACTTGAGGTTGAGCTg                      >  2:362790/1‑150 (MQ=255)
                                                                                                                              aGACCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGc                    >  1:372185/1‑150 (MQ=255)
                                                                                                                                aCCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCaa                  <  1:7200746/150‑1 (MQ=255)
                                                                                                                                aCCTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCaa                  >  1:8830866/1‑150 (MQ=255)
                                                                                                                                 ccTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAAt                 >  1:8667798/1‑150 (MQ=255)
                                                                                                                                 ccTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAAt                 <  2:5960092/150‑1 (MQ=255)
                                                                                                                                  cTTCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATc                <  1:10562192/150‑1 (MQ=255)
                                                                                                                                   ttCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATcg               >  1:5524777/1‑150 (MQ=255)
                                                                                                                                   ttCGCGCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGGAGAGCTGCGCCGCGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTCGCGCTGGCCGGCGCTGAAGTCCTTGAACGTGAGGGTGTGCGGGCAATcg               >  2:8464821/1‑150 (MQ=255)
                                                                                                                                     cgcgCAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATcgcg             <  1:6075038/150‑1 (MQ=255)
                                                                                                                                      gcgcAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATCgcgc            >  2:6172559/1‑150 (MQ=255)
                                                                                                                                        gcAGCGCCCCTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCATTGAACTTGAGGTTGAGCTGGCAATCGCGCCa          >  2:685813/1‑150 (MQ=255)
                                                                                                                                               cccTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGTAGTCCTTGAACGTGAGGTTGAGCTGGCAATCGCGCCAGAGCggg   >  1:6299329/1‑150 (MQ=255)
                                                                                                                                               cccTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATCGCGCCAGAGCggg   >  2:5484008/1‑150 (MQ=255)
                                                                                                                                                ccTCGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATCGCGCCAGAGCGGGc  <  2:9212436/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCAGCGGTCCGGCGCAGCCGGCGAGGAGCAGCAGCGCGAGCGGGACGAGACGGGTGAGTGGACGCATGGGGCGGGCTCCGGAGGGCGTTCAGCTCCTAGCATAGCGCCCTCGCCTTGCGCCGTTCAGACCTTCGCGCAGCGCCCCTTGCGGCCGCGCGCCAATTCCCTGTTGTCCTGGTCGTAGAGCTTCGCCCAGGGACGGAAGCCGATGTGCCCGGCCTCCAGGCGGTAGCGCTGGCCGGCGCTGAAGTCCTTGAACTTGAGGTTGAGCTGGCAATCGCGCCAGAGCGGGC  >  NC_008463/72293‑72586

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: