Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,319 T→A 6.2% intergenic (+28/‑48) purB → / → PP_4017 adenylosuccinate lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,530,3190TA 6.2% 500.1 / 15.0 156intergenic (+28/‑48)purB/PP_4017adenylosuccinate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (68/78);  new base A (6/4);  total (74/82)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

CAGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCCA  >  NC_002947/4530170‑4530467
                                                                                                                                                     |                                                                                                                                                    
 aGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTg                                                                                                                                                     <  2:2026730/150‑1 (MQ=255)
  gAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTgg                                                                                                                                                    <  1:3531754/150‑1 (MQ=255)
      cgcTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCgg                                                                                                                                                <  2:2351315/150‑1 (MQ=255)
          gctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGt                                                                                                                                            >  1:3061123/1‑150 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAGGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGtt                                                                                                                                           <  1:1561334/150‑1 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGtt                                                                                                                                           <  1:1561312/150‑1 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGtt                                                                                                                                           <  2:1954442/150‑1 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGtt                                                                                                                                           <  1:3656525/150‑1 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGtt                                                                                                                                           <  2:3814100/150‑1 (MQ=255)
             gACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGtttt                                                                                                                                         <  2:789964/150‑1 (MQ=255)
               ccTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                       <  1:1634776/150‑1 (MQ=255)
               ccTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGttgtac                                                                                                                                       <  1:244825/150‑5 (MQ=255)
                 ttaaTCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAtt                                                                                                                                     >  1:1223908/4‑150 (MQ=255)
                 ttCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAtt                                                                                                                                     >  1:2063781/1‑150 (MQ=255)
                 ttCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAtt                                                                                                                                     >  1:3403372/1‑150 (MQ=255)
                  tCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc                                                                                                                                    <  1:4372932/150‑1 (MQ=255)
                  tCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc                                                                                                                                    <  1:434171/150‑1 (MQ=255)
                    aTCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTccc                                                                                                                                  <  2:3572401/150‑1 (MQ=255)
                    aTCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTccc                                                                                                                                  <  1:3691161/150‑1 (MQ=255)
                        aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc                                                                                                                              >  2:1848815/1‑150 (MQ=255)
                        aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc                                                                                                                              <  1:4354917/150‑1 (MQ=255)
                                gACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTa                                                                                                                      <  2:312446/150‑1 (MQ=255)
                                 aCATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAc                                                                                                                     <  1:1411450/150‑1 (MQ=255)
                                 aCATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAc                                                                                                                     <  1:4999303/150‑1 (MQ=255)
                                     gCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAATCGCTTCTAAGCTGTGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTGCCGTTCGTAAATTAGGttt                                                                                                                 >  2:4319599/1‑150 (MQ=255)
                                     gCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGttt                                                                                                                 >  2:3788658/1‑150 (MQ=255)
                                     gCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGttt                                                                                                                 >  2:2824849/1‑150 (MQ=255)
                                     gCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGTACGAAAATTACGttt                                                                                                                 >  1:4047589/1‑150 (MQ=255)
                                       cAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGttttt                                                                                                               <  1:2836374/150‑1 (MQ=255)
                                        aGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                              <  2:748173/150‑1 (MQ=255)
                                         gCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                             <  1:2789473/150‑1 (MQ=255)
                                         gCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                             <  2:2613958/150‑1 (MQ=255)
                                          ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCa                                                                                                            >  1:5052076/1‑150 (MQ=255)
                                           cGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGttttttaca                                                                                                           >  1:3733826/1‑147 (MQ=255)
                                           cGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCaa                                                                                                           <  1:3128066/150‑1 (MQ=255)
                                           cGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAAGTTTTTTCaa                                                                                                           >  1:3733870/1‑150 (MQ=255)
                                               gCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                       >  1:4878182/1‑150 (MQ=255)
                                                ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGt                                                                                                      <  1:1496011/150‑1 (MQ=255)
                                                ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGt                                                                                                      <  1:1496069/150‑1 (MQ=255)
                                                 cAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                     >  2:1270534/1‑150 (MQ=255)
                                                  aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTg                                                                                                    <  2:1827623/150‑1 (MQ=255)
                                                  aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGTACGGAAATTAAGTTTTTTAAAAAGGGTg                                                                                                    >  1:1076628/1‑150 (MQ=255)
                                                     gCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAAc                                                                                                 <  1:4011396/150‑1 (MQ=255)
                                                            tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTGCCGTACGAAAATTAAGTTTTTTGAATAGGTTGAACATGAATc                                                                                          >  1:701629/1‑150 (MQ=255)
                                                            tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTACCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATc                                                                                          >  2:579778/1‑150 (MQ=255)
                                                              aaGCAGCTGACCCCCGCTACCTACATCTGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCt                                                                                        <  1:4826504/150‑1 (MQ=255)
                                                              aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGACCATTACGTTTTTTCAATAGGGTGAACATGAATCCt                                                                                        >  2:3928725/1‑150 (MQ=255)
                                                              aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCt                                                                                        >  1:2646047/1‑150 (MQ=255)
                                                               agcagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTg                                                                                       <  1:4260732/150‑1 (MQ=255)
                                                                gcagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa                                                                                      <  2:4593798/150‑1 (MQ=255)
                                                                 cagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGAt                                                                                     <  1:194740/150‑1 (MQ=255)
                                                                    cTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAACATTACGTTTTTTCAATAGGGTGAACATGAATCCTGATACt                                                                                  >  2:4923189/1‑150 (MQ=255)
                                                                    cTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATaat                                                                                  >  2:687435/1‑148 (MQ=255)
                                                                    cTGACCCCCGCTACCTACATCGGTAAAGCGGCAGCCCAGGCCAAACGCATCTAAGATGAGCGTCGCGTACAAAGCCCGGCCTGGACGGGCGTGTTTATTCCCGTAAGCAACTTAACTTTTTTAAAAAGGTGGAAAAAGTATCAAGATAca                                                                                  >  1:913208/1‑149 (MQ=255)
                                                                          cccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTg                                                                            <  1:1787289/150‑1 (MQ=255)
                                                                              ctacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGGACATGAATCCTGATACTCCACTGCAgc                                                                        >  1:412066/1‑150 (MQ=255)
                                                                               tacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgct                                                                       <  2:895099/150‑1 (MQ=255)
                                                                               tacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGTACGAAAATTACGTTTTTTGAAAAGGTTGAACATGAATCCTGATACTCAACTGCAgct                                                                       >  1:1559385/1‑150 (MQ=255)
                                                                                      aTCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACGGCAGCTGCTcggc                                                                >  1:4031900/1‑150 (MQ=255)
                                                                                          gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCa                                                            <  1:3232884/150‑1 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGGAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                       >  2:2502717/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                       >  1:3372578/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                       <  1:2560345/150‑1 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                       >  2:2502960/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                       >  1:3372549/1‑150 (MQ=255)
                                                                                                cgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTTTTCCCGTACGAAAATTACGTTTTTTAAATAGGTTGAAAATGAATCCTGATACTCCACGGCAGCTGCTCGGCGGCATCTCgg                                                      >  1:4753256/1‑150 (MQ=255)
                                                                                                     gCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACACCGCCCGGCCTGGTCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCGGCTCGGCGGCATCTCGGCCcgc                                                 <  2:2552318/150‑1 (MQ=255)
                                                                                                      cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgcg                                                >  2:4422785/1‑150 (MQ=255)
                                                                                                      cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgcg                                                >  2:3503918/1‑150 (MQ=255)
                                                                                                         aGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAAGAAAATTAAGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCCGCTGCTCGGCGGCATCTCGGCCCGCGact                                             >  1:1865544/1‑148 (MQ=255)
                                                                                                           gCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTc                                           >  2:3227393/1‑150 (MQ=255)
                                                                                                            ccAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCa                                          <  2:4650816/150‑1 (MQ=255)
                                                                                                             cAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCAt                                         <  2:4992326/150‑1 (MQ=255)
                                                                                                             cAAACGCAAATAAGCAGAGCGTCGCGTACAACGCAAGGCCTGGCCGGGCGTTTATATTCCCGTACGAAAATTAAGTTTTTTCAATAGGTTGAACATGAAAACTGATAATCCACAGCAGCTTCTCGGAGGCATATCGGCCCGCGAAATCAt                                         >  2:5103054/1‑150 (MQ=255)
                                                                                                               aaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                       <  1:3246350/150‑1 (MQ=255)
                                                                                                               aaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                       <  2:4143797/150‑1 (MQ=255)
                                                                                                               aaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                       <  1:2918250/150‑1 (MQ=255)
                                                                                                                  gCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGcgcg                                    >  2:1985016/1‑150 (MQ=255)
                                                                                                                  gCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAAAAGGTTTAACATGAATACTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGcgcg                                    >  2:954982/1‑150 (MQ=255)
                                                                                                                    aTCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGac                                  >  2:912714/1‑150 (MQ=255)
                                                                                                                    aTCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAAGAAAATTACGTTTTTTCTATAGGTTGAACATTAATCCTGTTACTCCACTGGAGCTGCTCGGCGGGATCTCGGCCCGGGAATTCATGCgagag                                  >  2:3007538/1‑149 (MQ=255)
                                                                                                                     tCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGCATGAATTcgcgggc                                 <  2:4723807/150‑4 (MQ=255)
                                                                                                                          gCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTgg                            <  2:48206/150‑1 (MQ=255)
                                                                                                                          gCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAAGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGAAACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTgg                            <  2:2022357/150‑1 (MQ=255)
                                                                                                                                   cgcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGaagaag                   <  1:991931/150‑1 (MQ=255)
                                                                                                                                    gcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGc                  >  2:1175018/1‑150 (MQ=255)
                                                                                                                                    gcgTACAAAGCCCGGCCTGGCCGGGCGTGTTTATTAACGTACGTAATTTACGTTTTTTAAATAGGTTTAACAAGAATACTTATACTCCACTGCAGCTGCTCGGCGGCATCACGGCCCGCGAATTAATGAGCGTCTACTGGGAGtataaga                  >  1:1397719/1‑149 (MQ=255)
                                                                                                                                      taaaaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCg                >  2:1235806/5‑150 (MQ=255)
                                                                                                                                      taaaaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCg                >  2:1235624/5‑150 (MQ=255)
                                                                                                                                      gTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGTAAATTACGTTTTTTCAATTGGTTGAACATGTATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGCATTCATGCGCGACTACTGGCAGAAGAAGCCg                >  2:1344050/1‑150 (MQ=255)
                                                                                                                                       tACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgc               <  2:4036318/150‑1 (MQ=255)
                                                                                                                                          aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGCAAATTACGTTTTTTCAATAGGTTGAACCTGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATATCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctgc            >  2:1719114/1‑150 (MQ=255)
                                                                                                                                          aaCGCCCGGCCAGGCCGGGCGTTTTTATTACCGGACGAGAATTACGTTTTTTCAATAGGTTGTACATGAATCCGGAGACGCCACAGCACCTGCGCGGCGGCCTCTCGGCCCGCCAATTCATGCGCGACAAATGGCAGAAGAATCCgcggc            >  1:3580379/1‑150 (MQ=255)
                                                                                                                                            cGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGctgctg          <  2:1386554/150‑1 (MQ=255)
                                                                                                                                            cGCCCGGCCTGGCCGGGCGATTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGctgctg          <  2:3403364/150‑1 (MQ=255)
                                                                                                                                            cGCCCGGCCTGGCCGGGCGATTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGctgctg          <  2:3403372/150‑1 (MQ=255)
                                                                                                                                             gCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAAAAGGTTGAACATGAATCCTGATACTCCACTGCAGCATCTCGGCGGCATCACGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTgg         >  2:3507842/1‑150 (MQ=255)
                                                                                                                                               ccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTa       <  1:2213746/150‑1 (MQ=255)
                                                                                                                                                  gCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGc    <  1:4423994/150‑1 (MQ=255)
                                                                                                                                                    cTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCCa  <  2:277757/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCCA  >  NC_002947/4530170‑4530467

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: