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.

AACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTT  >  NC_000913/599353‑599822
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
aaCCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTGCGGCAAAAACGc                                                                                                                                                                                                                                   >  1:518966/1‑251 (MQ=255)
  ccAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                                                                                                                                 <  1:515982/251‑1 (MQ=255)
     ggAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                                                                                                                              <  1:335804/251‑1 (MQ=255)
         gcgggctgTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                                                                                                                          <  2:130576/251‑1 (MQ=255)
               ggTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGTGGGGTGGTGATCCTGCGAt                                                                                                                                                                                                                                                 <  1:338830/223‑1 (MQ=255)
                gTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                                                                                                                                   <  2:27098/251‑1 (MQ=255)
                        atcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTACCGCGATCCGGCAAAAAAGccccggaagggaaagccgccctga                                                                                                                                                                                                            >  1:340993/1‑229 (MQ=255)
                        atcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGTCCGGCAAAAAAGCCCCGGGAGGGATCGCCGCCGTgg                                                                                                                                                                                                            >  1:169152/1‑249 (MQ=255)
                                   ggaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                                                                                                                              <  2:222372/251‑1 (MQ=255)
                                          gaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                                                                                                                      <  2:452317/251‑1 (MQ=255)
                                            atatatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                                                                                                                                    <  2:17201/251‑1 (MQ=255)
                                                atGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGTCCCGCAAAAACGCCCCAGGAGGGGTCGCCGCCCTGCAGGGCAAACCGGGAACGcaggaaa                                                                                                                                                                                  >  1:100430/1‑244 (MQ=255)
                                                      cgcgcCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGGGCCTTCGGTCCGGCCAAAACACCCCAGAAGAGATCGCCGCCGTGGAGGGCAAACGTGAAAACCCGAAAAGTAGTCt                                                                                                                                                                          >  2:280783/1‑251 (MQ=255)
                                                            aGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGACCGGGAAAAACGCCCCAGAAATGATGGCCGCCGTGAAGGACAAAATGGAAAccgtgcaacattacggtaccgc                                                                                                                                                                     >  1:335091/1‑230 (MQ=255)
                                                             gCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAAAGCTGGAAAGGAGTCTGGCggg                                                                                                                                                                    >  1:234400/1‑249 (MQ=255)
                                                                gatATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCTGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                                                                                                                                <  2:268109/249‑1 (MQ=255)
                                                                  tATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGGTCCCGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAAAGCTGGAAAGGAGTCTGCCGGGAggggg                                                                                                                                                              >  1:32560/1‑248 (MQ=255)
                                                                       gCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAAAGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCGCTTCGGGACAAACGCCCGGGAAAGGAACGCCGCCGGGCAGGGAAAACTGGAAACGACGGAAAGGAGTCTGCCGGGAGGGGTggtgg                                                                                                                                                         >  1:230061/1‑248 (MQ=255)
                                                                               tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTGCGGGAAAAACCCCCCCGAAGGGAACGGCGGCGTGGAGGGACAACTGGAAACGCTGACAAGCAACCTGCCGGGAGGCGGGGGGATAGATacc                                                                                                                                                >  1:116788/1‑250 (MQ=255)
                                                                                   cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                                                                                                                           <  2:343328/251‑1 (MQ=255)
                                                                                        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGAGGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                                                                                                                        <  1:139288/249‑1 (MQ=255)
                                                                                          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                                                                                                                     <  2:116788/250‑1 (MQ=255)
                                                                                               caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGGGGGGGGGGGGGGGGACCTGCGGACCGGGCAAAACGCCCGAGAAGTGGTCGCCGCCCGGCAGGACAAACTGGAAAAGCTGAAAAAGAGTATGCCGGGAGGCGGGGAGAGAGTTAAAACATACGATCGcagcc                                                                                                                               >  1:96583/1‑251 (MQ=255)
                                                                                                         aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGTCCCGCGAGCCGGCCAAAACGCCCGAGAAGCGGCAGCCCCCCGGAAGGGCAAAATGGGAAAGCTGAAAAAAAGTCTGCCCGGAGGGGTGGGGATAGTTAAAAACAAACATCGCAACCAGCTCCTTgt                                                                                                                     >  2:86156/1‑250 (MQ=255)
                                                                                                            gCCAGTGAAATGGGCGTGACCTTGTATCGGCGCGATGTTGCAAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACAAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCAATCCGGCAAAAACGCACGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTAATTGAcct                                                                                                                  <  2:76315/251‑2 (MQ=255)
                                                                                                                gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTGCCGCGATCCCGGAAAAACGCCCGGGAAGGGATCGCCGCCGGGGAGGAAAAACTGGAAACGCTGAAAAGGAGTTTGGCGGAAGGCGGGGAGGTAGTTACCACATACGAACGCAGCCAGCTCCTTGAACGCGcc                                                                                                              >  1:421171/1‑251 (MQ=255)
                                                                                                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt                                                                                                            <  1:292651/251‑1 (MQ=255)
                                                                                                                   aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc                                                                                                           <  1:403467/251‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCGGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTAGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc                                                                                                           <  2:353586/248‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc                                                                                                        <  1:203323/251‑1 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGTGAAGGGCAAACTAGAAAAGCTGAAAAGTAGTCTGCCGGGAGGGGTGGGGGTAGTTACAACATACGAACGCAGCCCGCTcctttacccccccctcaa                                                                                                         >  1:353586/1‑233 (MQ=255)
                                                                                                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGGAAAAACGCCCCAGAAGGGATCGCCGCCGGGAAGGGCAAACTGGAAACGCTGGAAAGAAGTCTGCCGGGAGGCGGGGGGATAGGTAAAACATACGAACGCAGCCCGCTCCTTGACCGCGCCCTCggc                                                                                                        >  1:228449/1‑249 (MQ=255)
                                                                                                                          cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGCGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaaac                                                                                                    >  1:361429/1‑249 (MQ=255)
                                                                                                                                 gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATCCGGCAAAAACGCCCGAGAAGTGGTCGCCGCCGGGACAGACAAAATGGGAAAGATGAAAAGAAGTCTGCCGGCAGGGGGGGGGAGATATACCAAATACGATAGCAGCCAGCTCATTGACCGCGCCATCGCAACCCTCAGCgg                                                                                            >  2:44610/1‑251 (MQ=255)
                                                                                                                                     aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTCCGATCCGGCCAAAACGCCCGAGAAGGGATGGCCGCCGCGAAGGACAAAATGGTAACACAGGCAAAGAGGCTGGCGGAAGGCGGGGGGAGAGTTAACACCAACGATCCCAGCCAGCTCCTTGAACGCGCCCTCCGCAACCTCCGCGGCa                                                                                          >  2:66588/1‑249 (MQ=255)
                                                                                                                                     aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGGAGGGCAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAAGGGGGGGGATAGGTACAACATACGAACGCAGCCAGCTCCTTTACCGCGCCCTCGAGACCCCCAGCgggc                                                                                          >  2:477623/1‑247 (MQ=255)
                                                                                                                                        tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCGGACCGGCAAAAAAGGCCGCGAAGTGATCGCCGCCGCGAAGGACAAACTGGAAAACCCGAAAAGAAGGATGGCGGAAGGGGGTGAGATAGTTACAAAATACGATCCCCGCCCGCCCATTTGCCGCGCCATCCACAACCTCAGCGGcccgttg                                                                                     >  2:403016/1‑245 (MQ=255)
                                                                                                                                          cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGGCAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGGCAAAATGGAAACGCCCAAAAGGAGTCTGGCGGGAGGCGGGGGGGTAGTTACAAAAAACGACAGCAGCCCGCCCATTGAACGCGCCacccaacacctcccccgggagtggcg                                                                                    >  2:295264/1‑226 (MQ=255)
                                                                                                                                                            cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTgtggt                                                                 <  1:120410/251‑1 (MQ=255)
                                                                                                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGc                                                              <  2:169152/251‑1 (MQ=255)
                                                                                                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAACGCCCGCGAAGGGATCGCCGCCGTGAAGGGCAAAACGGAAACGCTGAAAAGTAGTCTGCCGGGAGGGGTGGGGGTAGTTAAAACATACGACGGCCGCCAGCCCATTTACCGCGCCCTCGACCACCCCAACGGCAAAGTGGTGGAAAAGGTTATTGtggtggg                                                              >  2:219911/1‑250 (MQ=255)
                                                                                                                                                                    cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAg                                                         <  1:280783/251‑1 (MQ=255)
                                                                                                                                                                       ggAGATGCGCCGGGGCCTTGCCGACCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAACCGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAATGCGTGGAGATAGTTACAACATCCGATCGCAGCCAGCTCATTGACCGCGCCATCGACGACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCt                                                      <  1:261918/251‑1 (MQ=255)
                                                                                                                                                                           aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCCGGCAAAAAAGCCCGAGAAGGGAGCGCCGCCGGGAAGGACAAACCGGAAAAGCGGAAAAGGAGGCTGCCGGAAGGCGGGGGGATAGGTACAACAAACGGTCGCCGCCAGCTCATTTACCGCGCCCCCGACAAACTCCGCGGCCAGGTGCTGGGAGGGTTTTTTGTGGGGGCgggaggctgggg                                                  >  2:51484/1‑241 (MQ=255)
                                                                                                                                                                               gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGAACCCGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGGCAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGAAAGTTACAAAATACGAACGCCGCCCGCTCATTTACCGCGCCCTCGACAACCCCAGCGGCAAGGTGTTGGGAGGGTTTTTTTTGGGGGCgggggggtcgggggcg                                              >  2:335867/1‑237 (MQ=255)
                                                                                                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGAGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTgg                                     <  2:41067/250‑1 (MQ=255)
                                                                                                                                                                                          gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCCGTTTCTCTGGCa                                   <  1:390559/251‑1 (MQ=255)
                                                                                                                                                                                           ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAGCATACGATCGCAGCTAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAt                                  <  1:427944/251‑1 (MQ=255)
                                                                                                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCCGCGGTCCGGCAAAAAAGCCCGAGAAGGGATAGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCCGAAGGCGTGGAGAAAGTTACAACATACGATCGCCGCCAGCTCATTGACCGCGCCATCGACAACCTCCGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGGAGTCTGTGCGCTGTTTCTCTGGCAtgt                                >  1:217751/1‑251 (MQ=255)
                                                                                                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAtgtg                               <  2:361429/251‑1 (MQ=255)
                                                                                                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTgc                              <  1:172156/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCt                           <  1:410716/251‑1 (MQ=255)
                                                                                                                                                                                                   aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTc                          <  1:86156/251‑1 (MQ=255)
                                                                                                                                                                                                      ggcgaaggcgaaGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGGAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCggg                       >  1:319503/1‑250 (MQ=255)
                                                                                                                                                                                                        cgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGGGTGGAGATAGTTACAACAAACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGGTTATTGTGGGGGCGGTAGGCTTTGCGCTGTTTTTCTTGCATGTGCGCTCGgcgc                     >  1:550843/1‑251 (MQ=255)
                                                                                                                                                                                                            ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCtggt                 <  2:421171/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGattatt       >  1:948/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTAttt      >  1:156191/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggcgggGTGGAGATCCTGCGATCCGGCAAAAAAGCCCGAGCAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAACAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGTTCGGCGCTGGTCGCGATTAttt      >  2:459186/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt  <  1:475776/251‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
AACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTT  >  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: