Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG46.8% 104.8 / 34.8 80V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (5/37);  new base G (35/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.20e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725
                                                                                                                                                        |                                                                                                                                                     
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGggt                                                                                                                                                       <  1:705242/151‑1 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑Ggggggg                                                                                                                                                     >  1:465311/1‑150 (MQ=255)
    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑Gggggggg                                                                                                                                                    >  1:1423957/1‑149 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATc                                                                                                                                                <  2:828998/151‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATc                                                                                                                                                <  2:624493/151‑1 (MQ=255)
        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATcc                                                                                                                                               >  1:297705/1‑151 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTCCTg                                                                                                                                             >  1:242654/1‑150 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGc                                                                                                                                            <  2:649171/151‑1 (MQ=255)
             cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGTCCTTCg                                                                                                                                           >  2:1119493/1‑150 (MQ=255)
             cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGGCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGggggtgggcc                                                                                                                                               >  2:814152/1‑142 (MQ=255)
              tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCg                                                                                                                                           <  1:273639/149‑1 (MQ=255)
                tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGAA‑GGGGGGGGGagtctacgacc                                                                                                                                        >  1:934425/1‑140 (MQ=255)
                 aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTcccccgaacc                                                                                                                                       >  2:595117/1‑142 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCgg                                                                                                                                     <  2:33725/151‑1 (MQ=255)
                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGc                                                                                                                                    <  1:1014546/150‑1 (MQ=255)
                     acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCaa                                                                                                                                  <  1:1437870/151‑1 (MQ=255)
                      caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCaaa                                                                                                                                 <  1:1086297/151‑1 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCaaaa                                                                                                                                <  1:1199436/151‑1 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGTCCttaggccccggcaa                                                                                                                                 >  2:1429228/1‑137 (MQ=255)
                          gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                             <  1:1446620/151‑1 (MQ=255)
                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTCCTTCGGTCCGGCAAAAACGccc                                                                                                                          >  2:666971/1‑151 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                          <  1:677872/150‑1 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                          >  1:1187807/1‑149 (MQ=255)
                                aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGC‑GGGGGGGGGTCCCTCGGTCCCGGCAAAACCGcc                                                                                                                           >  2:702327/1‑148 (MQ=255)
                                aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCCAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGT‑GTGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                       <  2:938817/151‑1 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                   <  1:614841/151‑1 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGGCCCGCGATCCCGCCAAAAAGCCCgaga                                                                                                                      >  1:1365437/1‑148 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                       >  1:92689/1‑147 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGAAGTg                                                                                                                  >  1:273836/1‑151 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGA‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCCAGAAGGGATcg                                                                                                              >  1:1407645/1‑151 (MQ=255)
                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               >  1:111771/1‑149 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTCCTGGGATCCGGGAAAAAAGCCCggggagggaatg                                                                                                              >  2:125823/1‑138 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                            <  2:1105872/151‑1 (MQ=255)
                                            aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTTCTGCGATCCGGGAAAAAAGCCCgaggaggg                                                                                                                  >  1:453661/1‑142 (MQ=255)
                                                     cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGATCC‑GCGCTCCGGCAAAAAAACACGAGAAGTGTGCGccggcggg                                                                                                      >  2:1244090/1‑145 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGTGGGGTTCCTTCGATCCGGCAAAAAAGCCCGAGGAGGGGTCGGCGGCGGGAAgg                                                                                                  >  2:340150/1‑150 (MQ=255)
                                                       cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCACAAACGCCCGAGAAGTGATTgccgcc                                                                                                         >  2:123557/1‑142 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGGGATCCTGCTCTCCGGCAAAAACGCCCCGGAAGCGATGGGCGCCGTGAAGGAc                                                                                                >  1:534065/1‑150 (MQ=255)
                                                           tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCCAGAAGAGATCTCCGCCGGGAAGGAaaaa                                                                                             >  1:1088678/1‑150 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGCAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGCTCCGGCCAAAAAGCCCGCGAAGTGGATGCCGCCCTGAAGGACAAAc                                                                                            >  1:491067/1‑150 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                           <  1:1398867/147‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGTCCTGCGGTCCCGCAAAAAAGGCCGAGAAGTGGACGCCGGCCGGAAGGGGAAAcc                                                                                           >  2:621010/1‑146 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑CGGGCGGTGATCCTGCGGTCCCGCAAAAAACCCCCAGAGCGGGTCCCCCCCGGGAAGGACAAACtcg                                                                                         >  2:179219/1‑146 (MQ=255)
                                                                  gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGGCAAAAAGCCCGAAAAGGGATAGCCGCCCGGAAGGACAAACTGGaaa                                                                                      >  1:1379443/1‑150 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc                                                                                   <  2:372506/151‑1 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCCGGAAAAACGCCCGCGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGc                                                                                   >  1:972618/1‑150 (MQ=255)
                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGG‑GGGGGGGGGGTCCCGCGATCCGGCAAAAAACCCCGAGAAGGGATCGGCGCCGTGGAGGACAAACTGGAAAAGc                                                                                   >  2:1044690/1‑149 (MQ=255)
                                                                        ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                                <  1:20350/150‑1 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                               <  1:118900/150‑1 (MQ=255)
                                                                           cGAAGGTCCCGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGACGTGGC‑GGGC‑CGGGGGGTGGCCCTGGGGTGCGGCCAAAACACACCAGAAGAGATCGCCGCCCGGGAGGGCAAAACGGAAAAGCTGaaa                                                                              >  2:905469/1‑149 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCCTGAAGGACAAACTGGAAACGCTGaaa                                                                              >  1:670476/1‑149 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCC‑GGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGAGCCTGGGATCCGGCAAAAAAGCCCGAGAAGGGAACGCCGCCGTGAAGGGCAAAAAGGGAACCGTGaaa                                                                              >  1:1178368/1‑148 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:71536/151‑1 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                          <  1:517069/151‑1 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGAACGCCCCCGTGGAGGACAAACTGGAAAAGCTGGAAag                                                                            >  2:1281141/1‑149 (MQ=255)
                                                                              aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                         <  1:1059934/151‑1 (MQ=255)
                                                                                 tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                       >  1:1387122/1‑150 (MQ=255)
                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                     <  1:730301/151‑1 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                              <  1:1357151/151‑1 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCCGAAAAGTAGTCTGCCGGa                                                                >  2:1304767/1‑149 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  1:766133/151‑1 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  2:465148/150‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATa                                                   >  1:639011/1‑149 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGGCTGCCGGGAGGCGGGGGGGTAGt                                                 >  1:675222/1‑151 (MQ=255)
                                                                                                           ggaaTTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGGGCCTGCGATCCGGGAAAAACGCCCGCGAAGGGGTCGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAAGTAGTTTGCCGGGAGGCGTGGGGATAGTTaaa                                             >  2:1294579/4‑148 (MQ=255)
                                                                                                             cATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                          <  2:560697/151‑1 (MQ=255)
                                                                                                             cATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCCAGAAGGGATCGCCGCCGTGTAGGGCAAAATGGGAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGATAATTtcaac                                           >  1:168618/1‑150 (MQ=255)
                                                                                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGCGAAGGACAAAATGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGGGGGGAGATAGTTacaaca                                          >  1:944126/1‑150 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGGGGTGTTCCTGCGCTCCGGCGAAAACGCCCCAGAAGAGATCCGCGCCGTGAAGGACAAACCGGAAACGCTGAAAAGAAGGCTGCCGGAAGGCGGGGCGGTAGGTacaacaaa                                        >  1:667890/1‑148 (MQ=255)
                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  1:595302/151‑1 (MQ=255)
                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  1:353522/151‑1 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                    <  1:680331/145‑1 (MQ=255)
                                                                                                                         aaaCGGCGAAGACGAGGGGGC‑GGGA‑GGGGGGGTGATCCTGCGATCCGGCAGAAGCGCCCGAGAAGTGATCGGCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTATAGATAGTTACAGCATACGATCGca                                <  2:852902/149‑1 (MQ=255)
                                                                                                                           aCGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                            <  1:21218/151‑1 (MQ=255)
                                                                                                                             ggcgaaggcgaaGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccagc                          <  2:1101566/151‑1 (MQ=255)
                                                                                                                             ggcgaaggcgaaGTGGC‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATCGcagccag                           >  1:1040915/1‑150 (MQ=255)
                                                                                                                                 aaggcgaagTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                      <  1:360307/151‑1 (MQ=255)
                                                                                                                                   ggcgaagTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                    <  1:1247630/151‑1 (MQ=255)
                                                                                                                                     cgaagTGGG‑GGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAGCGCCCGAGAAGTGATGGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                  <  2:1249757/151‑1 (MQ=255)
                                                                                                                                     cgaagTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                  <  1:81381/151‑1 (MQ=255)
                                                                                                                                           ggc‑gggc‑gGGGTGGTGACCCTGCGATCCGGCAAAACCGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             <  2:1178009/150‑1 (MQ=255)
                                                                                                                                             c‑gggc‑gggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  2:111751/151‑1 (MQ=255)
                                                                                                                                               gggc‑gggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  2:187789/150‑1 (MQ=255)
                                                                                                                                                     gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     >  1:745248/1‑150 (MQ=255)
                                                                                                                                                        tggtgATCCTGCGATCCGGCAAAACCGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  1:397838/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725

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