Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT36.1% 30.3 / 17.0 36V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (4/19);  new base T (13/0);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803639‑804119
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               
ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTTTTTTTGTGTCCCCGCCTA‑TTGctggc                                                                                                                                                                                                                                            >  1:609070/1‑247 (MQ=255)
      ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCTcccccttatgtggtcattcct                                                                                                                                                                                                                                       >  1:350820/1‑232 (MQ=255)
       cGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCCCGCTATTACCTTTGTAGGTATTGTACAGACCGCTACA‑GCGTCCGGTCTGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTCTTTTTTTGCCTCCTCATA‑TTTCgggccttcctc                                                                                                                                                                                                                                     >  1:96096/1‑242 (MQ=255)
              atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAACTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATATTTTTGGCTCCGCattttttggggatacccttatattt                                                                                                                                                                                                                               >  1:589495/1‑228 (MQ=255)
                                           aTTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCCTA‑TTTCGGGCATCCCTCAATATA‑TTGGTGGGGTGATGGCAGGGGGGGGCaaa                                                                                                                                                                                                 >  1:452453/1‑250 (MQ=255)
                                                    tcACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTCTGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCCCCTA‑TTGCGGTGATCCATCCATATA‑TGAGGGGGGCGAGGGCAGGCGGGGGCCAAATGATccccc                                                                                                                                                                                       >  1:563373/1‑248 (MQ=255)
                                                     cACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑Cgg                                                                                                                                                                                     <  1:560560/251‑1 (MQ=255)
                                                        gATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTc                                                                                                                                                                                  <  1:72215/251‑1 (MQ=255)
                                                                  ttgttgGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGcc                                                                                                                                                                          <  1:146892/249‑1 (MQ=255)
                                                                       gggaTTGTACAGACC‑CTAAAAGCGTCCGGTATGGAAGTTTTCCCGTTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCg                                                                                                                                                                   <  2:53893/248‑1 (MQ=255)
                                                                          aTTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCCTA‑TTTCGGTCATCCCTCAATATA‑TTCGTGGGGTGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGGCAGACAGCCACACCCG‑ttt                                                                                                                                                               >  1:592061/1‑250 (MQ=255)
                                                                             gTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTCCTCCCCATT‑TTTCGGGCATCCATCCATATA‑TGCGTGGGGTGATGGCAGGGGGCGGGCAAATGAACCT‑CGGGGCAAACAggcccaaccgcttcgg                                                                                                                                                             >  1:363605/1‑236 (MQ=255)
                                                                                cAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGACGTATGTGCTGACCAACTGCCCTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTTTTTTTTTGCCTCCCCATACTTG‑GGGCATCAATCAATATA‑TTCGTGGGAAGATGGCAGGGGGCGGGAAAATGATCCT‑CGGGGGAGACAACCACACCCC‑TTACGGTgg                                                                                                                                                         >  1:543309/1‑250 (MQ=255)
                                                                                           a‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTCCTCCCCATA‑TTGCGGTCATCCCTCAATATA‑TTCGGGAGATGGTTGGAGGGGGCGGGAAAATGATCCTACGGG‑CAGACAGCCCCACCCG‑GTACGGGGCATTAggggac                                                                                                                                               >  1:125975/1‑247 (MQ=255)
                                                                                                 ccGGTAGGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCa                                                                                                                                         <  2:355667/251‑1 (MQ=255)
                                                                                                    gTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTc                                                                                                                                      <  1:412482/251‑1 (MQ=255)
                                                                                                           aCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGCTTTTTTTGCCTCCGCCTA‑TTGCGGGCATCCCTCAATATA‑TTCGTGAGGTGATGGCCGGGGGCGGCAAAAAGAACCT‑CGGGTCCGACAGCC‑CACCCCCTTTCGGTGCATTTGGGGCAATGGccggcggggggg                                                                                                                               >  1:125621/1‑240 (MQ=255)
                                                                                                                 cccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGggcggc                                                                                                                          <  1:319251/250‑1 (MQ=255)
                                                                                                                    gcggcgggaTGTGCGGACCAACTGCCATAACCCCCTCTGCGCCGTCGGCGGCACTAGTAACGGTGATGACCATGTTTGTGGGTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTg                                                                                                                      <  2:262162/243‑1 (MQ=255)
                                                                                                                         cGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATAGA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACGCCCG‑TTACGGTGCATTAGGGACAATGGCAGGCGGTGAGGGCGGCGGTGAGTT‑g                                                                                                                <  1:429938/251‑1 (MQ=255)
                                                                                                                                          cgCCATAACCCACTCGGCGCCGCCGGCGGCCCTATTAACGGGGAGGACCATGTTTTGGGTTTATGGGCGGCCCAGGGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGGGAGATGATGGCAGGCGGCGGGAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACATTGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTaa                                                                                               <  2:239327/250‑1 (MQ=255)
                                                                                                                                                 actcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACAcc                                                                                        <  1:160038/251‑1 (MQ=255)
                                                                                                                                                       cctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGat                                                                                  <  2:250104/250‑1 (MQ=255)
                                                                                                                                                          gcgcCGGCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGatat                                                                                <  2:383259/250‑1 (MQ=255)
                                                                                                                                                                gcggcggCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCAGCAATATA‑TGCGGGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAGCAGCTGCTTAAGGACACCTGGGATATCGACTAt                                                                         <  2:303326/250‑1 (MQ=255)
                                                                                                                                                                  ggcggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTA‑TTTAGGGCATCCATCAATATAATG‑GGGAGATGATGGGAGGCGGGGGGAAAATTATCCC‑CGGGTCAGAAAACCCCACCCG‑TTACGGGGCATTAGGGACAAAGGGAGTCGGTGAGGGGGGGGGTGAGTTTGGAAAAACA‑CAGCTTTATGACACCTGGGATATCGTCTATcc                                                                       >  1:149967/1‑251 (MQ=255)
                                                                                                                                                                  ggcggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATcc                                                                       <  1:55080/251‑1 (MQ=255)
                                                                                                                                                                      gcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGc                                                                   <  2:235854/251‑1 (MQ=255)
                                                                                                                                                                             tAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCCCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGtt                                                             <  1:150508/250‑1 (MQ=255)
                                                                                                                                                                                     atgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCAt                                                    <  1:444361/251‑1 (MQ=255)
                                                                                                                                                                                               ttttGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCATA‑TTTCGGGCATCCATCCATATA‑TGCGTGAGATGGTGGCAGGCGGGGGCCAAATGATCCT‑CGGGGCAGACAACCACACCCC‑TTACGGGGCATTTGGGACAATGGCAGTCGGTGAGGGCGGCGGGGAGGT‑GGTAAAACAGCTGCTTAATGACAACCGGGATATCGACTATCCGGGCgggggttggggggatcttaccggg                                           >  1:333773/1‑227 (MQ=255)
                                                                                                                                                                                                   ggTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCCCA‑TTTCGGGCATCCATCCATATA‑TGCGTGAGATGATGGCAGGGGGGGGCAAAATGAACCC‑CGGGGCAGACAACCACACCCC‑GTACGGGGGATTAGGGGCAATGGCAGGCGGGGGGGGGGGGGGGGGGGT‑TGTGAAAAAGGTGCTTAATGACACCTGGGGTATCGACTATTCGGGCGTGGGTGGGGGGGATCTGGACGGAAAAcc                                      >  1:458062/1‑251 (MQ=255)
                                                                                                                                                                                                    gTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTA‑TTGCGGGCATCCATCAATATA‑TGCGTGGGATGATGGCAGGCGGGGGCAAAATGATCCT‑CGGGGCAGACAACCACACCCG‑TTACGGGGCATTAGGGACAAAGGCAGTCGGGGGGGGGGGGGGGGAGTT‑GGGAAAAAAGCTGGTTAATGACACCCTGGATATCGACTATCCCGGCGTGGGTGCGGGGGATCTGACCGGAAAAccc                                     >  1:124650/1‑250 (MQ=255)
                                                                                                                                                                                                                cccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGtt                                                                                                                  <  2:503325/163‑1 (MQ=255)
                                                                                                                                                                                                                cccAGCGTTATGGCGGAATTTTTTTGCCTTCGCATA‑TTGCGGGCATCCCTCCATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATTCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTTGGGACAATGGCAGTCGGGGGGGGcggcgggggggt                                                                                                                  >  1:503325/1‑157 (MQ=255)
                                                                                                                                                                                                                           ggcggcATTTTTGTGCCTCAGCATA‑TTGCGGTCATGCTTCAATGTG‑TGCGTGAGATGCTGGCAGGCGGCGGCCAAATGATCAT‑CGGGTCAGACAGCCACACCCG‑CTACGGTGCATTAGGGACAATGGCAGTCGGTGCGGGCGGCGGTGAGTT‑GGTAAAACCGCTGCTTAATGACACCTGGGGTATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCCCAgg               >  2:41472/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gtgCCTCCGCATA‑TTGCGGTCATACATAAATATA‑TAAGCGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAAAGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTggagatggc  >  2:73396/1‑244 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               
ATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAA‑GCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATA‑TTGCGGTCATCCATCAATATA‑TGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCT‑CGGGTCAGACAGCCACACCCG‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTT‑GGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803639‑804119

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