Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG27.7% 68.6 / 10.3 47V288G (GTG→GGG) cusAcopper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base T (1/33);  new base G (12/1);  total (13/34)
Fisher's exact test for biased strand distribution p-value = 3.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCG  >  NC_000913/599327‑599826
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                              
aaGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                                                                                                                               <  1:671981/251‑1 (MQ=255)
      aaaGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGGTCGTCGATCGAACGGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                                                                                                                                                                                                                                                         <  2:45045/251‑1 (MQ=255)
         gcgcgcTGGATGCGTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAAGATATGGTGCGCGCCAGCGGCTATCTGCAAACGCGCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTAGCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTg                                                                                                                                                                                                                                                      <  2:71308/251‑1 (MQ=255)
          cgcgcTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                                                                                                                                     <  1:150239/251‑1 (MQ=255)
                 gATGCGTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATGGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg                                                                                                                                                                                                                                             <  2:337665/251‑1 (MQ=255)
                                       gcggTTCGTCGATCGAACTGGCGGAAGCGGAAAATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                                                                                                                      <  2:641438/251‑1 (MQ=255)
                                          gTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                                                                                                                                   <  1:439946/251‑1 (MQ=255)
                                             cgtcgATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCTATCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                                                                                                                                <  2:92673/251‑1 (MQ=255)
                                                              gaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                                                                                                                             <  1:367201/251‑1 (MQ=255)
                                                               aagcggaaTATATTGTGCGCGCCGGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCAGTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                                                                                                                            <  2:350342/251‑1 (MQ=255)
                                                               aagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                                                                                                                            <  1:45924/251‑1 (MQ=255)
                                                                ggcggaaTATAGGGTGCGCGCAAGCGGCTATCTACAAACCCTCGACGACTTTAATCACATCGTTTGAAAAGCCAGTGAAAAGGGCGTGCCCGTTTATCTGCGCGATGTGGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                                                                                                                           <  2:612834/250‑1 (MQ=255)
                                                                        atatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                                                                                                                                    <  1:132244/250‑1 (MQ=255)
                                                                            ggTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTGGCGAATGTCCAGATTGGCCCGGAGATGCGCCGGGGCAGTGCCGAACTAAACGGCGAGGGCGAGGTGACGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCGCGATAAGCTATCTCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                                                                                                                                <  1:398174/250‑1 (MQ=255)
                                                                            ggTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGACCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTGCTGCGATCCGGCAAAAACGCCCCAGGAGTGACTGGCGCCGGGGAGGACAAACCGGGAACGCTGTAAagg                                                                                                                                                                                >  1:67659/1‑249 (MQ=255)
                                                                               gggcgccagggtcaaacTGCAAACGCTCGCCTACTTTAATCGCATCGTTTTAAAAGCCAGTAAAAAGGGCGTGCCCGTTTATCTGCGTGATGTTGCGAAGGTCCAGATTGGCCCGTAGATGCGCAGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGTGGGGCGGGGTGGTGCTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                                                                                                                            <  2:184326/235‑1 (MQ=255)
                                                                                                   aaaCGCTCGACGCCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                                                                                                                        <  1:521034/251‑1 (MQ=255)
                                                                                                      cGCTCGACGACTTTAATCACATCGTTTTACAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAGCCGGCCACAACGCCCGGGAAGGGATCGGCGCCCTGGAGGGGAAACTGTGAAAGCTGGAAAATAGTCTggcgcggaggggggggggaaat                                                                                                                                                     >  2:667081/1‑230 (MQ=255)
                                                                                                           gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGGCAAAAAGCCCGGGAAGGGATAGGCGCCGGGGAGGGCAAAATGGGAACGCTGAAAAGTAGGCTGCCGGGCGGCGGGGGGATAGTTtcaac                                                                                                                                                >  1:71665/1‑251 (MQ=255)
                                                                                                             cgccttgactcacttcgctgtAAAAGCCAGTGAAAATGGCGTGCCCGTTTACCGGCGCCATGGTGCGAAGGTCCAGATGGGCGCGGAGAGGCGCCGGGGTATTGCGGAAGTAAACGGGGAAGGCGAAGTCGCGGCCGTGTTGCTGATCCTACGCTCCTGCAGAAAGGCCCTAGAAGTGATCGCCGCCGTGAAGGACACGCTGGAAACGCTGAAAAGTAGTCTGCCGGAATGCGTGTAGGTAGTTCCAACAt                                                                                                                                              <  2:672507/231‑1 (MQ=255)
                                                                                                                    aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                                                                                                                        <  1:590646/250‑1 (MQ=255)
                                                                                                                     aTCACAGCGTTTTAAAAGCCAGGGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATGGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAATTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGGGTGGAGATAGTTACAGCATACGATCGc                                                                                                                                      <  2:660695/251‑1 (MQ=255)
                                                                                                                        acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCGAGAAGTGGTCGCCGCCGCGAAGGGAAAACAGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGGTACAACATACGAACGccac                                                                                                                                   >  1:604277/1‑248 (MQ=255)
                                                                                                                            cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                                                                                                                                <  1:632621/250‑1 (MQ=255)
                                                                                                                                ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                                                                                                                            <  1:489507/250‑1 (MQ=255)
                                                                                                                                   aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGGACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATTCCGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGGAAAGGAATCTGGCGGGAGGCGGGGGGAGAGTTACAACATACGATCGcagacagccccttg                                                                                                                         >  1:333343/1‑244 (MQ=255)
                                                                                                                                     aGCCAGTGAAAGTGGCGTGCCCGTTTATCTGCGCGATGTGGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATGGCCGAACGAAACGGCGAAGGCGAAGTGGCGTGCGGGGTGGGGATCCTGCGATCCGTCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAGCGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                                                                                                                      <  2:329834/251‑1 (MQ=255)
                                                                                                                                      gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCGGGCAGGACAAACCGGAAACGCGGAAAAGTAGTCTGCCGGAAGGCGTGGAGGTAGTTTCAAAAAAAGAACGAAGCCAGCTccattaccc                                                                                                                     >  1:187542/1‑243 (MQ=255)
                                                                                                                                       ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc                                                                                                                    <  1:323396/251‑1 (MQ=255)
                                                                                                                                          gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc                                                                                                                 <  1:106718/251‑1 (MQ=255)
                                                                                                                                              aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGCGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGGTATGCGATCCGGCCAAAACGCCCGAGAAGTGAACGCCGCCGTGAAGGACAAACCGGGAACGCTGGAAAGTATTCTGCCGGCAGGGGGGGAGATAGTTACAAAATACGATCCCCGCCCGCTCAGTGACCGCGCCATCg                                                                                                             >  1:329341/1‑251 (MQ=255)
                                                                                                                                                  ggCGTGCCCGTTTATCTGCGCTATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGTCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACTCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa                                                                                                         <  1:643297/251‑1 (MQ=255)
                                                                                                                                                                  ggggcgagggtGCGAAGGTCGCGATTGGCCCGGAGAGGCGCCGGGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACTCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCGGAAAAGTAGTCGGCCGGAAGGCGTGGAGATAGTTACAACATACTATCGCAGCCAGGTCATTGACCGCGGCATCGACAACCGCAGCGGCAAGTTg                                                                                         <  2:315800/241‑1 (MQ=255)
                                                                                                                                                                   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAg                                                                                            <  1:459876/247‑1 (MQ=255)
                                                                                                                                                                             gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGt                                                                               <  1:140169/250‑1 (MQ=255)
                                                                                                                                                                                    tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGCTCCGGGCAAAAAACCCCCGAAGGGGTCGCCGCCGCGGAGGACAAACTGGGAAAGCCGGAAAGGAGTCTGCCGGGAGGCGGGGGGATAGTTAAAACATACGATCGCAGCCAGCTCATTGTCCGCGCCATCGCCAACCCCAGCGGCCAGCTGCTGGTAGAGGTTTTtgtg                                                                       >  2:500856/1‑251 (MQ=255)
                                                                                                                                                                                     ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGTCCGGCCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGGCTGCCGGGAGGCGTGGGGATAGGTACAACATACGGTCGCTGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGGTTATTgtgg                                                                      >  1:496450/1‑251 (MQ=255)
                                                                                                                                                                                                cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGTTTAGTGTGGTGGCGGTAGTc                                                           <  2:67659/251‑1 (MQ=255)
                                                                                                                                                                                                     gtGCGCCGGGGGTGTGCCGAATGAAACGGAGAAGGGGAAGGGGCGGGAGGGGGGGTGAGCCGGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGTTTAGTGTGGTGGCGGTAGTCTGTgc                                                      <  2:218202/250‑1 (MQ=255)
                                                                                                                                                                                                           cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGCTCCTGCGCTCCGGCAAAAACGCCCGAGAAGAGATCGCCGGCGCGAAGGGCAAACTGGGAAACCCGAAAAGTAGTCTGGCGGGAGGGGTGGGGATAGTTACAACATACGGTCCCAACCAGCTCATTGAACGCGCCACCCACAACCTCCGCGGCAAGCTGCTGGTAGAGATTATTGTGGTGGCGGGAGGCtgtggggtgtt                                                >  1:157783/1‑244 (MQ=255)
                                                                                                                                                                                                            ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCCGGCCAAAACGCCCGAGCAGTGATTGCCGCCGTGAAGGACCAACTTGGAACGCTGTAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCCTTTACCCCGCCCTCGACAACCTCAGCGGCAAGTTTCTGGGAGAGTTTATTGTGGGGGCGGGAGGCTGTGCGCTGttt                                               >  1:486049/1‑251 (MQ=255)
                                                                                                                                                                                                                     ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAt                                      <  1:436379/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGAGATCCGGCAAAAACGCCGAGGAGGGATCGCCGCCGTGAAGGACAAACTGGGAAAGCTGAAAAGTAGGCTGCCGGGAGGGGGGGGGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGGAGAGTTTTTTGTGGGGGCGGGAGGCTGTGCGCTGGTTTTCTgggctgtggggtc                              >  1:649108/1‑240 (MQ=255)
                                                                                                                                                                                                                                         gaagTGTCGGGCGGGGTGGTGATCGTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCGGTGCGCTGTTTCTCTGGCATGTGCGCTCGTCGCTGGTGGc                  <  1:437496/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gggggggTGGTGATCCGGCGGTCCGGGAAAAACGCCCGAGAAGTGATCGACGCCGAGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTc         <  2:379937/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGt       <  1:136210/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGCAGTGATCACCGCCGTGAAGGACACACTGGAAACGCTGACACGTCGTCTGCCGGAAGCCATGGAGATAGTTACAACATACGATCCCAGCTAGCCCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt      >  1:316803/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     ggggTGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt      <  2:491730/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        gtggtgATCCTGAGATCCGGCAAAAACGCCCGAGAAGTGATAGCCGCCGTGAAGGACAAACTGGAAACGCTGAACAGTAGTCTGCCGGCAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAACTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGCCTGTGCGCTGTTTCTCTGGCCTGTGCGCTCGGAGCTGGTGGCGATTATTTCGTTGcc   >  1:324821/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         tggtgCTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGCTCGCCGCCGTGAAGGACAAACTGGAAACGCTGAACAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTTGCGATTATTTCGTTGccc  >  2:536874/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                              
AAGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCG  >  NC_000913/599327‑599826

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: