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

AAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCC  >  NC_000913/599332‑599825
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
cgaaGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGAt                                                                                                                                                                                                                                                        <  2:121234/249‑1 (MQ=255)
  aaGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATcc                                                                                                                                                                                                                                                      <  1:66333/251‑1 (MQ=255)
         gtGGATGCTGCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAATCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATc                                                                                                                                                                                                                                               <  2:37763/250‑1 (MQ=255)
                     aaCCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGc                                                                                                                                                                                                                                  <  1:323154/251‑1 (MQ=255)
                                       tcgtcgATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGTGCCGGAGATGCGGCGGGGCATTGTCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATcgccg                                                                                                                                                                                                                 <  1:236353/250‑1 (MQ=255)
                                        cggcgatcgaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCgccgcc                                                                                                                                                                                                               <  2:381347/248‑1 (MQ=255)
                                                gaACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGG‑AAAAACACCCGCGAAGTGATCGCCGCCGTGAAggg                                                                                                                                                                                                       >  1:95768/1‑250 (MQ=255)
                                                       cggaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                                                                                                                                <  1:208018/251‑1 (MQ=255)
                                                            gcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGGGATCCGGC‑AAAAACGCCCGGGGAGGGATCGccgccggggaggggcaaatggaca                                                                                                                                                                                            >  1:219547/1‑232 (MQ=255)
                                                              ggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGCCCCGGAGATGCGCCGGGGCATTGCCGAACGAAACGGCGAAGGCTAAGTGGCGGGCGGTGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAATGACGAACTGGAAACGc                                                                                                                                                                                         <  2:334453/251‑1 (MQ=255)
                                                                    tatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGC‑AAAAACACCCGAGGAGTGATCGCCGCCGTGAAGGGCAAAATGGAAAAGCTGaaa                                                                                                                                                                                    >  1:34728/1‑250 (MQ=255)
                                                                      tGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTGGATCCGGC‑CAAAACGCCCGCGAAGGGCTCGCCGCCCGGAATGACAAACTGGAAAAGCTGAAAagg                                                                                                                                                                                 >  1:406646/1‑250 (MQ=255)
                                                                         tGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGTGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGTGATCCCGCGATCCCGC‑CAAAACCCCCGCGGAGGGATCCCCGCCGGGCAGGGAAAAATGGGAACGCTGGAAaggag‑t                                                                                                                                                                             >  1:205385/1‑251 (MQ=255)
                                                                             cgcCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGc                                                                                                                                                                         <  1:38886/251‑1 (MQ=255)
                                                                             cgcCAGCGGATATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGc                                                                                                                                                                         <  1:137995/251‑1 (MQ=255)
                                                                                    ggCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCTTTGGTCCGGC‑AAAAACTCCCAGGAAGGGTTTGCCGCCGCGAAGGAGAAACTTGAAACGCTGaaaaacagtactccggaagg                                                                                                                                                                   >  1:303728/1‑235 (MQ=255)
                                                                                            gCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTgg                                                                                                                                                             <  1:210314/248‑1 (MQ=255)
                                                                                              aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGGACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATTCGGC‑AAAAAAGACCGGGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGAAG‑GCTCCCGGAAGGCGGGgag                                                                                                                                                           >  1:78323/1‑248 (MQ=255)
                                                                                                  gCTCGACGACTTTAATCACAGCGTTTTAAAAGCCAGTGAAAAGGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCGAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTa                                                                                                                                                    <  2:416626/251‑1 (MQ=255)
                                                                                                      gacgacTTTAATCACATCGTTTTAAAAGCCAGAGAAAATGGCGTCCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGGGGGGGGGGGATCCTACGATTCGGC‑CAAAAAGCACGAGGAGAGGTCGCCGCCGTGAGGGGCAAACTGGAAACGCTGAGAAGAAG‑GATGCCGGAAGGGGGGGAGATCGttcgaa                                                                                                                                                 >  1:35378/1‑246 (MQ=255)
                                                                                                        cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTTGGGTCCGGG‑CAAAACCCCCGGGGAGGGGTCGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGAAG‑TCTGCCGGGAGGCgggggggtagtttaaacct                                                                                                                                              >  1:307573/1‑233 (MQ=255)
                                                                                                                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGCTGCGACGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCAGAAATAAACGCCGAAAGCGAAGTGGCGGGGGGGGGGGGGGTCCTACGGACCGGCCAAAAAC‑CCCGACGAGAGATCGGCGCCCTGCAGGGCAAACTGGAAACGCAGAAAAGGAGGTCGCCCGGAAG‑CGGGGAGATAGTTacaaaaatcgctcgccgacaacc                                                                                                                              >  1:406996/1‑232 (MQ=255)
                                                                                                                                aGCCAGTGAAAATGGCGTGCCCTTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGAGCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGGGGAGATATTTACAACATAC‑GATCGCAGCCAGCTCATTGAcc                                                                                                                     <  2:21706/251‑1 (MQ=255)
                                                                                                                                          aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGa                                                                                                           <  1:130831/251‑1 (MQ=255)
                                                                                                                                                          aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGCCCGCGATCCGGC‑AAAAACGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAACGCTGAAAAGGAG‑TCTGCCGGGAGGCGGGGGGGTAGTTACAACA‑ACAGATCGCAGCCCGCTCATTGACCGCGCCATCGACAACCTCCGCgggccg                                                                                           >  2:21985/1‑247 (MQ=255)
                                                                                                                                                             tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGATATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                                                                                                                                      <  1:374966/143‑1 (MQ=255)
                                                                                                                                                             tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCTCGATCCCGC‑CAAAACGCCCGGGGAGTGATTGCCGCCGGGGAGGACAAACTGGGAACGCTGAAAAGTAG‑TCTGCCGGAAGGGGGGGGGGTAGTTACAACAAAA‑GACCGCCGCCCGCCCCTTGACCGCGCCCTCGGCAACCTCAACGGCCAGTTg                                                                                        >  1:102924/1‑251 (MQ=255)
                                                                                                                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGt                                                                                          <  1:130049/248‑1 (MQ=255)
                                                                                                                                                                 cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTgg                                                                                    <  1:140689/251‑1 (MQ=255)
                                                                                                                                                                        gCGAAGGGCCAGATTGGCCCGGAGAGGCGCCGGGGCATAGCCGAACTAAACGGCAAAGGCGAAGTGGCGGGCGGGGTGTTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGtt                                                                             <  2:217509/251‑1 (MQ=255)
                                                                                                                                                                                 cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTgtggt                                                                    <  1:259799/251‑1 (MQ=255)
                                                                                                                                                                                           cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                                                                                                                                      <  2:322195/113‑1 (MQ=255)
                                                                                                                                                                                           cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                                                                                                                                      >  1:322195/1‑113 (MQ=255)
                                                                                                                                                                                                       ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGttt                                              <  1:50323/251‑1 (MQ=255)
                                                                                                                                                                                                       ggggCATTGCAGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGATCCGGC‑CAAAACGCACGGGGAGGGGTCGGCGCCGTGAAGGGCAAACTGGAAACGCTGTAAAGTAG‑TCTGGCGGAAGGCGGGGAGAGAGGTAAAACATAC‑GATCGAAGCCAGACCCTTGACCGCGCCATCGGCAACCTCAACGGCCAGCTGCTGGGAGAGTTTATTGTTGGGGCGGTAGTCGGTGCGCGGttt                                              >  2:331136/1‑251 (MQ=255)
                                                                                                                                                                                                           cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGAGAACCTGCGGTCCGGC‑AAAAACGCCCCAGAAGGGATCGCCGCCGTGAGGGGCAAACTGGATACGGTTAAAAGTCG‑TCTGCCGGGAGGGGTGGTGGGAGTTACAACATAA‑GATCGCCGCCAGCTCATTTACCGCGCCATCGGCCACCTCCGCGGCAAGGTGCTAGAGGGGGTTCTTTGGGTGGCGGGAGTCtgtgggcctgttcttt                                          >  1:68962/1‑239 (MQ=255)
                                                                                                                                                                                                                      tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGGGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGGAGAGTTTATTGTGGGGGCGGGAGTCTGTGCGCTGTTTCTCTGGCATGTgcgc                               >  1:267612/1‑251 (MQ=255)
                                                                                                                                                                                                                       aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCt                              <  2:102924/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGGTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGTTGGCg                <  2:380250/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ggcgggcgGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGatt             <  1:100743/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   gtggtgATCCTGCGCTCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAGGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGGGCGCTCGGCGCTGGTGGCGATTATTTCg       >  1:79517/1‑246 (MQ=255)
                                                                                                                                                                                                                                                   gtggtgATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAACCTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTTCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGTCATGTGCGCTCGGCGCTGGTGGCGATTATTTCg       <  2:79517/246‑1 (MQ=255)
                                                                                                                                                                                                                                                   gtggtgATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGcc  <  1:49062/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
AAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAG‑TCTGCCGGAAGGCGTGGAGATAGTTACAACATAC‑GATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCC  >  NC_000913/599332‑599825

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: ‑