Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG43.2% 95.6 / 39.0 81V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (6/40);  new base G (35/0);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 5.06e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722
                                                                                                                                                      |                                                                                                                                                    
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                      <  1:48107/151‑1 (MQ=255)
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                      <  2:602760/151‑1 (MQ=255)
 aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                     <  2:208023/151‑1 (MQ=255)
    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                  <  2:966561/151‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                 <  1:576390/149‑1 (MQ=255)
        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                    >  1:115193/1‑145 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGCTGTTGCGAAGGACCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGACGGgggggtattact                                                                                                                                             >  1:657684/1‑145 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg                                                                                                                                                  >  2:996060/1‑144 (MQ=255)
             cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCg                                                                                                                                          >  1:462825/1‑150 (MQ=255)
                     acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaa                                                                                                                                <  1:1304039/151‑1 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtggccccgtgttcc                                                                                                                                      >  2:403821/1‑132 (MQ=255)
                        tCGGTTTAAAAGCCAGTGAAAATGGCGTCCCCGTTTACCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATTCGCCGGGGCATGGCCGATATAAACGGCTAAGGCGAGGTTGCGGACGGGGTTGTGAACCTACGATCTCGCaaga                                                                                                                              <  2:1022784/150‑3 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGCCCTGCCATCCCGCgaaaa                                                                                                                             >  2:72533/1‑150 (MQ=255)
                           ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGTGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGCGGGGGGGGGGTTCTGCGATCCCGCcaaaa                                                                                                                             >  2:115726/1‑148 (MQ=255)
                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGTTACCGCCGTGCGGCAAAAAAGccc                                                                                                                        >  1:74218/1‑151 (MQ=255)
                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                         <  1:584906/150‑1 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                        <  2:922144/150‑1 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                                                                                       <  2:341982/150‑1 (MQ=255)
                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGTCTTCGAGCCCGGAAAAAAGccc                                                                                                                        >  1:702881/1‑147 (MQ=255)
                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGAGACCTCGATCCCGCCGAAAACGcc                                                                                                                         >  2:991375/1‑147 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCGGGAAAAAAGccc                                                                                                                        >  2:317670/1‑143 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGCGAAAACGCCCggggaggg                                                                                                                >  2:534008/1‑144 (MQ=255)
                                        tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTGTGCGGTCCCGCAAAAAAGCCCGCGAAGt                                                                                                                 >  2:827918/1‑147 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGGGGGGGGGGGAGCCTGTGATCCCGGCACAAAGGccggggag                                                                                                                   >  2:610432/1‑139 (MQ=255)
                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGGAAAAACGCCCCGGAAGGGATcg                                                                                                            >  1:241811/1‑150 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                       <  1:1095782/150‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCCGCAAAAAAGCCCGCGAAGGGGTCGCCGCCGTg                                                                                                    >  2:316724/1‑148 (MQ=255)
                                                       cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGCGATCCCGCGATCCCGCCAAAACGCCCCAGAAGTGATCGCCGGCGCGAAgg                                                                                                >  1:304018/1‑149 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGGCTGGGGTCAGGCAAAAAAGCCCGAGAAGGGATTGGCGCCGCGAAgg                                                                                                >  2:727704/1‑148 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAgg                                                                                                >  1:309825/1‑148 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAGCCGGCAAAAACGCCCGCGGAGTGATCGCCGCCGTGAAgg                                                                                                >  1:34195/1‑148 (MQ=255)
                                                         tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGGTGCTGCGATCCGGCAAAAAAGCCCCGGGAGGGATTGGCGCCGGTAAGGaaa                                                                                             >  2:74047/1‑148 (MQ=255)
                                                         tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAgg                                                                                                >  1:353109/1‑147 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCCAAAACGCCCGAGAAGAGATCCCCGCCGTGAAGGACAAAc                                                                                          >  1:910509/1‑150 (MQ=255)
                                                              cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCCAGAAGTGATAGCCGCCCTGAAGGACAAACAgg                                                                                       >  1:747633/1‑151 (MQ=255)
                                                              cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGTTCCGCAAAAACGCCCGAGAAGGGAGCGCCGCCGCGGAGGACAAACCgg                                                                                       >  2:580520/1‑151 (MQ=255)
                                                              cTGCGCGATGTTGCGAAGGCCCAGGTTGGCCCGGAGATGCCGCGGGGCATTGCCGACCTACCCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGCTCCGGCCATCCGGCAAAAACGCCCGATAATTGATCGCCGCCGTGAAGGACAAACTg                                                                                        <  2:231585/150‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                     <  1:1148555/151‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                     <  1:92079/151‑1 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:1348665/151‑1 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:550685/151‑1 (MQ=255)
                                                                  gcgATGTTGCGAAGGTCCAGATTGGCCTGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  2:198681/150‑1 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAAGACAAACTggga                                                                                     >  2:187342/1‑145 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                <  1:644946/151‑1 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGACCGGCCAAAAAGCCCCAGGAAGGATCGCCGCCGGGGAGGGCCAACAGGaacgcctgg                                                                              >  2:650537/1‑144 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGGGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCCTCGACGCGGTGCAAACCGCCGCGAAGCGAAACCCGGCGCGTAGGACAAACTggtgaatgcga                                                                              >  2:878190/1‑140 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGaa                                                                             >  1:850793/1‑150 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTTGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                            <  2:750087/150‑1 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGgaa                                                                            >  1:1169886/1‑150 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                            <  2:770576/150‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                         <  2:16302/151‑1 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGTTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                       <  1:875316/149‑1 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                     <  2:1116215/151‑1 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGGGAAGGGCAAAAAGGAAAAGCTGGAAagta                                                                        >  1:390410/1‑148 (MQ=255)
                                                                                     gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGTCCCGCAAAAACCCCCGAGAAGTGGTCGCCGCCGGGGAGGACAAAATGGGAACGCTGAAAAGGAGTCTGcc                                                                >  2:1165935/1‑150 (MQ=255)
                                                                                       ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                              <  1:1280703/150‑1 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                           <  1:911740/151‑1 (MQ=255)
                                                                                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  1:247146/151‑1 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAATCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                         <  1:837706/151‑1 (MQ=255)
                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGTAACCCGGGGTCCGGAAAAAACGCCCGAGGAGGGAGCGCCGCCGGGAAGGACAAAATTGGAAAGCTGGAAAATAGTCTGCCgggcgg                                                          >  1:730810/1‑144 (MQ=255)
                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGCCCCCGAAGAGATCGCCGCCCGGGAGGACAAAATGGAAAAGCTGAAAAGTAGTGTGCCCGGAAgcg                                                         >  2:1318355/1‑148 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  2:1083148/151‑1 (MQ=255)
                                                                                               gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTTAAAAGTAGTCTGCCGGAAGGCgggg                                                     >  1:450287/1‑151 (MQ=255)
                                                                                                agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                    <  2:832082/151‑1 (MQ=255)
                                                                                                agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGGGAAAAGGTCGCCGCCGTGGAGGACAAACTGGAAACGCCGAAAAGTAGTCTGCCGGGAggggggg                                                     >  2:1368303/1‑145 (MQ=255)
                                                                                                 gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                   <  1:1101239/151‑1 (MQ=255)
                                                                                                     cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                               <  1:70923/151‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                               <  2:747614/150‑1 (MQ=255)
                                                                                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGGCTGCCGGAAgggggggggaaaggt                                             >  1:336345/1‑138 (MQ=255)
                                                                                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                           <  2:647382/151‑1 (MQ=255)
                                                                                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                       <  1:934087/150‑1 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                    <  1:732031/150‑1 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                   <  2:944482/151‑1 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGGATCGGGAAAAAAGCCCCAGGAGTGATTGCCGCCGTGAAGGGCAAAATTGAAAAGGTGAAAAGGAGTTGGGCGGAAGGCGTGGAGATAATTACAACATa                                     >  1:1069350/1‑149 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                 >  1:360795/1‑151 (MQ=255)
                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:394997/151‑1 (MQ=255)
                                                                                                                     aaCTAAACGGTGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:532537/150‑1 (MQ=255)
                                                                                                                       cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGATACGGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGGAGGACAAACTGGGAACGCTGAAAAGTAGTCTGCCCGAAGGCGTGGAGATAGTTACAACATACGATc                                >  1:314060/1‑148 (MQ=255)
                                                                                                                        tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGATCCGGCAAAAACGCCCGAGGAGGGATCGGCGCCCGGAAGGGCAAACTGGAAACGCTGGAAAGGAGTCCGACGGAAGGGGGGGGGATAGGTAAAACATa                                     >  2:156483/1‑142 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                            <  2:787092/150‑1 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGATCCCGCAAAAACGCCCGGGAAGGGATCGCCCCCCTGTAGGACAAACTGGGAACGCTGAAAAGGAGGCTGCCGGCAGGGGGGGGGATAGTTATAACATAAGATCGccg                            >  1:285421/1‑148 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                    <  2:747526/151‑1 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAAATGGAAAAGCTGAAAAGAAGTCTGCCGGGAGGCGTGGGGGTAGTTTCAAAATACGATCGCAGCCAGCTcc                    >  2:929099/1‑150 (MQ=255)
                                                                                                                                   ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                 <  1:187348/151‑1 (MQ=255)
                                                                                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  1:1336335/151‑1 (MQ=255)
                                                                                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  1:474594/151‑1 (MQ=255)
                                                                                                                                             cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt       >  1:268176/1‑151 (MQ=255)
                                                                                                                                                 cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   >  2:159462/1‑151 (MQ=255)
                                                                                                                                                  ggggTGGTAATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAGCTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATGCGATCGCAGCCAGCTCATTGGCCGCGGCATAGACa  <  2:155337/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722

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