Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG39.5% 18.5 / 20.5 39V288G (GTG→GGG) cusAcopper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (14/1);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.72e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTC  >  NC_000913/599335‑599819
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aGCGCGCTGGATGCGTCAAACCATGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCTCCAGCGGCTATCTGCAAACGCTCGTCGACTTTAATCACATCGTGTTAAAAGGCAGTGAAAATGGCGTGCCCGTTTAGCTGCGCGATGTTGCGAAGGGCCAGATTGGCCCGGAGATGCGCCGGGTCATTGCCGAACTAAACGGCGAAGTCGAAGTGGCGGGCGGGGTGTTGATCCt                                                                                                                                                                                                                                             <  1:80897/251‑1 (MQ=255)
       gggATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGAGCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAAGCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc                                                                                                                                                                                                                                      <  2:123510/250‑1 (MQ=255)
             cTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTCGGATCCGGCac                                                                                                                                                                                                                                 >  1:329131/1‑249 (MQ=255)
                    ccAGTAAGCGGGCGGGTCGTCGATGGAACGGGCGGAAGAGGAATATATGGTGCGCGTCAGCGTCTAACTGCAAGCGCTCGACGACTGTAATCACATCGGTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCTATGTTGCTAAGGGCCAGATTGGCCCGGTTATGGGCCGGGGTATTGCCGAATTAAACGGCGAAGGCGAAGTGGAGGGCTGGGTGGGGTTCCTGCGATCCGGCAAAAACGccc                                                                                                                                                                                                                         <  2:318932/251‑1 (MQ=255)
                    ccAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGGACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATcacgcgcgccggcgaaaaatccc                                                                                                                                                                                                                         >  2:382195/1‑229 (MQ=255)
                             gggcggTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCCGTCCCGCAAAAAAGCCCggggagggg                                                                                                                                                                                                                >  1:436039/1‑243 (MQ=255)
                                        ggatcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                                                                                                                     <  2:220806/250‑1 (MQ=255)
                                                     ggaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg‑g                                                                                                                                                                                       <  1:61238/251‑1 (MQ=255)
                                                           ggaaTATATGGTGCGCGCCAGCGGTTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAATTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGTTGGTGATCCTGCGATCCTGCAAATACGCCCGAGGATGGATCGCCGCCGTGAATGACAAACTG‑GAAACGc                                                                                                                                                                                 <  2:268783/251‑1 (MQ=255)
                                                            gaaaatatGGTGCGCGCCAGGGGGTATCTGCAAACGCGCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTAGCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCt                                                                                                                                                                                <  2:490007/247‑1 (MQ=255)
                                                                atatgtgtcgcgaCAGCGGCTTTCTGCAAACGCTCGACGTCTTTAATCACATCGTTTTAAAAGACAGTGAAAATGGCGTGCCCGTTGATCTGCGCGAGGTTGCGAAGGTCCAGATTGGCCCGGGGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCTGGGGGGTGAGCCTGCGACCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTTAAGGACAAACTG‑GAAACGCTGaaa                                                                                                                                                                            <  2:530736/243‑1 (MQ=255)
                                                                atatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGaa                                                                                                                                                                             <  1:376609/250‑1 (MQ=255)
                                                                atatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTGCGGGCAAAACGCCCGCGAAGCGATCGCCCCCGTGGAgggggaaactgaaacgcatgc                                                                                                                                                                              >  2:94991/1‑231 (MQ=255)
                                                                    ggTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAagt                                                                                                                                                                         <  1:413499/250‑1 (MQ=255)
                                                                       gcgcgcCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAagtaggc                                                                                                                                                                     <  2:392890/251‑3 (MQ=255)
                                                                           gcCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGGAGGGCAAACAG‑GAAAAGCTGGAAAGTAGGCTGcc                                                                                                                                                                 >  1:375412/1‑251 (MQ=255)
                                                                            ccAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGCGAAGGGATGGCCGCCCGGAGGGACAAACTG‑GGAACGCTGGAAAGTAGTCTGCCg                                                                                                                                                                >  1:28227/1‑251 (MQ=255)
                                                                                           aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑Gaaa                                                                                                                                                                                    <  1:429562/216‑1 (MQ=255)
                                                                                           aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGGAGGACAAACTG‑GAAACGCTGGAAAGTAGTCTGCCGGGAGGCGGGgag                                                                                                                                                    >  1:423347/1‑248 (MQ=255)
                                                                                               gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGGGGTACGGCAAAAACGGCCGGGGAGGGGTCGCCGCCGTGTAGGGAAAACTGTGAAA‑GCTGGAAAGGAGTTTGCCGGAAggggggggggttattt                                                                                                                                             >  1:190752/1‑237 (MQ=255)
                                                                                                    acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGAAGTGATCGCCGCCGCGGAGGACAAAATG‑GGAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTacaaa                                                                                                                                         >  1:162483/1‑249 (MQ=255)
                                                                                                                                                        tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGGTACTGCGAGCCGGCGAAAACGCCAGCGGAGGGATCGCCCCCGGGAAAGGCGAACCG‑GAAACACCGGAAAGGAGTCTGCCGGAAGGCGGGGGGATAGTTAAAAAAAAAGATGGCACACCGGTCCTTGACCCCGCCATCGACAACCTCaaccggacat                                                                                    >  1:374004/1‑242 (MQ=255)
                                                                                                                                                            cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGAGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GGAACGCCGAAAAGAAATCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATCGCAGCCAGCCCATTGACCGCGCCATCGACAACCTCCGCGGCCAGTTGCt                                                                                >  1:67300/1‑251 (MQ=255)
                                                                                                                                                             gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTg                                                                               <  1:90695/251‑1 (MQ=255)
                                                                                                                                                                  gTTGCGAAGGTCCAGATTGGCCCGGAGATGGGCGGGGGCATTGCCGAGCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGGTCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAaga                                                                          <  2:329131/251‑1 (MQ=255)
                                                                                                                                                                    tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTTCAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCCCAGCGGCAAGGTGCTGGGagagg                                                                        >  1:182789/1‑250 (MQ=255)
                                                                                                                                                                             ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGGGATCCGGCAAAAACGCCCGAGAAGAGATCGCCGCCGGGAAGGGCAAACTG‑GGAACGCTGATAAGGAGTCTGCCGGGAGGGGTGGAGGTAGTTACAACAAACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCCGCGGCACGATGCTGGAGAAGTTTAtttggg                                                               >  1:410872/1‑247 (MQ=255)
                                                                                                                                                                               aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTgtggtg                                                             <  1:416536/251‑1 (MQ=255)
                                                                                                                                                                                   tGGCCGGGAGAGGCGCCGGGGGATTGCCGAACTAAACGGCGAAGGGGAAGTGGTGGGCGGGGTGGTGATCCTGCGATCCGTCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATTCGATCGCAGCCAGCTCATTGACCTCGTCATCGACAACCGCAGCGGCAAGTTGCTTGAAGAGTGTATTGTGGTGGCg                                                          <  2:34976/250‑1 (MQ=255)
                                                                                                                                                                                                    ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGttt                                        <  1:54904/251‑1 (MQ=255)
                                                                                                                                                                                                      ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGGCGTGAAGGACAAACTG‑GAAAAGCTGAAAAGTAGTCTGCCCGGAGGGGGGGGGATAGGTACAACACACGATCGCCGCCAGCTCATTGACCGCGCCATCGACAACCTCCGCGGCAAGGTGCTGGAAGAGGTTTTTGTGGGGGCGGGAGTCTGTGCGCTGttttt                                      >  1:504792/1‑249 (MQ=255)
                                                                                                                                                                                                                     aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTc                       <  1:515453/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGGGAAGGGCAAACTG‑GAAACGCTGGAAAGGAGTCCGCCGGGAGGGGTGGAGAAAGTTACAACATACGATCGCAGCCAGCTCATTTACCGCGCCATCGACAACCCCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGGAGTCTGTGCGCTGTTTTTCTGGCATGTGCGCttgggg                   >  1:102575/1‑246 (MQ=255)
                                                                                                                                                                                                                         gcgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACCCCCGGGAAGTGATCGCCCCCGTGAAGGACCAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGAGGTAGTTACAACATACGATCGCAGCCCGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGGTTATTGTGGTGGCGGGAGGCTGGGGGCTGTTTttcttggctttgggcccggcg                   >  1:206634/1‑234 (MQ=255)
                                                                                                                                                                                                                           gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCt                 <  2:102575/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGat        <  2:423347/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggcgggcgGGGTGGTTATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGatt       <  1:157225/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGattatt    >  1:438145/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTc  <  2:182789/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
AGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑GAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTC  >  NC_000913/599335‑599819

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