Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG40.0% 99.4 / 12.5 60V288G (GTG→GGG) cusAcopper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base T (6/30);  new base G (24/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 3.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTT  >  NC_000913/599353‑599822
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
aaCCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTGCGG‑CAAAAACGc                                                                                                                                                                                                                                   >  1:518966/1‑251 (MQ=255)
  ccAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGccc                                                                                                                                                                                                                                 <  1:515982/251‑1 (MQ=255)
     ggAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCgag                                                                                                                                                                                                                              <  1:335804/251‑1 (MQ=255)
         gcgggctgTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGt                                                                                                                                                                                                                          <  2:130576/251‑1 (MQ=255)
               ggTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGTGGGGTGGTGATCCTGCGAt                                                                                                                                                                                                                                                 <  1:338830/223‑1 (MQ=255)
                gTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATcgcc                                                                                                                                                                                                                   <  2:27098/251‑1 (MQ=255)
                        atcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTACCGCGATCCGG‑CAAAAAAGccccggaagggaaagccgccctga                                                                                                                                                                                                            >  1:340993/1‑229 (MQ=255)
                        atcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGTCCGG‑CAAAAAAGCCCCGGGAGGGATCGCCGCCGTgg                                                                                                                                                                                                            >  1:169152/1‑249 (MQ=255)
                                   ggaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑Tgg                                                                                                                                                                                              <  2:222372/251‑1 (MQ=255)
                                          gaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑t                                                                                                                                                                                      <  2:452317/251‑1 (MQ=255)
                                            atatatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGa                                                                                                                                                                                    <  2:17201/251‑1 (MQ=255)
                                                atGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGTCCCG‑CAAAAACGCCCCAGGAGGGGTCGCCGCCCTG‑CAGGGCAAAC‑CGGGAACGcaggaaa                                                                                                                                                                                  >  1:100430/1‑244 (MQ=255)
                                                      cgcgcCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGGGCCTTCGGTCCGG‑CCAAAACACCCCAGAAGAGATCGCCGCCGTG‑GAGGGCAAACGTGAAAAC‑C‑CGAAAAGTAGTCt                                                                                                                                                                          >  2:280783/1‑251 (MQ=255)
                                                            aGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGACCGG‑GAAAAACGCCCCAGAAATGATGGCCGCCGTG‑AAGGACAAAA‑TGGAAAccgtgcaacattacggtaccgc                                                                                                                                                                     >  1:335091/1‑230 (MQ=255)
                                                             gCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGG‑CAAAAAAGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGGAAAGC‑TGGAAAGGAGTCTGGCggg                                                                                                                                                                    >  1:234400/1‑249 (MQ=255)
                                                                gatATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCTGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGc                                                                                                                                                                <  2:268109/249‑1 (MQ=255)
                                                                  tATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGGTCCCG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGTG‑AAGGACAAAA‑TGGAAAAGC‑TGGAAAGGAGTCTGCCGGGAggggg                                                                                                                                                              >  1:32560/1‑248 (MQ=255)
                                                                       gCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAAAGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCGCTTCGG‑GACAAACGCCCGGGAAAGGAACGCCGCCGGG‑CAGGGAAAAC‑TGGAAACGA‑CGGAAAGGAGTCTGCCGGGAGGGGTggtgg                                                                                                                                                         >  1:230061/1‑248 (MQ=255)
                                                                               tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTGCGG‑GAAAAACCCCCCCGAAGGGAACGGCGGCGTGGAGGGACAA‑C‑TGGAAACGC‑TGACAAGCAACCTGCCGGGAGGCGGGGGGATAGATa‑cc                                                                                                                                                >  1:116788/1‑250 (MQ=255)
                                                                                   cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CAt                                                                                                                                           <  2:343328/251‑1 (MQ=255)
                                                                                        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGAGGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACg                                                                                                                                        <  1:139288/249‑1 (MQ=255)
                                                                                          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATc                                                                                                                                     <  2:116788/250‑1 (MQ=255)
                                                                                               caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGGGGGGGGGGGGGGGGACCTGCGGACCGG‑GCAAAACGCCCGAGAAGTGGTCGCCGCCCGG‑CAGGACAAAC‑TGGAAAAGC‑TGAAAAAGAGTATGCCGGGAGGCGGGGAGAGAGTTA‑AAA‑CATACGATCGcagcc                                                                                                                               >  1:96583/1‑251 (MQ=255)
                                                                                                         aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGTCCCGCGAGCCGG‑CCAAAACGCCCGAGAAGCGGCAGCCCCCCGG‑AAGGGCAAAA‑TGGGAAAGC‑TGAAAAAAAGTCTGCCCGGAGGGGTGGGGATAGTTA‑AAAACAAAC‑ATCGCAACCAGCTCCTTgt                                                                                                                     >  2:86156/1‑250 (MQ=255)
                                                                                                            gCCAGTGAAATGGGCGTGACCTTGTATCGGCGCGATGTTGCAAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACAAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCAATCCGG‑CAAAAACGCACGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTAATTGAcct                                                                                                                  <  2:76315/251‑2 (MQ=255)
                                                                                                                gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTGCCGCGATCCCG‑GAAAAACGCCCGGGAAGGGATCGCCGCCGGG‑GAGGAAAAAC‑TGGAAACGC‑TGAAAAGGAGTTTGGCGGAAGGCGGGGAGGTAGTTA‑CCA‑CATACGAACGCAGCCAGCTCCTTGAACGCGcc                                                                                                              >  1:421171/1‑251 (MQ=255)
                                                                                                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCAt                                                                                                            <  1:292651/251‑1 (MQ=255)
                                                                                                                   aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATc                                                                                                           <  1:403467/251‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCGGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TAGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATc                                                                                                           <  2:353586/248‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc                                                                                                        <  1:203323/251‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGG‑CAAAAAAGCCCGCGAAGGGATCGCCGCCGTG‑AAGGGCAAAC‑TAGAAAAGC‑TGAAAAGTAGTCTGCCGGGAGGGGTGGGGGTAGTTA‑CAA‑CATACGAACGCAGCCCGCTcctttacccccccctcaa                                                                                                         >  1:353586/1‑233 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑GAAAAACGCCCCAGAAGGGATCGCCGCCGGG‑AAGGGCAAAC‑TGGAAACGC‑TGGAAAGAAGTCTGCCGGGAGGCGGGGGGATAGGTA‑AAA‑CATACGAACGCAGCCCGCTCCTTGACCGCGCCCTCggc                                                                                                        >  1:228449/1‑249 (MQ=255)
                                                                                                                          cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGCG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaaac                                                                                                    >  1:361429/1‑249 (MQ=255)
                                                                                                                                 gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATCCGG‑CAAAAACGCCCGAGAAGTGGTCGCCGCCGGG‑ACAGACAAAA‑TGGGAAAGA‑TGAAAAGAAGTCTGCCGGCAGGGGGGGGGAGATATA‑CCA‑AATACGATAGCAGCCAGCTCATTGACCGCGCCATCG‑CAACCCTCAGCgg                                                                                            >  2:44610/1‑251 (MQ=255)
                                                                                                                                     aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTCCGATCCGG‑CCAAAACGCCCGAGAAGGGATGGCCGCCGCG‑AAGGACAAAA‑TGGTAACACAGGCAAAG‑AGGCTGGCGGAAGGCGGGGGGAGAGTTAACAC‑CA‑ACGATCCCAGCCAGCTCCTTGAACGCGCCCTCCGCAACC‑TCCGCGGCa                                                                                          >  2:66588/1‑249 (MQ=255)
                                                                                                                                     aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGG‑CAAAAACGCCCGGGAAGGGATCGCCGCCGTG‑GAGGGCAAAC‑TGGAAACGC‑TGAAAAGGAGTCTGCCGGAAAGGGGGGGGATAGGTA‑CAA‑CATACGAACGCAGCCAGCTCCTTTACCGCGCCCTCGAGACCC‑CCAGCgggc                                                                                          >  2:477623/1‑247 (MQ=255)
                                                                                                                                        tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCGGACCGG‑CAAAAAAGGCCGCGAAGTGATCGCCGCCGCG‑AAGGACAAAC‑TGGAAAACC‑CGAAAAGAAGGATGGCGGAAGGGGGTGAGATAGTTA‑CAA‑AATACGATCCCCGCCCGCCCATTTGCCGCGCCATCCACAACC‑TCAGCGGcccgttg                                                                                     >  2:403016/1‑245 (MQ=255)
                                                                                                                                          cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGGCAAAAA‑GCCCGAGAAGGGATCGCCGCCGGG‑GAGGGCAAAA‑TGGAAACGC‑CCAAAAGGAGTCTGGCGGGAGGCGGGGGGGTAGTTA‑CAA‑AAAACGACAGCAGCCCGCCCATTGAACGCGCCacccaacacctcccccgggagtggcg                                                                                    >  2:295264/1‑226 (MQ=255)
                                                                                                                                                            cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTgtggt                                                                 <  1:120410/251‑1 (MQ=255)
                                                                                                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGc                                                              <  2:169152/251‑1 (MQ=255)
                                                                                                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCG‑CAAAAACGCCCGCGAAGGGATCGCCGCCGTG‑AAGGGCAAAA‑CGGAAACGC‑TGAAAAGTAGTCTGCCGGGAGGGGTGGGGGTAGTTA‑AAA‑CATACGACGGCCGCCAGCCCATTTACCGCGCCCTCGACCACC‑CCAACGGCAAAGTGGTGGAAAAGGTTATTGtggtggg                                                              >  2:219911/1‑250 (MQ=255)
                                                                                                                                                                    cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAg                                                         <  1:280783/251‑1 (MQ=255)
                                                                                                                                                                       ggAGATGCGCCGGGGCCTTGCCGACCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAACCGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAATGCGTGGAGATAGTTA‑CAA‑CATCCGATCGCAGCCAGCTCATTGACCGCGCCATCGACGACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCt                                                      <  1:261918/251‑1 (MQ=255)
                                                                                                                                                                           aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCCGG‑CAAAAAAGCCCGAGAAGGGAGCGCCGCCGGG‑AAGGACAAAC‑CGGAAAAGC‑GGAAAAGGAGGCTGCCGGAAGGCGGGGGGATAGGTA‑CAA‑CAAACGGTCGCCGCCAGCTCATTTACCGCGCCCCCGACAAAC‑TCCGCGGCCAGGTGCTGGGAGGGTTTTTTGTGGGGGCgggaggctgggg                                                  >  2:51484/1‑241 (MQ=255)
                                                                                                                                                                               gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGAACCCG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGGG‑AAGGGCAAAC‑TGGAAAAGC‑TGAAAAGTAGTCTGCCGGGAGGCGGGGGGAAAGTTA‑CAA‑AATACGAACGCCGCCCGCTCATTTACCGCGCCCTCGACAACC‑CCAGCGGCAAGGTGTTGGGAGGGTTTTTTTTGGGGGCgggggggtcgggggcg                                              >  2:335867/1‑237 (MQ=255)
                                                                                                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGAGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTgg                                     <  2:41067/250‑1 (MQ=255)
                                                                                                                                                                                          gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCCGTTTCTCTGGCa                                   <  1:390559/251‑1 (MQ=255)
                                                                                                                                                                                           ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAG‑CATACGATCGCAGCTAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAt                                  <  1:427944/251‑1 (MQ=255)
                                                                                                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCCGCGGTCCGG‑CAAAAAAGCCCGAGAAGGGATAGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGGAGTCTGCCCGAAGGCGTGGAGAAAGTTA‑CAA‑CATACGATCGCCGCCAGCTCATTGACCGCGCCATCGACAACC‑TCCGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGGAGTCTGTGCGCTGTTTCTCTGGCAtgt                                >  1:217751/1‑251 (MQ=255)
                                                                                                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAtgtg                               <  2:361429/251‑1 (MQ=255)
                                                                                                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTgc                              <  1:172156/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCt                           <  1:410716/251‑1 (MQ=255)
                                                                                                                                                                                                   aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTc                          <  1:86156/251‑1 (MQ=255)
                                                                                                                                                                                                      ggcgaaggcgaaGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGGAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCggg                       >  1:319503/1‑250 (MQ=255)
                                                                                                                                                                                                        cgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGGGTGGAGATAGTTA‑CAA‑CAAACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGGTTATTGTGGGGGCGGTAGGCTTTGCGCTGTTTTTCTTGCATGTGCGCTCGgcgc                     >  1:550843/1‑251 (MQ=255)
                                                                                                                                                                                                            ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCtggt                 <  2:421171/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgggcgggGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGattatt       >  1:948/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggcgggGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTAttt      >  1:156191/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggcgggGTGGAGATCCTGCGATCCGG‑CAAAAAAGCCCGAGCAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAACAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGTTCGGCGCTGGTCGCGATTAttt      >  2:459186/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt  <  1:475776/251‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
AACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTG‑AAGGACAAAC‑TGGAAACGC‑TGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTA‑CAA‑CATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC‑TCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTT  >  NC_000913/599353‑599822

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