Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG51.3% 99.1 / 132.8 76V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (39/0);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 1.21e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-12

ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAA  >  NC_000913/599428‑599723
                                                                                                                                                     |                                                                                                                                                    
aCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggg                                                                                                                                                    >  2:577012/1‑151 (MQ=255)
aCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggg                                                                                                                                                    >  2:813210/1‑151 (MQ=255)
   cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggg‑ggggggggg                                                                                                                                                  >  1:1095159/1‑148 (MQ=255)
          acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGc                                                                                                                                           >  1:724583/1‑150 (MQ=255)
                    caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGCGATACTGCGGTCC‑GGCaa                                                                                                                                >  2:1091401/1‑150 (MQ=255)
                       cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAAc                                                                                                                            <  1:1164169/151‑1 (MQ=255)
                       cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTTCTGCGGTCC‑GGCcaaaa                                                                                                                             >  2:645359/1‑150 (MQ=255)
                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGATCC‑CGCaaaaa                                                                                                                             >  2:819523/1‑149 (MQ=255)
                         ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGGTCC‑GGCAAAAAAGc                                                                                                                          >  2:235487/1‑151 (MQ=255)
                             aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCga                                                                                                                      <  1:240903/151‑1 (MQ=255)
                              aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGGTCC‑GGCCAAAAAGccc                                                                                                                        >  1:827195/1‑148 (MQ=255)
                              aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGAGCCTGCGATCC‑GGCAAAAAAGccc                                                                                                                        >  2:967020/1‑148 (MQ=255)
                               aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCgcg                                                                                                                     >  2:19854/1‑148 (MQ=255)
                                  ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCgag                                                                                                                     >  1:341335/1‑147 (MQ=255)
                                    aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGTCCTGGGAACC‑GGCCAAAACGcccccgaa                                                                                                                   >  1:1041632/1‑142 (MQ=255)
                                     gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGAt                                                                                                              <  2:695455/151‑1 (MQ=255)
                                      tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGAGTCTGCGATCC‑GGGAAAAAAGCCCggggaggg                                                                                                                >  2:160598/1‑141 (MQ=255)
                                       gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGTGTTCGT‑CGCACC‑CGCAAAAACGCCCCCGAAggg                                                                                                                >  2:1107113/1‑145 (MQ=255)
                                        aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATcgc                                                                                                           <  2:407041/151‑1 (MQ=255)
                                         aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑CGCAAAAACGCCCGCGAAgaga                                                                                                               >  2:412040/1‑146 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGATCC‑GGCCAAAACGCCCGGGAAGGGATcgcc                                                                                                          >  1:742842/1‑147 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                       <  2:532631/150‑1 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                       <  2:354657/150‑1 (MQ=255)
                                              gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGt                                                                                                     <  2:33976/151‑1 (MQ=255)
                                                gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGTCCCGCGGGCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTg                                                                                                    >  1:532498/1‑150 (MQ=255)
                                                gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGGGAAGGGATCgccgcc                                                                                                       >  2:730391/1‑147 (MQ=255)
                                                 tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGAGCCGCGATCT‑GGCAAAAACCACCGAGGAGTGGGcgccggccgcg                                                                                                   >  1:350668/1‑144 (MQ=255)
                                                    ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGTCCTTCGATCC‑GGCAAAAACGCCCGAGAAGTgcacggccccgcgaa                                                                                                  >  2:767357/1‑134 (MQ=255)
                                                      gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGGGATCCTCCGGTCC‑GGCTCAAAAGCCCGAGAAGTGATAgccgcc                                                                                                       >  1:645043/1‑141 (MQ=255)
                                                        ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGTGATCCTGGGACCC‑GGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGaa                                                                                                  >  2:939747/1‑144 (MQ=255)
                                                         tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                           <  2:999851/150‑1 (MQ=255)
                                                          aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                          <  1:164534/150‑1 (MQ=255)
                                                          aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATTGCCGCCGCGAAGGACaaa                                                                                           >  2:345166/1‑149 (MQ=255)
                                                           tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGGCAAAAAACCCCG‑GAAGGGGGCGCCGCCGGGGAGGACAAAATg                                                                                        >  1:532183/1‑151 (MQ=255)
                                                             tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGGTCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGCGAAGGGCAAACTgg                                                                                       >  2:511920/1‑150 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                         <  2:402101/147‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATACTGGGGTCC‑GGCCAAAACGCCCGAGGAGTGATCGCCGGCGGGAAGGAGAAAATgg                                                                                       >  2:424879/1‑149 (MQ=255)
                                                               cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:864255/151‑1 (MQ=255)
                                                                gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑CGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    >  2:1172624/1‑150 (MQ=255)
                                                                 cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                   <  1:377894/150‑1 (MQ=255)
                                                                   aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                <  2:1166670/151‑1 (MQ=255)
                                                                   aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                <  2:419715/151‑1 (MQ=255)
                                                                   aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                <  1:302139/151‑1 (MQ=255)
                                                                   aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAAATGGaaa                                                                                    >  1:776730/1‑147 (MQ=255)
                                                                    tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCCGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACCAACTGGAAAAGCTg                                                                               >  2:994860/1‑151 (MQ=255)
                                                                    tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCCGCGCTCC‑GGCAAAAACACCCGGGAAGCGGTTCCCCCCGCGCAGGGCAAACTGGgaa                                                                                    >  1:114992/1‑146 (MQ=255)
                                                                       tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCGATCC‑GGCCAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                             >  2:19497/1‑150 (MQ=255)
                                                                         cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                          <  2:1142216/151‑1 (MQ=255)
                                                                         cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCC‑GGGAAAAAAGCCCGAGAAATGATCGCCGCCGGGGAGGACAAACTTGAAAAGCTGaaaa                                                                           >  2:681670/1‑150 (MQ=255)
                                                                                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                    <  1:979571/150‑1 (MQ=255)
                                                                                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                    <  1:399393/150‑1 (MQ=255)
                                                                                   gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGGGAAGTGGTCGCCCCCCTGGAGGACAAACTGGAAAAGCTGGAAAGGAGTCTGcc                                                                 >  1:442413/1‑150 (MQ=255)
                                                                                    aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                               <  2:1112471/151‑1 (MQ=255)
                                                                                        gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                           >  2:400621/1‑151 (MQ=255)
                                                                                        gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGGAAAAACGGCCGGGAAGGGATCGCCGCCGGGAAGGGCAAAATGGAAACGCTGAAAAGGAGTCTGCCGGAAgg                                                           >  2:424728/1‑151 (MQ=255)
                                                                                         cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                           >  2:1004037/1‑150 (MQ=255)
                                                                                              agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCGCGGTCC‑GGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGGAGGACAAAATGGGAACGCCGGAAAGGAGGCTGCCGGGAGGGGTgg                                                      >  1:231179/1‑150 (MQ=255)
                                                                                               gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                     <  2:998385/150‑1 (MQ=255)
                                                                                                   cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGGCCTGCGCTCC‑CGCGAAAACGCCCGAGAAGGGATCCCCGCCGCGGAGTACAAACATGAAACGCGGAAAAACAGTCTGCCGGGAAGCgggggg                                                    >  1:1093152/1‑145 (MQ=255)
                                                                                                    gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                                <  2:1158970/150‑1 (MQ=255)
                                                                                                           cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                          >  2:340004/1‑149 (MQ=255)
                                                                                                            aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                        >  2:1020837/1‑150 (MQ=255)
                                                                                                               gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGGTCCTGTGATCC‑GGGAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGGAGGATGCCGGGAGGGGGGgagggaggttaaaaatac                                     >  1:666766/1‑135 (MQ=255)
                                                                                                               gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                    <  1:142059/151‑1 (MQ=255)
                                                                                                               gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                      >  2:887245/1‑149 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGGGGGGATAGTTACAACATACg                                    >  2:733055/1‑149 (MQ=255)
                                                                                                                  gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                 <  2:341339/151‑1 (MQ=255)
                                                                                                                         aCGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                          <  1:813194/151‑1 (MQ=255)
                                                                                                                            gcgaaggcgaagTGGCGGGC‑GGGGGGGGGATCCCGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                          >  2:910951/1‑148 (MQ=255)
                                                                                                                             cgaaggcgaagTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                      <  1:223013/151‑1 (MQ=255)
                                                                                                                                  gcgaagTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGa                 <  1:900584/151‑1 (MQ=255)
                                                                                                                                    gaagTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc               <  2:249521/151‑1 (MQ=255)
                                                                                                                                         ggcgggc‑gGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc           <  1:290018/150‑1 (MQ=255)
                                                                                                                                           cgggc‑gggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt        >  1:213631/1‑151 (MQ=255)
                                                                                                                                            gggc‑gggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc       <  1:970070/151‑1 (MQ=255)
                                                                                                                                            gggc‑gggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc       <  1:852092/151‑1 (MQ=255)
                                                                                                                                            gggc‑gggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc       >  1:425736/1‑151 (MQ=255)
                                                                                                                                                  gggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa  <  2:12515/151‑1 (MQ=255)
                                                                                                                                                  gggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa   >  2:79556/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAA  >  NC_000913/599428‑599723

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