Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,275,9570AG25.0% 110.3 / 10.2 53S163G (AGT→GGT) pdeCc‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base A (8/31);  new base G (13/0);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGCGAATGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTT  >  NC_000913/4275710‑4276199
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aaaGCGAATGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATTAAC‑CCCCTCTTCTGGtgtg                                                                                                                                                                                                                                                   >  1:303918/1‑251 (MQ=255)
       aTGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTTTTAAC‑CCTCTCTTCTGGGGTGAAGTGAt                                                                                                                                                                                                                                            >  1:620335/1‑251 (MQ=255)
               tattatTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGATCTATATGCAGCGGGGAAATTATGTGGCGGTTTTCAAA‑CCTCCATTCTGGGGTGAAGTGGTGTCttgtg                                                                                                                                                                                                                                    >  1:124976/1‑247 (MQ=255)
                  tattCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTT‑ATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTTTCAAC‑CCTCTCTTCTGGGGGGAAGTGGTGTCTTATGAccc                                                                                                                                                                                                                                >  1:306042/1‑251 (MQ=255)
                   attCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTTGGGTGAAGGGATGttttgtgaccc                                                                                                                                                                                                                                >  1:210211/1‑241 (MQ=255)
                       gTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTATCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACa                                                                                                                                                                                                                            <  1:467372/251‑1 (MQ=255)
                                         aTGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATTAAC‑CCTCTCTTCTGGAGGGAAGGGATGTCTGATGGCCCGCCATTGCAATGGGgggtgtgt                                                                                                                                                                                                          >  1:339442/1‑249 (MQ=255)
                                         aTGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGGTTTCACA‑CCTCTCTTTTGGGGGGTAGGGGTGTCTGGTGTCCCGCCACTGCAATGGGgggtgtcg                                                                                                                                                                                                          >  2:82424/1‑249 (MQ=255)
                                           gAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATg                                                                                                                                                                                                         <  2:601137/250‑1 (MQ=255)
                                                gATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTTTTAAC‑CCTCTCTTCTGGGGGGAAGGGGTGTCTGTTGACCCGACATTGCAATGGGGTGTGGATGATACgg                                                                                                                                                                                                   >  1:20195/1‑251 (MQ=255)
                                                     aTGCCCGTTATAACCATTTTTTATGATCAGCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACgg                                                                                                                                                                                                   <  2:451514/246‑1 (MQ=255)
                                                       gCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGaa                                                                                                                                                                                           <  2:66869/251‑1 (MQ=255)
                                                                aaCCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTTTCAAA‑CCTCTCTTCTGGAGTGGAGTGGTGTCTGGTGACCCCGCATTGCAATGGGGGGTGTATGATAAGGTG‑ACGAAAAACttttt                                                                                                                                                                                  >  1:527152/1‑251 (MQ=255)
                                                                             gCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCCCCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCaa                                                                                                                                                                      <  1:598855/250‑1 (MQ=255)
                                                                                cATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAagag                                                                                                                                                                  <  1:300774/251‑1 (MQ=255)
                                                                                 aTCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGGGTGAAGGGATGTCTGGTGACCCCACATTGCAATGGGGGGGGGATGATACGGGG‑ACGAAAACCTTTTTTTCGTTAAGCAAAGAgg                                                                                                                                                                 >  1:60625/1‑251 (MQ=255)
                                                                                    gCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCt                                                                                                                                                              <  1:412299/251‑1 (MQ=255)
                                                                                         tAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCCCTCTTCTGGGGTGAAGGGATGGCTGATGACCCCACATTGCCATGGGGGGTGTTTGATACGGGG‑ACGAAAAACTTTTTTTCGTTTAGCaaaaaaggctcggcc                                                                                                                                                         >  1:39216/1‑239 (MQ=255)
                                                                                             gcCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTTGGGTGGAGGGGTGTCTGATGACCCGACCTTGCACTGGGGGGGGTTTTATAAGGTGCAGGAAAAC‑TTTTTTTCGTTTAGCAAAGAGGCCCCTgcagcc                                                                                                                                                      >  1:379568/1‑249 (MQ=255)
                                                                                                    aaCGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTtc                                                                                                                                              <  2:303918/251‑1 (MQ=255)
                                                                                                      cGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTTACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CATCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTtctc                                                                                                                                            <  2:399204/251‑1 (MQ=255)
                                                                                                        gCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATGATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACATATATGCAGCGTGGAAAGTATGTGGCGTTTATCAAC‑CCTCGCTTCTGGATTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTAGGATACGGTG‑ATGAAAGCCGTTTACTCGTTAAGCAAATAGGCCACGGCAGCAACGTATtctcag                                                                                                                                          <  2:225190/251‑3 (MQ=255)
                                                                                                           ataCGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCt                                                                                                                                        <  2:60625/250‑1 (MQ=255)
                                                                                                              cGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTCTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGGGGGAAGGGATGTCTGGTGACCCGACATTGCAATGGGGTGTGTATGATACGGGG‑GCGAAAACCTTTTTTTCGTTAAGCAAAAAGGCCTCGGCAGCAAAGTTTTCTCCGCTggtt                                                                                                                                    >  1:539068/1‑248 (MQ=255)
                                                                                                               gATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGGGGGAAGTGATGTCTGATGAACCGACATTGCAATGGGGGGTGTATGATACGGGG‑ACGAAAACCCTTTTTTCGTTAAGCAAAGAGGGCCCGGCAGCAACGTTTTCTCCGCCGATTc                                                                                                                                   >  1:61733/1‑251 (MQ=255)
                                                                                                                          ccGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCGCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGAt                                                                                                                        <  1:57908/251‑1 (MQ=255)
                                                                                                                                ataAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGGTACGCCTTTTGTCTCGGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGGGGCGGTTATCAAC‑CCTGTGGTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTGTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATGCACTTGAAGGATTTAAcc                                                                                                                  <  2:379568/251‑1 (MQ=255)
                                                                                                                                ataAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAAc                                                                                                                   <  1:432912/250‑1 (MQ=255)
                                                                                                                                      gTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACaa                                                                                                            <  1:625140/251‑1 (MQ=255)
                                                                                                                                                cGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCt                                                                                                  <  2:339442/251‑1 (MQ=255)
                                                                                                                                                   ttCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAt                                                                                                <  1:206089/250‑1 (MQ=255)
                                                                                                                                                       atctatTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtat                                                                                            <  2:20195/250‑1 (MQ=255)
                                                                                                                                                       atctatTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCCCTCTTTTGGAGGGAAGTGATGTCTGGTGAACCGCCATTGCAATGGGGGGGGTATGATACGGTG‑ACGCAAAACTTTTTTTCGTTAAGCAAAAAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAAAGGCTATTtatt                                                                                           >  1:510779/1‑250 (MQ=255)
                                                                                                                                                                 aCCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATTAAA‑CCTCTCTTCTTGGGTGAAGGGGTGTCTGATGACCCGACATTGCCATGGGGGGGGTGTGATACGGTG‑ACGAAAACCTTTTTTTCGTTTAGCAAAGAGGGCTCGGCAGCAACGTTTTTTCCGCTGATTCTTTTGTAAGATTTTACCGCACAAAGAAAAGGCTATTTttatacggcaggt                                                                                 >  2:168473/1‑242 (MQ=255)
                                                                                                                                                                   cgcgATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTa                                                                               <  1:248095/251‑1 (MQ=255)
                                                                                                                                                                    gcgATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTAT‑AACACCCCTCTTCTGGAGGGGAGAGATGTCTGGTGAACCGCCATTGCAATGGGGGGGGGATGATAAGGGG‑ACGAAAAACTTTTTTTCGTTAAGCAAAGAGGCCTCGGCAACAAAGTTTTCTTCGCCGATTTATTTGAAGGGTTTAACCGTACAAAGAAAAGGGTATTTTTATTTGGCAGTTTAt                                                                              >  1:517113/1‑251 (MQ=255)
                                                                                                                                                                      gATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCCCTCTTTTGGGGGGAAGTGATGTCTGATGACCCGCCCTTGCAATGGGGGGTGTATGATACGGGG‑GCGAAAACCCTTTTTTCGTTTAGAAAAACGGGCCCGGCAGCAACGATTTCTCCGCTGATTTATTTGAAGGATTTTACCCCACAAAAAAAACGCTATTTTTATGCCACAGTTTAttt                                                                            >  1:290158/1‑250 (MQ=255)
                                                                                                                                                                           gCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACaa                                                                       <  2:527152/251‑1 (MQ=255)
                                                                                                                                                                                        ggCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTg                                                           <  1:232571/250‑1 (MQ=255)
                                                                                                                                                                                          cTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGGGAAGGGATGTCTGATGACCCGGCATTGCAATTGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAAAGGCCTCGGCAGCAACGTTTTCTTCGCTGTTTCATTTGGAGGATTTTACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAAAGCCCAATTGc                                                          >  1:280360/1‑249 (MQ=255)
                                                                                                                                                                                              aaaaTGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAt                                                    <  1:221789/251‑1 (MQ=255)
                                                                                                                                                                                               aaaTGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGGGGTGTCTGATGACCCGACATTGTCATGGGGGGGGTATGATACGGGG‑ACGAAAACCTTTTTCTTGTTAAGCAAAGAGGCCTCGGCAGCCACGATTTTTCCTCTGGTTCATTTGAAGGATTTAACCGGACAAAAAAATGGCTATTTTTATGCGACCGTTTTTTCGACAAAAAGCCCCATAGCAGccctt                                                   >  2:580997/1‑248 (MQ=255)
                                                                                                                                                                                                        atatGAAGAGGGGAAAGTATGTGGCGGTTATCAAC‑CCTCTCTTCTGTTATGAAGTCATGTCTGATGACTCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACTAAAACCTTTTTATCTTTAAGCAAAGAGGCCTCGGCAGAAACGGTTTCTTCGCTGATTCATTTGAAGGATTTATCCGTTCAACGAAATGGCTATTTAGATGCGACAGTTTATTCGAGAAATCGCTCTATTGCAgggagtggtgcgatt                                          <  2:308957/251‑15 (MQ=255)
                                                                                                                                                                                                               gCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATc                                  <  1:82424/251‑1 (MQ=255)
                                                                                                                                                                                                               gCGGGGAAATTAGGGGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGATGCCTCGGCAGCAACGTTTTCTCCGCTGAGTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGGTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATc                                  <  2:424520/251‑1 (MQ=255)
                                                                                                                                                                                                                cGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCa                                 <  1:168473/251‑1 (MQ=255)
                                                                                                                                                                                                                  gggAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTTCAATTGGTTGTGTAGGATACGGGG‑ACGAACACCTTTTTCTCGTTAAGCAAAGAGTCCTCGGCAACAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACCAAGAAATGGCTATTTATATGCGACATTTTATTCTACAAAACGCCCAATAGCAGCCAT‑GTGTGCAACTTCATATCAAc                               <  2:512496/251‑1 (MQ=255)
                                                                                                                                                                                                                          aTGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGGATGATACGGTG‑ACGAAAACCTTTTGCACGTTAAGCAAAGAGGCCGCGGCAGCAACGTTTTATCCGCTGATGAATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTtat                        <  1:199142/250‑1 (MQ=255)
                                                                                                                                                                                                                             gggCGGTTATCAAC‑CCTCGCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAAcc                    <  2:388380/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAACCCAttt               >  2:177746/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAACCCAttt               >  2:601775/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCGTATAACCCAttt               <  2:539068/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAACCCAttttt             <  1:500118/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   gctGGAGTGAAGTGATGTCTGATGACCCGACATTTCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACTCCCAATTGCAGCCATTGT‑TTCTACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCtt  <  2:515326/247‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AAAGCGAATGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAAC‑CCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTG‑ACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT‑TGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTT  >  NC_000913/4275710‑4276199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑