Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,6530CA30.6% 36.7 / 11.3 36L30F (TTG→TTT) pdeRcyclic di‑GMP phosphodiesterase PdeR
Reads supporting (aligned to +/- strand):  ref base C (23/2);  new base A (1/10);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 2.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGT  >  NC_000913/1346412‑1346897
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aaaGCCCGTGGAAAAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACTATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCg                                                                                                                                                                                                                                              <  1:345623/251‑1 (MQ=255)
       gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAATACATCGCTGCttt                                                                                                                                                                                                                                       >  2:87031/1‑251 (MQ=255)
        ggtaataaaccccgaaccaCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCGGCTTTc                                                                                                                                                                                                                                      <  2:373484/233‑1 (MQ=255)
                          aCCGCCGGTGTGTCATGCCATGCAGAAGCGCTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTAACAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATgg                                                                                                                                                                                                                    >  2:242090/1‑251 (MQ=255)
                          aCCGCCGGTCTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGCCTCATCATCCGGCTGGCGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCACAACCTCGCTGCTTTCCGTCAACCGCCAGTATgg                                                                                                                                                                                                                    >  2:250375/1‑251 (MQ=255)
                                                      cgTTGCCAGCCCATTCACGTTTCTT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGCTTGTGCGATCCGAGAAAGTTGTACacc                                                                                                                                                                                        >  2:257234/1‑249 (MQ=255)
                                                          tcccaccccatCCCGTTTTTTTGATTTTGCGCCCTACCAGATGCAAAAAAAGCGCGCTT‑CACTCGCCGTTCAGACTCCACAGCAGGCGGGAGGTGAAGACCCTCATTCCCC‑TG‑ATGAGAGCACCTTGT‑TCGGCAAATAACAGCAATGTACGATCAGTTTTTGGGTTTTGGGTGAAAAAAAAGC‑AAAAAACACCGCTATTTTCCAAAACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGt                                                                                                                                                                                      <  2:10244/240‑1 (MQ=255)
                                                          gCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCCCAGAAAGATCGCTTTCA‑TCGCCGATCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCAC‑TG‑ATGCGAGCAGCCTGT‑TCGGGAGATAACTGCACTGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGA‑GAAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtg                                                                                                                                                                                    >  1:339897/1‑251 (MQ=255)
                                                                  aTTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTcc                                                                                                                                                                            >  2:197478/1‑251 (MQ=255)
                                                                      ccGTTTATT‑GATTTTTCGCCCAACCAGATGCACAGAAAGACCGTTTTCA‑CCGACGGTCAAACTCATCATCAGGCGGGGGGTGAGGCCCTTCATCTCCC‑TG‑ATGCGAGCAGCCTGT‑TCCGCAAAAATCTGCAATGTACGATCAGGTTCTGTGTTTTGGGTGGAAAAAAAGC‑GAAAAACATCGCTGCTTTCCGCCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTaa                                                                                                                                                                        <  1:556880/251‑1 (MQ=255)
                                                                        gTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATAGCTTTCA‑TCGACGGTCAGACTCATCAAAAGGCTGGAGGCGATGACCGTCATTTCCA‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATACCTGCAATGTACGATCAGGTTCTGTGGTTTCGGAGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTAAACCGCCAGTATGGATTGTGAGATCAGAGAAAGATGTACAACGTTGTGGACTCCCTAACg                                                                                                                                                                      >  1:77271/1‑251 (MQ=255)
                                                                        gTTTACT‑GATTTTTCGTCCTAAAAGATGAACAGCAAGACCGCTTTCC‑TCGCCGGTCAGACTCATCATCCGGCTGGAGGTGATGACCGTCATTTCAC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTCCGATCAGGTTCTGTGGTTTAGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACg                                                                                                                                                                      >  2:272923/1‑251 (MQ=255)
                                                                              t‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTc                                                                                                                                                                >  1:440216/1‑251 (MQ=255)
                                                                                          ccTACCAGATGCACAGAAAGAACGCTTTCA‑TCGCCGGTCAGACTCATCATCCGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑CTGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCCTGTTGCatat                                                                                                                                                     >  2:55583/1‑251 (MQ=255)
                                                                                                aGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCg                                                                                                                                               >  1:174837/1‑251 (MQ=255)
                                                                                                              aTCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGt                                                                                                                                 >  2:200171/1‑251 (MQ=255)
                                                                                                                cGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTcc                                                                                                                               >  1:320072/1‑251 (MQ=255)
                                                                                                                  cTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgtg                                                                                                                                                                                  >  1:343642/1‑198 (MQ=255)
                                                                                                                   tttCA‑TCGACGGTCAGACTCATCATCAGGCTGGCGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCACAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcgcg                                                                                                                            >  2:91570/1‑251 (MQ=255)
                                                                                                                        accGACGGTCAGCCTACTCACCAGGTTGGGGGTGATGACCGTC‑TTTCCCGTG‑AGGCGAGCAGCCTGT‑TCGGAAGATAACTGCAATGTACGATCAGTTTTTGTGTTTTGGGTGGAAGAAAAGC‑GAAAACCATCGCTGCTTGCCGGCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTGTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta                                                                                                                        <  1:242090/249‑1 (MQ=255)
                                                                                                                            aCGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaat                                                                                                                     >  1:521006/1‑250 (MQ=255)
                                                                                                                                 cAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAgg                                                                                                               >  2:472881/1‑251 (MQ=255)
                                                                                                                                   caCTCATCATCAGGGTGGGGGGGA‑GGCCGTTCTTTCCC‑TGGAGGCGAGCAGCTTGC‑TGGGAAGATAACGACAATTTACGAGCGGTTTTTGTGGTTTCGGTGGAAAAAAAGC‑GAAAAACACCGCTTCTTCCCTTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaa                                                                                                             <  2:171995/250‑1 (MQ=255)
                                                                                                                                    ccTCATCATCGGGCTGGGGGTGAGCACCGTCTTCTCCC‑GG‑AGGCGAGCGGCCTGT‑TCGGAAAATAACTGCATTGCACGATTGTTTTCGGTGTTTGGGGTTG‑AAAAAAGCCGAAAAATATCGCGTCTTTCCGTCAGCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaac                                                                                                            <  1:250375/250‑1 (MQ=255)
                                                                                                                                    aCTCATCACCAGGCTGGAGGTTATGCCCGTCTTTCCCC‑TG‑AGGGGAGCA‑TCTGTGTCTGGAGATAACGGCAATGAACTATCGGGTTGTGATGTTTGGGGGGAAAAAAAGC‑AAAAAACTCCGCTGCTTCCCGTCATCCGCCAGTATGTATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGTCTCCCTAACGGTTTTCAGGATGTATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaac                                                                                                            <  1:272923/251‑1 (MQ=255)
                                                                                                                                        atcctcgggtggggggggatgcCCGTCTTTCCCC‑TG‑ATGCGAGCAGCCTGT‑CCGGAAGATACCTGCAATGTACGATCAGGTTCTGTGGTTTGGGTGGTAAAAAAGC‑CAAAAACCTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAacaaca                                                                                                        <  1:200171/234‑1 (MQ=255)
                                                                                                                                                                gTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcggcgg                                                                                >  1:3047/1‑251 (MQ=255)
                                                                                                                                                                gTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcggcgg                                                                                >  1:357743/1‑251 (MQ=255)
                                                                                                                                                                gTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcggcgg                                                                                >  1:114717/1‑251 (MQ=255)
                                                                                                                                                                 tctttcccc‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCATTGTCCGATCAGTTTCTGTGTTTTGGGTGGAAAAAAAGC‑GAAAAACACCGCTGTTTCCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcggcgg                                                                                <  1:197478/244‑1 (MQ=255)
                                                                                                                                                                                     cagcCTGT‑GCGGCAGATAACTGAAATGCACGATCAGGTTCGGTAGGTTCGATGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCAAGTATGGATGGTGCGATCCGAGAAAGTTGTACAACGTTGTGGAGTCCCTAACGGTTTTCGTGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATCAAAGGAAAAACATGACTTACATAAAATTAACGGCGGC‑TCAAAGCAGTGgtttttt                                                            >  1:438345/1‑251 (MQ=255)
                                                                                                                                                                                              ggggaagaaatcgaaaaTGAACGAGCAGTTTCTGTGTTTTGGGTGGAAAAAAAGC‑GAAAAACTCCGCTGTTTCCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TAAACGCATTGTTTTTTTTATAtatt                                                    <  1:312063/236‑1 (MQ=255)
                                                                                                                                                                                                                            tCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGATGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGACAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCGTAAA‑GCCTTGTTTTTTTTTTATATATATTTTTAATCCAGTTTTTGTAAAACGcc                       >  2:483429/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ggggggtAGAAAAGC‑AAAAAACATCGCTGCTTTCCGCCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCaaa             <  1:87031/246‑1 (MQ=255)
                                                                                                                                                                                                                                         tggtAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TTAACGCATTGTTTTTTTTTTTTATTTATTTTTAATCCAGTTTTGGAAAAACGCCCGTTTTTAAACGa          >  2:318526/1‑250 (MQ=255)
                                                                                                                                                                                                                                                aaaGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TAAACGCATTGTTTTTTTTTTATATTTATTTGTAATCAAGTTTTGGAAAAACGACCGTTTTCAAACGAAAGTCAGt  >  2:490699/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATT‑GATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCA‑TCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCC‑TG‑ATGCGAGCAGCCTGT‑TCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGC‑GCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGC‑TAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGT  >  NC_000913/1346412‑1346897

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