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

GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599430‑599722
                                                                                                                                                     |                                                                                                                                                 
gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                  <  2:966561/151‑1 (MQ=255)
   cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                 <  1:576390/149‑1 (MQ=255)
     acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGCTGTTGCGAAGGACCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGACGGgggggtattact                                                                                                                                             >  1:657684/1‑145 (MQ=255)
         cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCg                                                                                                                                          >  1:462825/1‑150 (MQ=255)
                 acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCaa                                                                                                                                <  1:1304039/151‑1 (MQ=255)
                   aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtggccccgtgttcc                                                                                                                                      >  2:403821/1‑132 (MQ=255)
                    tCGGTTTAAAAGCCAGTGAAAATGGCGTCCCCGTTTACCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATTCGCCGGGGCATGGCCGATATAAACGGCTAAGGCGAGGTTGCGGACGGGGTTGTGAACCTACGAT‑CTCGCaaga                                                                                                                              <  2:1022784/150‑3 (MQ=255)
                     cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGCCCTGCCAT‑CCCGCgaaaa                                                                                                                             >  2:72533/1‑150 (MQ=255)
                       ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGTGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGCGGGGGGGGGGTTCTGCGAT‑CCCGCcaaaa                                                                                                                             >  2:115726/1‑148 (MQ=255)
                         ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGTTACCGCCGT‑GCGGCAAAAAAGccc                                                                                                                        >  1:74218/1‑151 (MQ=255)
                         ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGcc                                                                                                                         <  1:584906/150‑1 (MQ=255)
                          tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGccc                                                                                                                        <  2:922144/150‑1 (MQ=255)
                           aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCg                                                                                                                       <  2:341982/150‑1 (MQ=255)
                             aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGTCTTCGAG‑CCCGGAAAAAAGccc                                                                                                                        >  1:702881/1‑147 (MQ=255)
                             aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGAGACCTCGATCCCGCCGAAAACGcc                                                                                                                         >  2:991375/1‑147 (MQ=255)
                                 cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGT‑CCGGGAAAAAAGccc                                                                                                                        >  2:317670/1‑143 (MQ=255)
                                 cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGAT‑CCGGCGAAAACGCCCggggaggg                                                                                                                >  2:534008/1‑144 (MQ=255)
                                    tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTGTGCGGT‑CCCGCAAAAAAGCCCGCGAAGt                                                                                                                 >  2:827918/1‑147 (MQ=255)
                                     gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGGGGGGGGGGGAGCCTGTGATCCCGGCACAAAGGccggggag                                                                                                                   >  2:610432/1‑139 (MQ=255)
                                      aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGGAAAAACGCCCCGGAAGGGATcg                                                                                                            >  1:241811/1‑150 (MQ=255)
                                           ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                       <  1:1095782/150‑1 (MQ=255)
                                                gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCCGCAAAAAAGCCCGCGAAGGGGTCGCCGCCGTg                                                                                                    >  2:316724/1‑148 (MQ=255)
                                                   cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGCGATCCCGCGAT‑CCCGCCAAAACGCCCCAGAAGTGATCGCCGGCGCGAAgg                                                                                                >  1:304018/1‑149 (MQ=255)
                                                    gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGGCTGGGGT‑CAGGCAAAAAAGCCCGAGAAGGGATTGGCGCCGCGAAgg                                                                                                >  2:727704/1‑148 (MQ=255)
                                                    gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAgg                                                                                                >  1:309825/1‑148 (MQ=255)
                                                    gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAG‑CCGGCAAAAACGCCCGCGGAGTGATCGCCGCCGTGAAgg                                                                                                >  1:34195/1‑148 (MQ=255)
                                                     tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGGTGCTGCGAT‑CCGGCAAAAAAGCCCCGGGAGGGATTGGCGCCGGTAAGGaaa                                                                                             >  2:74047/1‑148 (MQ=255)
                                                     tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAgg                                                                                                >  1:353109/1‑147 (MQ=255)
                                                        aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCCAAAACGCCCGAGAAGAGATCCCCGCCGTGAAGGACAAAc                                                                                          >  1:910509/1‑150 (MQ=255)
                                                          cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCCAGAAGTGATAGCCGCCCTGAAGGACAAACAgg                                                                                       >  1:747633/1‑151 (MQ=255)
                                                          cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGT‑TCCGCAAAAACGCCCGAGAAGGGAGCGCCGCCGCGGAGGACAAACCgg                                                                                       >  2:580520/1‑151 (MQ=255)
                                                          cTGCGCGATGTTGCGAAGGCCCAGGTTGGCCCGGAGATGCCGCGGGGCATTGCCGACCTACCCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGCTCCGGCCAT‑CCGGCAAAAACGCCCGATAATTGATCGCCGCCGTGAAGGACAAACTg                                                                                        <  2:231585/150‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                     <  1:1148555/151‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                     <  1:92079/151‑1 (MQ=255)
                                                             cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:550685/151‑1 (MQ=255)
                                                             cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:1348665/151‑1 (MQ=255)
                                                              gcgATGTTGCGAAGGTCCAGATTGGCCTGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  2:198681/150‑1 (MQ=255)
                                                                gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAAGACAAACTggga                                                                                     >  2:187342/1‑145 (MQ=255)
                                                                 aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                <  1:644946/151‑1 (MQ=255)
                                                                   gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGA‑CCGGCCAAAAAGCCCCAGGAAGGATCGCCGCCGGGGAGGGCCAACAGGaacgcctgg                                                                              >  2:650537/1‑144 (MQ=255)
                                                                     tGCGAAGGTCCAGATTGGCCCGGGGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCCTCGAC‑GCGGTGCAAACCGCCGCGAAGCGAAACCCGGCGCGTAGGACAAACTggtgaatgcga                                                                              >  2:878190/1‑140 (MQ=255)
                                                                     tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGaa                                                                             >  1:850793/1‑150 (MQ=255)
                                                                      gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTTGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                            <  2:750087/150‑1 (MQ=255)
                                                                      gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGgaa                                                                            >  1:1169886/1‑150 (MQ=255)
                                                                      gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                            <  2:770576/150‑1 (MQ=255)
                                                                        gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                         <  2:16302/151‑1 (MQ=255)
                                                                            gTCCAGATTGGCCCGGAGTTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                       <  1:875316/149‑1 (MQ=255)
                                                                            gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                     <  2:1116215/151‑1 (MQ=255)
                                                                            gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGT‑CCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGGGAAGGGCAAAAAGGAAAAGCTGGAAagta                                                                        >  1:390410/1‑148 (MQ=255)
                                                                                 gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGT‑CCCGCAAAAACCCCCGAGAAGTGGTCGCCGCCGGGGAGGACAAAATGGGAACGCTGAAAAGGAGTCTG‑cc                                                                >  2:1165935/1‑150 (MQ=255)
                                                                                   ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCgg                                                              <  1:1280703/150‑1 (MQ=255)
                                                                                     ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAg                                                           <  1:911740/151‑1 (MQ=255)
                                                                                      gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAgg                                                          <  1:247146/151‑1 (MQ=255)
                                                                                       cccGGAGATGCGCCGGGGCATTGCCGAACTAATCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGc                                                         <  1:837706/151‑1 (MQ=255)
                                                                                         cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGTAACCCGGGGT‑CCGGAAAAAACGCCCGAGGAGGGAGCGCCGCCGGGAAGGACAAAATTGGAAAGCTGGAAAATAGTCTG‑CCgggcgg                                                          >  1:730810/1‑144 (MQ=255)
                                                                                         cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCCGCAAAAAAGCCCCCGAAGAGATCGCCGCCCGGGAGGACAAAATGGAAAAGCTGAAAAGTAGTGTGCCCGGAAgcg                                                         >  2:1318355/1‑148 (MQ=255)
                                                                                          ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTg                                                      <  2:1083148/151‑1 (MQ=255)
                                                                                           gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTTAAAAGTAGTCTG‑CCGGAAGGCgggg                                                     >  1:450287/1‑151 (MQ=255)
                                                                                            agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGga                                                    <  2:832082/151‑1 (MQ=255)
                                                                                            agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCGGGAAAAGGTCGCCGCCGTGGAGGACAAACTGGAAACGCCGAAAAGTAGTCTG‑CCGGGAggggggg                                                     >  2:1368303/1‑145 (MQ=255)
                                                                                             gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGgag                                                   <  1:1101239/151‑1 (MQ=255)
                                                                                                 cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAg                                               <  1:70923/151‑1 (MQ=255)
                                                                                                  gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAg                                               <  2:747614/150‑1 (MQ=255)
                                                                                                   ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGGCTG‑CCGGAAgggggggggaaaggt                                             >  1:336345/1‑138 (MQ=255)
                                                                                                     ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTac                                           <  2:647382/151‑1 (MQ=255)
                                                                                                          aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTacaaca                                       <  1:934087/150‑1 (MQ=255)
                                                                                                             gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATAc                                    <  1:732031/150‑1 (MQ=255)
                                                                                                             gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACg                                   <  2:944482/151‑1 (MQ=255)
                                                                                                             gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGGA‑TCGGGAAAAAAGCCCCAGGAGTGATTGCCGCCGTGAAGGGCAAAATTGAAAAGGTGAAAAGGAGTTGG‑GCGGAAGGCGTGGAGATAATTACAACATa                                     >  1:1069350/1‑149 (MQ=255)
                                                                                                               cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGAt                                 >  1:360795/1‑151 (MQ=255)
                                                                                                                gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:394997/151‑1 (MQ=255)
                                                                                                                 aaCTAAACGGTGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:532537/150‑1 (MQ=255)
                                                                                                                   cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGAT‑ACGGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGGAGGACAAACTGGGAACGCTGAAAAGTAGTCTG‑CCCGAAGGCGTGGAGATAGTTACAACATACGATc                                >  1:314060/1‑148 (MQ=255)
                                                                                                                    tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑CCGGCAAAAACGCCCGAGGAGGGATCGGCGCCCGGAAGGGCAAACTGGAAACGCTGGAAAGGAGTCCG‑ACGGAAGGGGGGGGGATAGGTAAAACATa                                     >  2:156483/1‑142 (MQ=255)
                                                                                                                     aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                            <  2:787092/150‑1 (MQ=255)
                                                                                                                     aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGAT‑CCCGCAAAAACGCCCGGGAAGGGATCGCCCCCCTGTAGGACAAACTGGGAACGCTGAAAAGGAGGCTG‑CCGGCAGGGGGGGGGATAGTTATAACATAAGATCGccg                            >  1:285421/1‑148 (MQ=255)
                                                                                                                            gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                    <  2:747526/151‑1 (MQ=255)
                                                                                                                            gaaggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAAATGGAAAAGCTGAAAAGAAGTCTG‑CCGGGAGGCGTGGGGGTAGTTTCAAAATACGATCGCAGCCAGCTcc                    >  2:929099/1‑150 (MQ=255)
                                                                                                                               ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                 <  1:187348/151‑1 (MQ=255)
                                                                                                                                  gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  1:1336335/151‑1 (MQ=255)
                                                                                                                                  gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  1:474594/151‑1 (MQ=255)
                                                                                                                                         cgggcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt       >  1:268176/1‑151 (MQ=255)
                                                                                                                                             cgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   >  2:159462/1‑151 (MQ=255)
                                                                                                                                              ggggTGGTAATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAGCTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATGCGATCGCAGCCAGCTCATTGGCCGCGGCATAGACa  <  2:155337/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTG‑CCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599430‑599722

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