Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,319 T→A 5.8% 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 5.8% 916.5 / 26.3 275intergenic (+28/‑48)purB/PP_4017adenylosuccinate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (151/107);  new base A (8/9);  total (159/116)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

CAGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACG  >  NC_002947/4530170‑4530464
                                                                                                                                                     |                                                                                                                                                 
 aGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTg                                                                                                                                                  <  2:9008526/150‑1 (MQ=255)
   aaGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGc                                                                                                                                                >  1:4239164/1‑150 (MQ=255)
      cgcTGCTGACCGTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCgg                                                                                                                                             <  2:22447/150‑1 (MQ=255)
       gctgctGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAACGCCCGGCCAGGCCggg                                                                                                                                            <  2:6620942/150‑1 (MQ=255)
        ctgctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGc                                                                                                                                           <  2:8880887/150‑1 (MQ=255)
         tgctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCg                                                                                                                                          >  2:7075291/1‑150 (MQ=255)
         tgctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCg                                                                                                                                          >  1:4122621/1‑150 (MQ=255)
          gctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGt                                                                                                                                         >  2:2360988/1‑150 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGtt                                                                                                                                        >  2:8757700/1‑150 (MQ=255)
           ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGtt                                                                                                                                        <  2:6530566/150‑1 (MQ=255)
            tgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGttt                                                                                                                                       <  2:3060670/150‑1 (MQ=255)
            tgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGttt                                                                                                                                       >  1:7070362/1‑150 (MQ=255)
            tgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGttt                                                                                                                                       <  1:4707935/150‑1 (MQ=255)
               ccTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                    <  1:6679146/150‑1 (MQ=255)
                cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAt                                                                                                                                   >  2:3582938/1‑150 (MQ=255)
                cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAt                                                                                                                                   >  1:2246663/1‑150 (MQ=255)
                cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGttgtacg                                                                                                                                   <  2:1436178/150‑6 (MQ=255)
                 ttCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTAtt                                                                                                                                  >  1:5085652/1‑150 (MQ=255)
                  tCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc                                                                                                                                 <  2:7584031/150‑1 (MQ=255)
                    aTCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTccc                                                                                                                               <  2:1201356/150‑1 (MQ=255)
                     tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCg                                                                                                                              <  2:1766541/150‑1 (MQ=255)
                     tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCg                                                                                                                              >  1:5647561/1‑150 (MQ=255)
                        aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGTAc                                                                                                                           >  1:8567210/1‑150 (MQ=255)
                           gCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaa                                                                                                                        >  1:7363835/1‑150 (MQ=255)
                             cTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaaaa                                                                                                                      >  1:2754151/1‑150 (MQ=255)
                             cTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaaaa                                                                                                                      >  2:846177/1‑150 (MQ=255)
                                gACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTa                                                                                                                   <  2:5092281/150‑1 (MQ=255)
                                   aTGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGt                                                                                                                >  1:3190495/1‑150 (MQ=255)
                                      ccAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttt                                                                                                             >  2:4451738/1‑150 (MQ=255)
                                        aGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                           <  1:3430522/150‑1 (MQ=255)
                                          ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCa                                                                                                         <  2:1619686/150‑1 (MQ=255)
                                             aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGCAAATTACGTTTTTTCAATa                                                                                                      >  1:5457123/1‑150 (MQ=255)
                                             aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATa                                                                                                      <  1:1398422/150‑1 (MQ=255)
                                             aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATa                                                                                                      <  1:3172264/150‑1 (MQ=255)
                                               gCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                          >  2:7427234/1‑144 (MQ=255)
                                               gCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAAGAAAATTAAGTTTTTTCAATAgg                                                                                                    >  2:7643417/1‑150 (MQ=255)
                                                ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGt                                                                                                   <  2:2064495/150‑1 (MQ=255)
                                                ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTAAATAGGt                                                                                                   >  1:4916134/1‑150 (MQ=255)
                                                 cAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                  >  1:6715321/1‑150 (MQ=255)
                                                 cAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                  >  2:7927614/1‑150 (MQ=255)
                                                 cAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGATTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                  >  2:1331502/1‑150 (MQ=255)
                                                  aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                    >  2:6425517/1‑147 (MQ=255)
                                                  aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                    <  1:6425517/147‑1 (MQ=255)
                                                  aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                    <  1:1826715/147‑1 (MQ=255)
                                                       cGAGCTCAAGCAGCTGTCCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCCATCGCATCTAAGCTGCGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGGGCGGACCTTCCGTTTTTTCACTAGGTTGaccct                                                                                            >  2:5678993/1‑146 (MQ=255)
                                                       cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt                                                                                            <  1:3007647/150‑1 (MQ=255)
                                                       cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt                                                                                            >  1:5261975/1‑150 (MQ=255)
                                                         aGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGa                                                                                          <  2:4661670/150‑1 (MQ=255)
                                                            tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCATACGCATCTAAGCTGAGCGTCTCGTACTACGACCGGACTGGCCGTTAGTTTTTATTCCCGTAAGAATTTTACTTTTTTTCAATATTTTGAACATGAAtg                                                                                       <  2:7902532/150‑2 (MQ=255)
                                                            tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATc                                                                                       <  1:935092/150‑1 (MQ=255)
                                                              aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCt                                                                                     >  2:6025289/1‑150 (MQ=255)
                                                              aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCt                                                                                     <  1:38751/150‑1 (MQ=255)
                                                                   gcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATAc                                                                                >  2:2083481/1‑150 (MQ=255)
                                                                           ccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGc                                                                        >  2:8405598/1‑150 (MQ=255)
                                                                            cGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCa                                                                       <  1:5643241/150‑1 (MQ=255)
                                                                                      aTCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc                                                             <  2:1938848/150‑1 (MQ=255)
                                                                                       tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCg                                                                                                                              >  2:4819580/1‑84 (MQ=255)
                                                                                       tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcg                                                            >  2:2708166/1‑150 (MQ=255)
                                                                                       tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCg                                                                                                                                          <  1:4390776/72‑1 (MQ=255)
                                                                                        cGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcgg                                                           <  2:6996102/150‑1 (MQ=255)
                                                                                        cGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGGTTTTATTCCCGGACGAAAATTACGTTTTTTGAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcgg                                                           >  1:4031942/1‑150 (MQ=255)
                                                                                          gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGCAAATTACGTTTTTTCACTCGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCa                                                         >  2:8528720/1‑150 (MQ=255)
                                                                                          gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACCCCCGGCCTGGCCGGGCGTTTTTATTACCCTACGAAAAGTAAGTTTTTTAAAAAGGGTGAACATTAAACCAGATACGCCCCGGCACCAGCGCggcggcc                                                         >  1:814981/1‑149 (MQ=255)
                                                                                             aCGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTTTTGCCGTTGGAAAATTTCGTTTTTTAAATAGGTTTTTAATTAATAGTGTTAATCCACTGGATCTTATAGGCGGCAtat                                                      >  1:3258099/1‑148 (MQ=255)
                                                                                             aCGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCAtct                                                      >  2:5180919/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                    >  2:3302360/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                    <  1:2727541/150‑1 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAAGTTTTTTAAAAAGGTTGAACAATAATCCAGATACACCACAGCAGCTGCTCGGCGGCATCTCg                                                    >  1:7156765/1‑150 (MQ=255)
                                                                                               gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGGTTTTATTCCCCGACGCAAATTACGTTTTTTCAATAGGGTGTACATGAATCCAGATACTCCACAGCAGCTGCGCGGCGGCCTCTCg                                                    >  1:1132591/1‑150 (MQ=255)
                                                                                                cgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCgg                                                   <  1:1936970/150‑1 (MQ=255)
                                                                                                  gCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGTAAAATACGTTTTTTCAATAGGTTGTACATGGATCCGGGTACTCCACTGCAGCTGCTCGGCGGCATCTCGGcc                                                 >  1:1870994/1‑150 (MQ=255)
                                                                                                      cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTATTCCCGGACGAACAATACGTTTTTTCAATAGGTTGAACAAGAATCCAGATACACCACGGCACCAGCTCGGCGGCAACGCGGCCcgcc                                             >  1:914456/1‑149 (MQ=255)
                                                                                                       ccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACAGCAGCTGCTCGGCGGCATCTCGGCCCGCGa                                            >  2:5704994/1‑150 (MQ=255)
                                                                                                        cAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGaa                                           >  2:756138/1‑150 (MQ=255)
                                                                                                            ccAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCa                                       <  1:8684766/150‑1 (MQ=255)
                                                                                                             cAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCAt                                      <  1:7952130/150‑1 (MQ=255)
                                                                                                               aaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGACAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                    <  2:4119181/150‑1 (MQ=255)
                                                                                                               aaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                    <  2:4584089/150‑1 (MQ=255)
                                                                                                                 cGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGATTTTTCAATAGGATGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgcgc                                  <  2:5408815/150‑1 (MQ=255)
                                                                                                                  gCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTAAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGcgcg                                 >  2:912603/1‑150 (MQ=255)
                                                                                                                     tCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTACCGGACGAAAAATAAGTTTTTTAAAAATGTTGAACAAGAAAACAGATACACCACAGCAGCTGCTCGGCGGCAACTCGGCCCGCGAAATCATGCGCGaca                              >  1:1751663/1‑149 (MQ=255)
                                                                                                                      cTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCCTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATACAGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGacta                             >  1:528483/1‑150 (MQ=255)
                                                                                                                       tAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGactac                            >  2:4335582/1‑150 (MQ=255)
                                                                                                                       tAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGactac                            >  1:6087426/1‑150 (MQ=255)
                                                                                                                           cTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGc                        <  1:3882483/150‑1 (MQ=255)
                                                                                                                            tGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCa                       <  1:3749063/150‑1 (MQ=255)
                                                                                                                             gatcGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCag                      <  2:6269886/147‑1 (MQ=255)
                                                                                                                                   cgcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGTCCGCGAATTCATGCGCGACTACTGGCAGaagaag                <  2:7214346/150‑1 (MQ=255)
                                                                                                                                    gcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGc               <  1:7793296/150‑1 (MQ=255)
                                                                                                                                     cgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGcc              >  1:6747143/1‑150 (MQ=255)
                                                                                                                                       tACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgc            >  2:2865375/1‑150 (MQ=255)
                                                                                                                                        aCAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGaa                                           >  2:1764964/1‑118 (MQ=255)
                                                                                                                                         aaaCGCCCGGCCAGGCCGGGCGTTTTTATTACCCTACGAAAATTACGTTTTTTCAATAGGTTGAAAATTAATACTGATACACCACTGCAGCTGCTCGGCGGCATCTCGGCCCCCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctg          >  1:2563601/2‑150 (MQ=255)
                                                                                                                                          aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGACAATTACGTTTTGTCACTAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAAGTCATGCGCGACTACTGGCAGAAGAAGCCgctgc         >  2:8972863/1‑150 (MQ=255)
                                                                                                                                          aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACCTGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctgc         >  2:2983381/1‑150 (MQ=255)
                                                                                                                                          aaCGCCCGGCCTGGCCGGGCGTTTTTATTACCGTACGAAAATTACGTTTTTTGAATAGGTTGAACATGAATACTGTTACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTAATGCGCGACTCCTGGCAGAAGAAGCCgctgc         >  1:6528657/1‑150 (MQ=255)
                                                                                                                                          aaCGCCCGGCCAGGCCGGGCGTTTTAATTCCCGGAAGAAAATTACGTTTTTTCAATAGCATGAACATGAATCATGATACTCCACAGCAGATGCTCGGCGGCATCTAGGCCCGCGAAATCAAGCTCGACAACTGGCAGAAGAAGCCgctgc         >  2:5495187/1‑150 (MQ=255)
                                                                                                                                             gCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATCCTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCAAGCGCGACTACTGGCAGAAGAAGCCGCTGCTgg      >  2:8913670/1‑150 (MQ=255)
                                                                                                                                             gCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTgg      >  1:7518457/1‑150 (MQ=255)
                                                                                                                                               ccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAAGTTTTTTAAATAGGTTGTACATGTATACATATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAACCCGCTGCTGGTa    >  2:3880575/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAGAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACG  >  NC_002947/4530170‑4530464

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: