Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,644,364 Δ1 bp intergenic (‑115/‑120) sohB ← / → PA14_RS16560 protease SohB/histidine phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,644,3610G.99.8% 2878.6 / ‑3.2 583intergenic (‑112/‑123)sohB/PA14_RS16560protease SohB/histidine phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (249/332);  minor base C (0/1);  total (249/334)
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 612 total aligned reads displayed.

CTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTC  >  NC_008463/3644212‑3644505
                                                                                                                                                     |                                                                                                                                                
 tGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgg                                                                                                                                                 >  2:7786982/1‑150 (MQ=255)
  ggCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAggg                                                                                                                                                <  2:3370163/150‑1 (MQ=255)
  ggCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAggg                                                                                                                                                >  2:4967565/1‑150 (MQ=255)
   gCCCGAAAGCCCGCATCTTCGGCAATGACCTCCACACAACCCCCACTCTCCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggc                                                                                                                                               <  1:5535308/150‑2 (MQ=255)
   gCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggc                                                                                                                                               <  1:90737/150‑2 (MQ=255)
   gCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggc                                                                                                                                               <  1:2033719/150‑2 (MQ=255)
    ccAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                              >  1:3308580/1‑148 (MQ=255)
     cAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggcac                                                                                                                                             >  1:3715882/1‑147 (MQ=255)
        aaaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGcc                                                                                                                                         <  1:2745381/150‑1 (MQ=255)
                  tatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCAGGAGCAGcc                                                                                                                               >  2:7062503/1‑150 (MQ=255)
                   atTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCg                                                                                                                              >  1:10064105/1‑150 (MQ=255)
                   atTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCg                                                                                                                              >  2:10332454/1‑150 (MQ=255)
                    ttCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGa                                                                                                                             <  1:10216959/150‑1 (MQ=255)
                     tCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAt                                                                                                                            >  2:9498004/1‑150 (MQ=255)
                     tCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAt                                                                                                                            <  2:3011457/150‑1 (MQ=255)
                       ggCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                          >  1:4229227/1‑150 (MQ=255)
                          aaTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                       >  2:10833178/1‑150 (MQ=255)
                          aaTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                       >  1:3989734/1‑150 (MQ=255)
                           aTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGc                                                                                                                      <  1:7030194/150‑1 (MQ=255)
                            tGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCa                                                                                                                     >  2:7974794/1‑150 (MQ=255)
                              aaCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGc                                                                                                                   <  2:6055816/150‑1 (MQ=255)
                               aCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                  <  2:9637207/150‑1 (MQ=255)
                                 tCCACACACCCCCCTCGCTACTGTGGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                <  1:159825/150‑1 (MQ=255)
                                 tCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                >  1:8031477/1‑150 (MQ=255)
                                 tCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                <  1:6538778/150‑1 (MQ=255)
                                  ccACCCAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCAGACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTGTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                               <  2:6203271/150‑1 (MQ=255)
                                  ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                               <  2:9021376/150‑1 (MQ=255)
                                  ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                               >  1:890556/1‑150 (MQ=255)
                                  ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                               <  1:10676559/150‑1 (MQ=255)
                                  ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                               >  1:3256730/1‑150 (MQ=255)
                                    acacaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGa                                                                                                             <  1:216671/150‑1 (MQ=255)
                                    acacaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGa                                                                                                             <  1:5784816/150‑1 (MQ=255)
                                     cacaACCCCCTCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGAc                                                                                                            <  2:4175962/150‑1 (MQ=255)
                                       caACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAt                                                                                                          <  1:1179482/150‑1 (MQ=255)
                                          cccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCa                                                                                                       <  2:2832303/150‑1 (MQ=255)
                                            cccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaaa                                                                                                     <  2:384108/150‑1 (MQ=255)
                                             ccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAAc                                                                                                    <  1:734569/150‑1 (MQ=255)
                                              ctctctACTGTTGACGCCCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                   <  2:8628708/150‑1 (MQ=255)
                                              ctctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                   >  1:5079690/1‑150 (MQ=255)
                                              ctatctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                   >  2:6853857/4‑150 (MQ=255)
                                                   tACTGTTGACGCACCCGTCTACGGGCGACACGAATCCGCAATCATACCCAGGCCACCGACGGGCGGCCCGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAATCAAGCGt                                                                                              >  2:2258827/1‑150 (MQ=255)
                                                    aCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTa                                                                                             <  1:9821725/150‑1 (MQ=255)
                                                     cTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAt                                                                                            <  1:9278192/150‑1 (MQ=255)
                                                     cTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAt                                                                                            >  1:10784182/1‑150 (MQ=255)
                                                      tGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATg                                                                                           >  2:745336/1‑150 (MQ=255)
                                                      tGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATg                                                                                           <  1:1860888/150‑1 (MQ=255)
                                                       gTTGACGCACCCGTCTACGGGCGACCCGATTGCGCAAGCATACCCCTGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTTCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGt                                                                                          <  1:1256826/150‑1 (MQ=255)
                                                       gTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGt                                                                                          >  1:8760678/1‑150 (MQ=255)
                                                        ttGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtt                                                                                         >  1:8368823/1‑150 (MQ=255)
                                                          gACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttt                                                                                       <  2:8586786/150‑1 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGTCCACCGAAGGGCGGCCAGCGATGCCACTGCCTGGAGATCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                     <  2:7931532/150‑1 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                     <  2:9333125/150‑1 (MQ=255)
                                                             gCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTc                                                                                    <  2:1025656/150‑1 (MQ=255)
                                                               aCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGt                                                                                  >  1:10943096/1‑150 (MQ=255)
                                                                 ccGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTg                                                                                <  2:5504318/150‑1 (MQ=255)
                                                                  cGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGa                                                                               <  2:9850252/150‑1 (MQ=255)
                                                                   gTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGGATGACATTCAAACAAGCGTATGTTTTTTCGTTGAc                                                                              >  2:239466/1‑150 (MQ=255)
                                                                    tCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                             >  2:485241/1‑150 (MQ=255)
                                                                    tCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                             >  1:7098851/1‑150 (MQ=255)
                                                                    tCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                             >  1:5237643/1‑150 (MQ=255)
                                                                     cTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAg                                                                            >  1:854741/1‑150 (MQ=255)
                                                                      tACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGc                                                                           <  1:10430467/150‑1 (MQ=255)
                                                                      tACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGc                                                                           >  2:1754704/1‑150 (MQ=255)
                                                                        cGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGccc                                                                         <  2:671832/150‑1 (MQ=255)
                                                                         gggCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCt                                                                        <  2:3964592/150‑1 (MQ=255)
                                                                           gCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGTAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCttg                                                                      <  2:8265352/150‑3 (MQ=255)
                                                                             gACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGt                                                                    >  2:145525/1‑150 (MQ=255)
                                                                               cacGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGa                                                                  <  2:10001309/150‑1 (MQ=255)
                                                                                  gATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATcc                                                               >  1:8448572/1‑150 (MQ=255)
                                                                                   aTTGCGCAAGCATACCCCGGCCACCGACGGGCGGCCAGCGCTGCCACTGCCTGGAGAGCCTTTGCCGGGCACGCCCGGAGCAGCCGATACAGGCAGCGAATGACATTCAAACAAGCGTATGATTTTTCGTTGACAGCCCTGGGTGATCCt                                                              <  2:4574076/150‑1 (MQ=255)
                                                                                    ttGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTc                                                             <  1:7823764/150‑1 (MQ=255)
                                                                                     tGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCa                                                            <  1:4298486/150‑1 (MQ=255)
                                                                                      gcgcAAGCATACCCCGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAt                                                           <  1:8121922/150‑1 (MQ=255)
                                                                                      gcgcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAt                                                           <  2:4692057/150‑1 (MQ=255)
                                                                                       cgcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATa                                                          <  2:4116197/150‑1 (MQ=255)
                                                                                         cAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAcc                                                        >  2:10580323/1‑150 (MQ=255)
                                                                                         cAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAcc                                                        >  1:4464335/1‑150 (MQ=255)
                                                                                           aGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCt                                                      <  2:11114066/150‑1 (MQ=255)
                                                                                              aTACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgc                                                   >  1:47078/1‑150 (MQ=255)
                                                                                               tACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgcg                                                  >  2:2137890/1‑150 (MQ=255)
                                                                                               tACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgcg                                                  >  2:9140780/1‑150 (MQ=255)
                                                                                                 ccccGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGa                                                <  2:6400965/150‑1 (MQ=255)
                                                                                                        cACCGACGGGCGGCCAGCGATGCCCCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCa                                         >  2:2629704/1‑150 (MQ=255)
                                                                                                        cACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCa                                         <  1:6083264/150‑1 (MQ=255)
                                                                                                         aCCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAc                                        <  1:3846472/150‑1 (MQ=255)
                                                                                                           cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAccc                                      >  1:6700936/1‑150 (MQ=255)
                                                                                                           cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAccc                                      >  1:48167/1‑150 (MQ=255)
                                                                                                           cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAccc                                      >  1:7617758/1‑150 (MQ=255)
                                                                                                           cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAGACCCTCGCGGAACTTTCAccc                                      >  2:1248395/1‑150 (MQ=255)
                                                                                                              cGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACccgcc                                   >  1:3596675/1‑150 (MQ=255)
                                                                                                              cGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCGCATACCCTCGCGGAACTGTCACccgcc                                   >  2:444062/1‑150 (MQ=255)
                                                                                                                 gcggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAg                                <  1:6339834/150‑1 (MQ=255)
                                                                                                                        gCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGcc                         <  1:5880553/150‑1 (MQ=255)
                                                                                                                        gCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGcc                         <  2:5739259/150‑1 (MQ=255)
                                                                                                                              ccACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTggg                   <  1:2964978/150‑1 (MQ=255)
                                                                                                                               cACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgc                  <  1:1488675/150‑1 (MQ=255)
                                                                                                                                aCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgca                 <  1:11109963/150‑1 (MQ=255)
                                                                                                                                    ccTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCAt             >  1:2486206/1‑150 (MQ=255)
                                                                                                                                         agagCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTAcc        <  2:7649700/150‑1 (MQ=255)
                                                                                                                                          gagCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCt       <  2:509833/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTC  >  NC_008463/3644212‑3644505

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: