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

GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  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)
                 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:550685/151‑1 (MQ=255)
                                                             cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                    <  1:1348665/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)
                                                                                                                                                     |                                                                                                                                                 
GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  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: