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% 2424.4 / ‑0.9 492intergenic (‑112/‑123)sohB/PA14_RS16560protease SohB/histidine phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (230/260);  minor base A (1/0);  total (231/261)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 526 total aligned reads displayed.

CCTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCGA  >  NC_008463/3644211‑3644507
                                                                                                                                                      |                                                                                                                                                  
 cTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgg                                                                                                                                                   <  2:6732529/151‑1 (MQ=255)
 cTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgg                                                                                                                                                   >  1:1581751/1‑151 (MQ=255)
   ggCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggc                                                                                                                                                 >  2:10391243/1‑150 (MQ=255)
    gCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                                >  1:3952267/1‑149 (MQ=255)
      cAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggcacg                                                                                                                                              <  1:8782793/151‑5 (MQ=255)
      cAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggcacg                                                                                                                                              <  2:6982066/151‑5 (MQ=255)
         aaaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGccc                                                                                                                                          >  2:8883012/1‑151 (MQ=255)
         aaaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGccc                                                                                                                                          <  1:3443075/151‑1 (MQ=255)
         aaaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGccc                                                                                                                                          <  2:3915904/151‑1 (MQ=255)
         aaaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGccc                                                                                                                                          <  2:2780049/151‑1 (MQ=255)
          aaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCg                                                                                                                                         >  2:3018397/1‑151 (MQ=255)
             cccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGag                                                                                                                                      >  2:5542421/1‑151 (MQ=255)
               cGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagca                                                                                                                                    <  1:660779/151‑1 (MQ=255)
                gCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagcag                                                                                                                                   >  1:2813739/1‑151 (MQ=255)
                  atatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGcc                                                                                                                                 >  2:8004146/1‑151 (MQ=255)
                  atatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGcc                                                                                                                                 <  2:9376004/151‑1 (MQ=255)
                    atTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGa                                                                                                                               >  2:1994927/1‑151 (MQ=255)
                     ttCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAt                                                                                                                              <  2:6575898/151‑1 (MQ=255)
                      tCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATa                                                                                                                             >  1:7565876/1‑151 (MQ=255)
                        cgCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaaa                                                                                                                           <  2:2197698/150‑1 (MQ=255)
                         gCAATGAACTCCACACACCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAg                                                                                                                          <  2:5878027/151‑1 (MQ=255)
                          cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                         <  2:4617776/151‑1 (MQ=255)
                          cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                         <  1:2084986/151‑1 (MQ=255)
                           aaTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGc                                                                                                                        >  1:10260562/1‑151 (MQ=255)
                            aTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCa                                                                                                                       <  2:5015843/151‑1 (MQ=255)
                             tGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAg                                                                                                                      >  1:6125130/1‑151 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGTCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  >  1:4560604/1‑151 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  >  1:10005747/1‑151 (MQ=255)
                                 cTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGaa                                                                                                                  <  1:4726580/151‑1 (MQ=255)
                                  tCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                                 >  2:6336771/1‑151 (MQ=255)
                                   ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATg                                                                                                                >  2:7388464/1‑151 (MQ=255)
                                   ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATg                                                                                                                <  1:10756318/151‑1 (MQ=255)
                                    cacaccCCCCCCTCTCGACTGCTGACGCCCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACAGACGGGCTGACAGCTATGCCACGTCCTGGAGAGCCTTTGCATTGCACGCCCGGAGCAGCCGTTAAAGGCAGCGAATGa                                                                                                               <  2:5656272/151‑1 (MQ=255)
                                     acaccACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGAc                                                                                                              <  1:290171/151‑1 (MQ=255)
                                       acaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAt                                                                                                            <  1:1440639/151‑1 (MQ=255)
                                        caACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAtt                                                                                                           <  1:7259264/151‑1 (MQ=255)
                                         aaCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTc                                                                                                          >  2:31397/1‑151 (MQ=255)
                                           cccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaa                                                                                                        >  2:482384/1‑151 (MQ=255)
                                            ccccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCaaa                                                                                                       >  1:2464194/1‑151 (MQ=255)
                                              ccTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACa                                                                                                     <  1:8527699/151‑1 (MQ=255)
                                               ctctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACaa                                                                                                    >  2:8720274/1‑151 (MQ=255)
                                                tctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAg                                                                                                   >  1:8092831/1‑151 (MQ=255)
                                                 ctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGc                                                                                                  >  2:811834/1‑151 (MQ=255)
                                                  tctACTGATGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGACCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAATCg                                                                                                 >  2:7201467/1‑151 (MQ=255)
                                                   ctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGTCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATATAGGCAGCGAATGACATTCAAACAAGCGt                                                                                                <  2:545604/151‑1 (MQ=255)
                                                   ctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGt                                                                                                >  2:8364960/1‑151 (MQ=255)
                                                     aCTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAt                                                                                              <  1:329262/151‑1 (MQ=255)
                                                      cTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATg                                                                                             >  1:6810786/1‑151 (MQ=255)
                                                       tGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGACTGACATTCAAACAAGCGGATGt                                                                                            >  2:10497610/1‑151 (MQ=255)
                                                       tGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGt                                                                                            <  2:6911996/151‑1 (MQ=255)
                                                           gACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGttttt                                                                                        <  1:9881849/151‑1 (MQ=255)
                                                                aCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGtt                                                                                   <  1:644541/151‑1 (MQ=255)
                                                                aCCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGtt                                                                                   >  2:4265594/1‑151 (MQ=255)
                                                                  ccGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGa                                                                                 >  2:7740703/1‑151 (MQ=255)
                                                                   cGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACACGCGTATGTTTTTTCGTTGAc                                                                                <  2:1216434/151‑1 (MQ=255)
                                                                   cGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGAc                                                                                >  2:10699850/1‑151 (MQ=255)
                                                                    gTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                               <  1:1166929/151‑1 (MQ=255)
                                                                    gTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACa                                                                               <  2:9521443/151‑1 (MQ=255)
                                                                     tCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAg                                                                              >  1:9392630/1‑151 (MQ=255)
                                                                      cTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGc                                                                             >  1:1031604/1‑151 (MQ=255)
                                                                             cGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGt                                                                      >  1:7823846/1‑151 (MQ=255)
                                                                                cacGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGAt                                                                   <  1:4232608/151‑1 (MQ=255)
                                                                                cacGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGAt                                                                   >  2:6138260/1‑151 (MQ=255)
                                                                                 acGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATc                                                                  <  2:5358978/151‑1 (MQ=255)
                                                                                  cGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATcc                                                                 >  1:9980298/1‑151 (MQ=255)
                                                                                    agtGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTGTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTGGTTCGTTGACAGCCCTGGGTGATCCTc                                                               >  2:823259/3‑151 (MQ=255)
                                                                                        cgcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAc                                                           <  1:1318489/151‑1 (MQ=255)
                                                                                        cgcAAGCAGACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAc                                                           <  1:56098/151‑1 (MQ=255)
                                                                                          cAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAccc                                                         <  2:3567031/151‑1 (MQ=255)
                                                                                             gCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcg                                                      <  1:9507400/151‑1 (MQ=255)
                                                                                               aTACCCCGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgcg                                                    <  1:2222886/151‑1 (MQ=255)
                                                                                                     aGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACtt                                              >  1:4381768/1‑151 (MQ=255)
                                                                                                       gCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTc                                            >  1:9161907/1‑151 (MQ=255)
                                                                                                        ccACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTAGGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCa                                           >  2:1920633/1‑151 (MQ=255)
                                                                                                              aCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACccgcc                                     >  2:1424372/1‑151 (MQ=255)
                                                                                                                   cggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGt                                <  1:5299408/151‑1 (MQ=255)
                                                                                                                    ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTc                               >  2:1551764/1‑151 (MQ=255)
                                                                                                                     gcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTcc                              <  1:8631799/151‑1 (MQ=255)
                                                                                                                       cAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGc                            >  2:10109157/1‑151 (MQ=255)
                                                                                                                         gCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGccc                          <  2:1719532/151‑1 (MQ=255)
                                                                                                                            aTGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTg                       >  1:5992264/1‑151 (MQ=255)
                                                                                                                              gCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTggg                     >  2:8510359/1‑151 (MQ=255)
                                                                                                                               ccACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATCCCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgc                    <  2:7593443/151‑1 (MQ=255)
                                                                                                                                cACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgca                   >  1:8308828/1‑151 (MQ=255)
                                                                                                                                cACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgca                   <  1:4884021/151‑1 (MQ=255)
                                                                                                                                 aCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgcag                  >  1:909301/1‑151 (MQ=255)
                                                                                                                                    gCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCAt               >  1:10470961/1‑151 (MQ=255)
                                                                                                                                       tGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTa            <  1:700734/151‑1 (MQ=255)
                                                                                                                                       tGGAGAGCCTGTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTa            <  2:6056390/151‑1 (MQ=255)
                                                                                                                                        ggAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTAc           <  1:1762602/151‑1 (MQ=255)
                                                                                                                                        ggAGAGCCTTGGCATGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAGGCCATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGAGCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTAc           <  2:4446498/151‑1 (MQ=255)
                                                                                                                                          agagCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTTACAGCACTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCt         >  2:9606282/1‑151 (MQ=255)
                                                                                                                                          agagCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTACGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCt         <  2:2534144/151‑1 (MQ=255)
                                                                                                                                           gagCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACACGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTg        <  1:6339065/151‑1 (MQ=255)
                                                                                                                                             gCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGAt      <  2:2391298/151‑1 (MQ=255)
                                                                                                                                              ccTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATTACATTCAAACAAGCGAATGTATTAACGTTGACAGCCCTGGCTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGACCGCCCGTGGGCAGCAACTACCAGAtt     >  1:3516134/1‑151 (MQ=255)
                                                                                                                                              ccTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGAtt     >  2:3421454/1‑151 (MQ=255)
                                                                                                                                              ccTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGAtt     >  1:5405255/1‑151 (MQ=255)
                                                                                                                                               cTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTc    >  2:3774760/1‑151 (MQ=255)
                                                                                                                                                tttGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCg   <  1:6221827/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: