Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG41.4% 148.8 / 38.3 87V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (4/47);  new base G (36/0);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 2.50e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599426‑599726
                                                                                                                                                      |                                                                                                                                                         
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                           <  2:510452/151‑1 (MQ=255)
    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                                                                                        <  2:1037437/150‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                      <  1:374109/149‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGAt                                                                                                                                                     <  2:714907/150‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                                                                                                                                                    <  1:694369/151‑1 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                                  <  1:570784/151‑1 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                                  <  1:440060/151‑1 (MQ=255)
          cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                                   >  2:647338/1‑149 (MQ=255)
               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTccctcggt                                                                                                                                             >  2:234768/1‑144 (MQ=255)
               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGa                                                                                                                                              <  2:984673/149‑1 (MQ=255)
                tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑c                                                                                                                                           <  2:321184/150‑1 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCgg                                                                                                                                        <  2:1080381/151‑1 (MQ=255)
                      caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCgg                                                                                                                                        >  1:1014744/1‑147 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCaaaa                                                                                                                                   <  1:1039046/151‑1 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGT‑TCgg                                                                                                                                        >  2:719280/1‑146 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCaaaa                                                                                                                                   <  2:582443/149‑1 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAAc                                                                                                                                 <  2:986005/151‑1 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGGT‑CCGGCaaaaa                                                                                                                                  >  1:416418/1‑150 (MQ=255)
                           ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGAT‑CCGGCaaaaa                                                                                                                                  >  1:72937/1‑148 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTTCGGTCCCGGCaaaaaaccc                                                                                                                              >  2:737581/1‑145 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGT‑TCCGGAAAAAAGCCCgag                                                                                                                          >  2:13750/1‑149 (MQ=255)
                                   gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGT‑CCGGCAAAAACGccc                                                                                                                             >  1:668514/1‑145 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGt                                                                                                                      <  1:308282/151‑1 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGt                                                                                                                      <  2:1153894/151‑1 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCCGT‑CCGGGCAAAAAGccc                                                                                                                             >  2:1058553/1‑144 (MQ=255)
                                      aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGa                                                                                                                    <  1:343730/151‑1 (MQ=255)
                                        tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGATCCTGCGGA‑CCGGCAAAAACGCCCGGGAAGTGATc                                                                                                                  >  1:403434/1‑151 (MQ=255)
                                        tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGGTCCTGCGGT‑CCGGCAAAAAACCCCGAGAAGGGAt                                                                                                                   >  1:486903/1‑150 (MQ=255)
                                        tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTg                                                                                                                     >  1:15269/1‑148 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGAT‑CCGGCCAAAAACCCCgcgga                                                                                                                        >  2:1169639/1‑140 (MQ=255)
                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCCGCAAAAAAGCCCgaggaggg                                                                                                                     >  1:972378/1‑144 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCCAAAAAGCCCCAGAAGTGATcgc‑c                                                                                                              >  1:643190/1‑151 (MQ=255)
                                             aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTTCGGTCCCGGAAAAAACCCCCggaagggagaggc                                                                                                                >  2:58163/1‑137 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgc‑cgcc                                                                                                           <  1:1007746/150‑1 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgc‑cgcc                                                                                                           <  2:977859/150‑1 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgc‑cgcc                                                                                                           <  2:710423/150‑1 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGGAAAAACGCCCGAGAAGGGATCgc‑cgcc                                                                                                           >  1:303239/1‑150 (MQ=255)
                                                gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑Gt                                                                                                        <  2:606948/151‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGaa                                                                                                     <  1:805935/150‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAg                                                                                                    <  1:723108/151‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAACCTTGGGT‑CCGGGGAAAACAGCCGGGAAGTGGTcgc‑cgc                                                                                                            >  2:611720/1‑144 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGC‑CGCC‑GCGAAgg                                                                                                   >  2:1093014/1‑150 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGATCGC‑CGCC‑GCGAAgg                                                                                                   >  2:368234/1‑150 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGC‑CGCC‑GTGGAgg                                                                                                   >  1:693299/1‑150 (MQ=255)
                                                       cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGAc                                                                                                 <  1:1173063/151‑1 (MQ=255)
                                                           tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGGAGGGGTCGC‑CGCC‑GCGAAgggaaaa                                                                                              >  2:720591/1‑145 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAAc                                                                                             <  1:643321/150‑1 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGGAAGGATCGC‑CGCC‑GGGGAGGACaaa                                                                                              >  2:628887/1‑149 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGaaa                                                                                       <  2:698157/151‑1 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑TCGGCGAAAAACCCCGCGAAGGGGTCGC‑CGCC‑GTGAAGGGCAAAATGGaa                                                                                        >  2:1011495/1‑150 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCCGCCAAAAAGCCCGGGGAGGGAATGC‑CGCCGGGGAGGGACAAA‑TGGaa                                                                                        >  2:227943/1‑150 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGC‑CGCC‑GTGAAGGACAAACTGGaaa                                                                                       >  1:827157/1‑147 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGAT‑CCCGCAAAAAAGCCCGCGGAGGGGTCGC‑CGCC‑GGGGAGGACAAACTGGGAACGcggga                                                                                >  2:48795/1‑146 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGGTCCTCCGGT‑CCGGCCAAAAACCCCCACAAGAGATCGC‑CGCC‑GTCAAGGACAAACTGGAAACGCCGaaaa                                                                              >  2:245674/1‑150 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAag                                                                             <  2:1122884/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGaaaa                                                                              <  1:283278/149‑1 (MQ=255)
                                                                                gtacAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGGTCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAagtagt                                                                         <  1:527455/147‑1 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                        <  2:36925/151‑1 (MQ=255)
                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                       <  2:447289/150‑1 (MQ=255)
                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                    <  2:297983/151‑1 (MQ=255)
                                                                                        tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGCG‑CCGGCCAAAAAGCCCGGGGAGGGGTCGC‑CGCC‑GTGGAGGGCAAAACGGGAACGCCGAAAAGAAGTCTGCCgg                                                                  >  2:510648/1‑149 (MQ=255)
                                                                                        tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCCAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCgg                                                                  >  1:955270/1‑149 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                                <  1:79552/150‑1 (MQ=255)
                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGC‑CGCC‑GGGAAGGACAAAATGGGAACGCTGAAAAATAGTCTGCCGGGAGGCg                                                            >  2:1020766/1‑150 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCCGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGGAGGCggg                                                          >  1:708289/1‑149 (MQ=255)
                                                                                                agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                        <  1:296620/151‑1 (MQ=255)
                                                                                                 gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                       <  2:1018363/151‑1 (MQ=255)
                                                                                                  aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                      <  1:822088/151‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                                  <  2:233586/151‑1 (MQ=255)
                                                                                                        cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                                <  1:1159587/151‑1 (MQ=255)
                                                                                                        cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCCC‑CGCC‑GTGAAGGGCAAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGGGGAGATa                                                    >  1:281425/1‑147 (MQ=255)
                                                                                                         ggggCATTGCCGAACTAAACGGCGAAGGCGCAGTGGCGGGCGGGGGGGTGATTACGCGAA‑ACCGGAAAAACGCCCGGGAAATGATCACAGGCC‑GGTAGGGAC‑AACCGGAAACGCCGGAAAGGAGTCAGCCGGGAGGCGGGGAGGTa                                                    >  1:880789/1‑146 (MQ=255)
                                                                                                            gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                            <  1:1099876/151‑1 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGTGGGGGGGGGATTCTTTCGG‑CCGGCAAAAAAAACCGAGGAGAGGGCGC‑CGCC‑GGGAAGGACAAAAAAGAAAAACTGAAAAGAAGTCTGTCCGGAGGCGGGGGGATAATTacaaac                                           >  2:562073/1‑146 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑CCGGCAAAAAAGCCCCGGAAGGGATCGG‑CGCC‑GTGGAGGGCAAAATGGGAAAGCTGAAAAAGAGTCTGCCGGGAAGGGGGGGGATAGTTacaaca                                           >  2:854516/1‑148 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                        <  2:869291/150‑1 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                       <  1:287750/151‑1 (MQ=255)
                                                                                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                    <  2:503504/150‑1 (MQ=255)
                                                                                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                   <  1:283094/151‑1 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                         <  2:532079/150‑1 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                        <  2:1164062/151‑1 (MQ=255)
                                                                                                                                    gcgaagTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGGTCGC‑CGCC‑GGGGAGGACAAACTGGAAACGCTGAAAAGAAGGCTGCCGGAAGGCGGGGGGATAGTTACAACATAAGAACCCAGCCAGCTCCTTg                     >  2:305705/1‑150 (MQ=255)
                                                                                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                  <  1:716051/151‑1 (MQ=255)
                                                                                                                                              gggcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc          <  1:110949/151‑1 (MQ=255)
                                                                                                                                              gggcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc          <  2:1083467/151‑1 (MQ=255)
                                                                                                                                               ggcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg         <  1:386206/151‑1 (MQ=255)
                                                                                                                                                gcgggGTGGAGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTAATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCGGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGTGCCATCGa        <  2:516621/151‑1 (MQ=255)
                                                                                                                                                 cgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc       >  1:68807/1‑151 (MQ=255)
                                                                                                                                                   gggTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa      >  1:768245/1‑150 (MQ=255)
                                                                                                                                                      tggtgATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  1:882290/1‑151 (MQ=255)
                                                                                                                                                      tggagATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc    >  2:400276/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                         
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGC‑CGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599426‑599726

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