Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,878,974 C→A 10.1% intergenic (+18/+87) rpiA → / ← PP_5151 ribose‑5‑phosphate isomerase A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,878,9740CA10.1% 913.3 / 19.5 301intergenic (+18/+87)rpiA/PP_5151ribose‑5‑phosphate isomerase A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (110/151);  new base A (37/2);  total (148/153)
Fisher's exact test for biased strand distribution p-value = 7.67e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

AACCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCC  >  NC_002947/5878825‑5879123
                                                                                                                                                     |                                                                                                                                                     
gcccaggggttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctaggtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAAc                                                                                                                                                       >  1:6277674‑M1/138‑150 (MQ=255)
aaCCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAC                                                                                                                                                       <  1:442480‑M1/150‑1 (MQ=255)
aaCCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAC                                                                                                                                                       <  1:795380‑M1/150‑1 (MQ=255)
 tccaggggttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAAcc                                                                                                                                                      <  2:6897392‑M1/14‑1 (MQ=255)
 aCCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACC                                                                                                                                                      >  2:5191978‑M1/1‑150 (MQ=255)
 aCCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACC                                                                                                                                                      >  2:1519396‑M1/1‑150 (MQ=255)
  cccGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCC                                                                                                                                                     <  2:3191102‑M1/150‑1 (MQ=255)
  cccGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCC                                                                                                                                                     <  2:304303‑M1/150‑1 (MQ=255)
   caggggttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCg                                                                                                                                                    <  2:5045641‑M1/16‑1 (MQ=255)
     ggggttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCt                                                                                                                                                  >  1:1083587‑M1/133‑150 (MQ=255)
        ggAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCACCT                                                                                                                                               >  2:7077482‑M1/1‑150 (MQ=255)
         ttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTa                                                                                                                                              >  2:6708526‑M1/129‑150 (MQ=255)
         ttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTa                                                                                                                                              >  2:3005795‑M1/129‑150 (MQ=255)
         ttgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagcagaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTa                                                                                                                                              >  2:7825433‑M1/129‑150 (MQ=255)
                 aGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCAACAAGGTAc                                                                                                                                      >  2:5643196/1‑150 (MQ=255)
                      aGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCAc                                                                                                                                 <  1:1015895/150‑1 (MQ=255)
                           aaCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTACAAGGTACCCCACCCGcc                                                                                                                            >  2:8670079/1‑150 (MQ=255)
                            aCGCCATCGTCGGCGTGGTGACCACCGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                           <  2:2099808/150‑1 (MQ=255)
                                  tcgtcgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAAcc                                                                                                                     <  2:1209064/150‑1 (MQ=255)
                                  tcgtcgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAAcc                                                                                                                     <  1:8855414/150‑1 (MQ=255)
                                    gtcgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCACCTACAACGTAACCCACCGGCCGGGACCCGt                                                                                                                   >  2:2631797/1‑150 (MQ=255)
                                      cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGCAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                 <  2:6091499/150‑1 (MQ=255)
                                      cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                 <  1:1304071/150‑1 (MQ=255)
                                       ggCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                >  1:7425192/1‑150 (MQ=255)
                                       ggCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                <  2:7425192/150‑1 (MQ=255)
                                          gtggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGAt                                                                                                             <  2:4764316/150‑1 (MQ=255)
                                            ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                           <  1:1851388/150‑1 (MQ=255)
                                               gACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtg                                                                                                        >  1:5872964‑M2/1‑148 (MQ=255)
                                               gACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCACGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtg                                                                                                        <  1:94708‑M2/150‑3 (MQ=255)
                                                 atgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTa                                                                                                      >  1:3402857‑M1/89‑150 (MQ=255)
                                                  cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtag                                                                                                     >  1:7545622‑M2/1‑145 (MQ=255)
                                                      gcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGc                                                                                                 >  1:3056155‑M1/84‑150 (MQ=255)
                                                        ggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcccgtcagggcctgTTTGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCgg                                                                                               >  1:3316318‑M1/82‑150 (MQ=255)
                                                        ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                               >  1:7803910‑M2/1‑139 (MQ=255)
                                                         gggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCggg                                                                                              >  1:151524‑M1/81‑150 (MQ=255)
                                                         cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcggg                                                                                              <  1:898482‑M1/150‑47 (MQ=255)
                                                          tGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCCGCCGGGAACCGGGGCGCTCCCgggaggcgcgggg                                                                                             <  2:4324241‑M2/150‑14 (MQ=255)
                                                            ttCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCTCCAAAAATGTACCCCACCCGCCTGGAACCCGGGCGATCCCtgtaggagcgggttt                                                                                           <  1:7536578‑M2/150‑16 (MQ=255)
                                                            accatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAg                                                                                                  <  1:9022030‑M1/66‑1 (MQ=255)
                                                                cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                       >  1:4220359‑M2/1‑131 (MQ=255)
                                                                  gcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgc                                                                                     >  1:9114111‑M2/1‑129 (MQ=255)
                                                                     ctggtcagtgagccgaacctgttctttgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCTATCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                  >  2:1913279‑M1/69‑150 (MQ=255)
                                                                       gccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACa                                                                                <  2:9021574/150‑1 (MQ=255)
                                                                       gccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACa                                                                                <  2:6211290/150‑1 (MQ=255)
                                                                          ggcggATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGc                                                                             >  2:8034799/1‑150 (MQ=255)
                                                                            cggATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCGCCTACAAGGTACCCCACCCGCCTGGCACCGTGGCGATCCCTGGAGGAGCGGGTTTAACCGCGAACACGccc                                                                           >  2:1067814/1‑150 (MQ=255)
                                                                                tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGTAACCCCCTCCAACAAGGGACCCCACCCGCCGGGGACCGGGGCGCTCCCGGGAGGGGCGGGGTTACCCGCGGACAGGcccaggc                                                                       >  1:3593556/1‑146 (MQ=255)
                                                                                 ctgctgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACa                                                                      <  2:1615180/150‑1 (MQ=255)
                                                                                   gctgctgctgGGCACCGCTTAAGGCGTGAAGAGCCTGAAGGCCGAGTAAACCCCTTCGGGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTCGGCGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                    <  1:2604883/150‑1 (MQ=255)
                                                                                    ctgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGc                                                                   <  2:8420692/150‑1 (MQ=255)
                                                                                     tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTGAAACCCCGCCGACAAGGTACCCCACCCCCCCGGGACCCTGGCGCTCCCAGGAGGGGCGGGGGTTCCCGCGCACACGCCCGGACAGGCg                                                                  >  1:4830919/1‑150 (MQ=255)
                                                                                      gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCCCCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGa                                                                 <  2:276542/150‑1 (MQ=255)
                                                                                        tgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGGAGGAGCGGGTTTACCCCCGAACAGGCCCTGACAGGCGCtc                                                               >  2:3537064/1‑150 (MQ=255)
                                                                                         gctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                              <  1:5184362/150‑1 (MQ=255)
                                                                                         gctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                              <  1:1335510/150‑1 (MQ=255)
                                                                                           tgGGCACCGCTGAAGGCGTGAAGAGCCTGACGGCCGAGTAAACCCCTGCGCGGGTAAACCCGCTCCCACCCGGTACCCCACCGGCCTGGACCCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagt                                                            <  1:3435066/150‑1 (MQ=39)
                                                                                             ggCCCCGCTGAAGGCGTGAAGAGCCTGACGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                          <  2:238675/150‑1 (MQ=255)
                                                                                              gCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCTCCAACAAGGTACCCCACCCGCCTGGAACCGGGGCGATCCCAGTAGGAGCGGGTTTACCCCCGAAAAGGCCCAGACAGGCGCtcagtcag                                                         >  2:1257918/1‑150 (MQ=255)
                                                                                               cACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATcagtcagt                                                        <  2:8740351/150‑1 (MQ=255)
                                                                                                     tGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttac                                                                                         <  2:1305607‑M2/111‑18 (MQ=255)
                                                                                                        aGGCGTGAAGAGCCTGAAGGCCGAGTAAACCCCTGCGCGGGTAAACCCGCTCCTACACGGTCCCCCACCGGCCTGGAACCGGGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCtt                                               <  1:462748/150‑1 (MQ=39)
                                                                                                        aGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCtt                                               <  1:5036315/150‑1 (MQ=255)
                                                                                                           cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCACCAACAAGCTACCCCACCCGCCGGGAACCCTGGCGATCCCAGGAGGAGCGGGGTTAACCGCGCACAGGCCCTGAAACGCGATCAGTCAGTCGGCTGCGTGaa                                            >  2:8637017/1‑150 (MQ=255)
                                                                                                            gTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGATTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAac                                           <  2:487996/150‑1 (MQ=255)
                                                                                                            gTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCAACAACGGACACCACCGGCCGGGAACCCGGGCGATACCGGTCGGGGCCGGTTTACCCGCCAACAGGCCCAGACAGGCGCTCCCTCAGTCGGCTGCTTGaaa                                           >  2:5705803/1‑149 (MQ=255)
                                                                                                              gAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                               <  2:6509400‑M2/96‑12 (MQ=255)
                                                                                                                agagCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACAt                                       <  2:8882664/150‑1 (MQ=255)
                                                                                                                 gagCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCAACTACAAGGTACCCCACCCGCCTGGGACCGTGGCGATCCCGGTAGGGGCGGGTTTACCCGCGAACAGGCCCGGACAGGCGATCAGTCAGTCCGCTGCTTGAACACATa                                      >  2:2051270/1‑150 (MQ=255)
                                                                                                                 gagCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCTCCAAAAAGGTACCCCACCCGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTAACCGCGAACAGGCCCTGACAGGGGATCAGTCAGTCGGCGGCTTGAAcacaga                                      >  2:5827062/1‑148 (MQ=255)
                                                                                                                  tgCCTGAAGGCCGTGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAg                                     <  1:7366218/149‑1 (MQ=255)
                                                                                                                  agCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAg                                     <  1:8372510/150‑1 (MQ=255)
                                                                                                                    ccTGAAGGCCGATTAAAACCCTTCGCGGTTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGaa                                   <  1:1096654/150‑1 (MQ=255)
                                                                                                                    ccTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGaa                                   <  2:2705710/150‑1 (MQ=255)
                                                                                                                       gAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAg                                <  2:7287967/150‑1 (MQ=255)
                                                                                                                       gAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCTCCAACAAGGTACCCCACCGGCCTGGAAACGTGGGGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGCTCAGTCAGTCGGCTGCTTGAACACGTAGAAAAg                                >  2:534811/1‑150 (MQ=255)
                                                                                                                        aaGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAgg                               <  2:4905135/150‑1 (MQ=255)
                                                                                                                           gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                            <  1:576975/150‑1 (MQ=255)
                                                                                                                             cGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGGCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAGCAGGCCGTCGGCGGCTTGAACACATAGGACAGGTTCgg                          >  2:8940319/1‑150 (MQ=39)
                                                                                                                             cGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                          <  1:6550406/150‑1 (MQ=255)
                                                                                                                                   aaaCCCTTGGCGGTTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTTTTGGAGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACt                    <  1:2661148/150‑1 (MQ=255)
                                                                                                                                         TTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACGGTGGCGATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAg              >  1:5384306‑M1/1‑150 (MQ=255)
                                                                                                                                          TCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAgg             >  2:3363451‑M1/1‑150 (MQ=255)
                                                                                                                                          TCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAgg             >  1:5320781‑M1/1‑150 (MQ=255)
                                                                                                                                          TCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAgg             >  1:2099198‑M1/1‑150 (MQ=255)
                                                                                                                                            gcgGGTTTACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgc           >  1:7500434‑M1/1‑21 (MQ=255)
                                                                                                                                            GCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGACCACATAGAACCGGTTCGGCTCACTGACCAGGTa           <  1:6953944‑M1/150‑1 (MQ=255)
                                                                                                                                            GCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTa           <  2:3803874‑M1/150‑1 (MQ=255)
                                                                                                                                             CGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGGAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTGCGGCTCACTGACCAGGtac          <  2:4117487‑M1/150‑2 (MQ=255)
                                                                                                                                               ggTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgc        >  2:8476108‑M1/1‑18 (MQ=255)
                                                                                                                                               GGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGACCCGGGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAAt        >  2:1931433‑M1/1‑150 (MQ=255)
                                                                                                                                                tTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCTATCCCTGTAGGTGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATg       >  1:4132196‑M1/2‑150 (MQ=255)
                                                                                                                                                gTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcg       >  1:2030149‑M1/1‑17 (MQ=255)
                                                                                                                                                GTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATg       >  1:2863467‑M1/1‑150 (MQ=255)
                                                                                                                                                 TAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATgg      >  1:2449091‑M1/1‑150 (MQ=255)
                                                                                                                                                    ACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcc   >  1:7077482‑M1/1‑150 (MQ=255)
                                                                                                                                                     cccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag  <  2:7500394‑M1/150‑139 (MQ=255)
                                                                                                                                                     cccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag  >  1:8023526‑M1/1‑12 (MQ=255)
                                                                                                                                                     CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc  <  1:6806784‑M1/150‑1 (MQ=255)
                                                                                                                                                     CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc  <  1:7176449‑M1/150‑1 (MQ=255)
                                                                                                                                                     CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc  <  1:3470093‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACCCGGTGGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCC  >  NC_002947/5878825‑5879123

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: