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% 2562.9 / ‑3.5 517intergenic (‑112/‑123)sohB/PA14_RS16560protease SohB/histidine phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (234/282);  minor base T (1/0);  total (235/282)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 546 total aligned reads displayed.

CCTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCGA  >  NC_008463/3644211‑3644507
                                                                                                                                                      |                                                                                                                                                  
ccTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAg                                                                                                                                                    >  2:3091002/1‑151 (MQ=255)
 cTGGCCAGACAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgg                                                                                                                                                   <  2:11037600/151‑1 (MQ=255)
  tGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAggg                                                                                                                                                  >  1:780001/1‑151 (MQ=255)
   ggCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggc                                                                                                                                                 <  1:1938508/151‑2 (MQ=255)
    gCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                                >  2:7938073/1‑149 (MQ=255)
       aGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGc                                                                                                                                            <  2:1207149/151‑1 (MQ=255)
             cccGCATATTCGGCAATGAACTCCACACACCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGag                                                                                                                                      <  1:395427/151‑1 (MQ=255)
             cccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCACGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGag                                                                                                                                      >  2:1452804/1‑151 (MQ=255)
             cccGCATAGTCGGCAATGAACTCCACCCCCCCCCCTCTCTCCGGGAGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCCGGCCACCGCCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGag                                                                                                                                      <  2:7613483/151‑1 (MQ=255)
              ccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagc                                                                                                                                     <  2:10320015/151‑1 (MQ=255)
              ccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagc                                                                                                                                     >  2:2431168/1‑151 (MQ=255)
              ccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGACTACGGGCGACACGAAAGCGCAAGCATACCCAGTCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGATAGCCTTACCAGGGCACGACCGCagc                                                                                                                                     >  1:551533/1‑151 (MQ=255)
               cGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagca                                                                                                                                    <  1:10046054/151‑1 (MQ=255)
                gCATCTTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTATACGGGCGACACGATTGCGCTAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagcag                                                                                                                                   <  2:4629189/151‑1 (MQ=255)
                  atatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGcc                                                                                                                                 <  1:10170960/151‑1 (MQ=255)
                   tatTCGGCAATGAACTCCACACACCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCg                                                                                                                                <  2:5658653/151‑1 (MQ=255)
                   tatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCg                                                                                                                                <  1:5959178/151‑1 (MQ=255)
                     ttCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAt                                                                                                                              <  1:5796596/151‑1 (MQ=255)
                       cGGCAATGAACTCCACCCAACCCCCTCTCTACTTTTGACGCACCCGTCTACGGGCGACAAGATTGCGCACGCATAACCAGGCCAGCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCATGGCACGCCCTGAGCAGCCGATaa                                                                                                                            <  2:2328268/151‑1 (MQ=255)
                       cGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                            <  1:2411839/151‑1 (MQ=255)
                          cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                         <  2:4580265/151‑1 (MQ=255)
                          cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                         >  1:4939790/1‑151 (MQ=255)
                             tGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAg                                                                                                                      >  1:4781030/1‑151 (MQ=255)
                               aaCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCTGCg                                                                                                                    >  2:5223776/1‑151 (MQ=255)
                               aaCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                    <  2:303711/151‑1 (MQ=255)
                               aaCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                    >  2:6468022/1‑151 (MQ=255)
                               aaCTCCACACAACCCCCCCCCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                    <  2:10659367/151‑1 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  <  2:4157122/151‑1 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  >  2:434762/1‑151 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  >  1:7253735/1‑151 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  >  1:3744200/1‑151 (MQ=255)
                                   ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATg                                                                                                                <  2:3448975/151‑1 (MQ=255)
                                     acacaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGAc                                                                                                              >  2:1753116/1‑151 (MQ=255)
                                       acaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAt                                                                                                            >  1:5300036/1‑151 (MQ=255)
                                         aaCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTc                                                                                                          <  1:6073519/151‑1 (MQ=255)
                                          aCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAATGGCAGCGAATGACATTCa                                                                                                         <  2:2055901/151‑1 (MQ=255)
                                           cccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaa                                                                                                        >  2:3865680/1‑151 (MQ=255)
                                            ccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaaa                                                                                                       >  1:3257958/1‑151 (MQ=255)
                                            ccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaaa                                                                                                       >  2:7385224/1‑151 (MQ=255)
                                            ccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaaa                                                                                                       <  1:10315994/151‑1 (MQ=255)
                                             cccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAAc                                                                                                      >  2:9844010/1‑151 (MQ=255)
                                             cccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAAc                                                                                                      >  1:2812737/1‑151 (MQ=255)
                                              ccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                     >  2:9502138/1‑151 (MQ=255)
                                              ccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                     <  2:6195860/151‑1 (MQ=255)
                                               ctctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACaa                                                                                                    >  1:867265/1‑151 (MQ=255)
                                                  tcgaCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCg                                                                                                 <  2:7055419/148‑1 (MQ=255)
                                                     aCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAt                                                                                              <  1:7544713/151‑1 (MQ=255)
                                                     aCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAt                                                                                              <  2:6739418/151‑1 (MQ=255)
                                                       tGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGt                                                                                            <  1:9598835/151‑1 (MQ=255)
                                                        gTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtt                                                                                           >  2:7541413/1‑151 (MQ=255)
                                                          tGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttt                                                                                         <  1:11209303/151‑1 (MQ=255)
                                                            aCGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                       >  1:4473403/1‑151 (MQ=255)
                                                             cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTc                                                                                      <  2:4921739/151‑1 (MQ=255)
                                                              gCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCg                                                                                     <  1:5090947/151‑1 (MQ=255)
                                                              gCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCg                                                                                     <  2:9766564/151‑1 (MQ=255)
                                                              gCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAGAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCg                                                                                     >  2:10564670/1‑151 (MQ=255)
                                                                aCCCGTCTACGGGCGCCACGATTGCGCAAGCATACCCAGGCCACCGCCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGATTTTTCGtt                                                                                   <  2:4457359/151‑1 (MQ=255)
                                                                aCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGtt                                                                                   >  1:8439125/1‑151 (MQ=255)
                                                                aCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGtt                                                                                   >  1:10565335/1‑151 (MQ=255)
                                                                 cccGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTg                                                                                  <  1:3406372/151‑1 (MQ=255)
                                                                  ccGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGa                                                                                 <  2:6900491/151‑1 (MQ=255)
                                                                    gTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                               <  2:1521258/151‑1 (MQ=255)
                                                                      cTACGGGCGACACGATTGCGCACGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGc                                                                             <  2:10180374/151‑1 (MQ=255)
                                                                      cTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGc                                                                             <  1:2951562/151‑1 (MQ=255)
                                                                       tACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGcc                                                                            >  1:9469291/1‑151 (MQ=255)
                                                                          gggCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTg                                                                         <  1:2159822/151‑1 (MQ=255)
                                                                           ggCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTgg                                                                        <  2:8213412/151‑1 (MQ=255)
                                                                           ggCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTgg                                                                        >  1:6211841/1‑151 (MQ=255)
                                                                             cGACACGAGTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGt                                                                      <  2:8060405/151‑1 (MQ=255)
                                                                                    aTTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTAGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTc                                                               <  2:554239/151‑1 (MQ=255)
                                                                                      tGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAt                                                             <  1:5745857/151‑1 (MQ=255)
                                                                                        cgcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAc                                                           <  2:1989623/151‑1 (MQ=255)
                                                                                            aGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTc                                                       <  2:7264506/151‑1 (MQ=255)
                                                                                              cATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgc                                                     <  1:7187019/151‑1 (MQ=255)
                                                                                               aTACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgcg                                                    >  1:7988669/1‑151 (MQ=255)
                                                                                                tACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCgg                                                   <  2:878284/151‑1 (MQ=255)
                                                                                                 aCCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGa                                                  >  1:4057849/1‑151 (MQ=255)
                                                                                                      ggCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACttt                                             <  2:5823490/151‑1 (MQ=255)
                                                                                                       gCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTc                                            <  2:3353307/151‑1 (MQ=255)
                                                                                                       gCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTc                                            <  1:1304714/151‑1 (MQ=255)
                                                                                                        ccACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCa                                           <  1:6534416/151‑1 (MQ=255)
                                                                                                         cACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAc                                          >  1:9785026/1‑151 (MQ=255)
                                                                                                           ccGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAccc                                        <  2:10091459/151‑1 (MQ=255)
                                                                                                            cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACccg                                       <  2:2266147/151‑1 (MQ=255)
                                                                                                                  gcggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAgg                                 >  2:9215891/1‑151 (MQ=255)
                                                                                                                    ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTc                               >  2:2739742/1‑151 (MQ=255)
                                                                                                                    ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTc                               >  1:1598294/1‑151 (MQ=255)
                                                                                                                    ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTc                               <  1:2604282/151‑1 (MQ=255)
                                                                                                                     gcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTcc                              <  2:2587996/151‑1 (MQ=255)
                                                                                                                       cAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGc                            >  2:8571939/1‑151 (MQ=255)
                                                                                                                         gCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGccc                          <  1:1105230/151‑1 (MQ=255)
                                                                                                                               ccACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgc                    <  2:6385053/151‑1 (MQ=255)
                                                                                                                                 aCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgcag                  >  1:9272846/1‑151 (MQ=255)
                                                                                                                                  cTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgcagc                 <  2:137622/151‑1 (MQ=255)
                                                                                                                                   tGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgcagca                >  2:3573961/1‑151 (MQ=255)
                                                                                                                                     ccTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATc              >  1:5667660/1‑151 (MQ=255)
                                                                                                                                      cTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCt             >  1:4126441/1‑151 (MQ=255)
                                                                                                                                         gagaGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTAcc          <  1:7004398/151‑1 (MQ=255)
                                                                                                                                         gagaGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTAcc          >  2:5992881/1‑151 (MQ=255)
                                                                                                                                          agagCCTTTTCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCt         <  1:3512488/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCGA  >  NC_008463/3644211‑3644507

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: