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

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599426‑599726
                                                                                                                                                      |                                                                                                                                                        
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                          <  2:223111/151‑1 (MQ=255)
 taCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGACCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGggtg                                                                                                                                                         >  1:972026/2‑151 (MQ=255)
 aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                         <  1:632192/151‑1 (MQ=255)
  aCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg                                                                                                                                                        <  1:1180465/151‑1 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                                                                                       >  1:667640/1‑150 (MQ=255)
      tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggggt                                                                                                                                                    >  1:1016023/1‑147 (MQ=255)
        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggggggcc                                                                                                                                                  >  2:762390/1‑145 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATcc                                                                                                                                                  >  1:68247/1‑148 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                               <  1:388864/151‑1 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                               <  2:546422/151‑1 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTg                                                                                                                                                >  1:685583/1‑150 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCgg                                                                                                                                       <  1:680094/151‑1 (MQ=255)
                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCa                                                                                                                                     <  2:1354956/151‑1 (MQ=255)
                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGCT‑cc                                                                                                                                         >  2:1292395/1‑147 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAG‑CCGGCgaaaa                                                                                                                                 >  2:568749/1‑150 (MQ=255)
                            tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGc                                                                                                                              <  1:847999/150‑1 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCCCGGCaaaa                                                                                                                                  >  2:488952/1‑145 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGGAAAAAAGccc                                                                                                                            >  1:934465/1‑149 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCCAGAT‑CCCGCaaaa                                                                                                                                  >  2:692731/1‑143 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCgcg                                                                                                                         >  2:779405/1‑147 (MQ=255)
                                   gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAg                                                                                                                      <  2:321278/151‑1 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACCGGCAT‑CCCGCAAAAAAGccc                                                                                                                            >  2:971581/1‑144 (MQ=255)
                                     cAGTGAAAATGGCGTGCGCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCGCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGt                                                                                                                     <  2:206274/150‑1 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGG‑CCCGGCAAAACAACCGCGAAGTg                                                                                                                    >  2:760912/1‑151 (MQ=255)
                                       gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                  <  1:896591/151‑1 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  2:66047/150‑1 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  1:229207/150‑1 (MQ=255)
                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAACGCCCGGGAAGGGAAcg                                                                                                                >  1:1096413/1‑150 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  1:829495/150‑1 (MQ=255)
                                                gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGT‑CCGGCCAAAAAGCCCGAGGAGGGATCgccgcc                                                                                                           >  2:1069805/1‑149 (MQ=255)
                                                 cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGATCgccgcc                                                                                                           >  2:1466685/1‑148 (MQ=255)
                                                     cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GTGAAgg                                                                                                   >  2:716912/1‑151 (MQ=255)
                                                       cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGAc                                                                                                 <  2:93262/151‑1 (MQ=255)
                                                          ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACaa                                                                                               <  1:158984/150‑1 (MQ=255)
                                                          ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACaa                                                                                               <  1:1423670/150‑1 (MQ=255)
                                                          ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCCCC‑GGGAAGGGCaa                                                                                               >  2:1136918/1‑150 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GGGAAgggaaaa                                                                                              >  2:667084/1‑144 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGACTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATAGCCGCC‑GTGAAGGACaaa                                                                                              >  2:202110/1‑149 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGGGAT‑CCCGCAAAAACGCCCGAGAAGGGAGCGCCGCC‑CTGAAGGGCaaa                                                                                              >  2:29173/1‑149 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGTCGCC‑GTGAAGGACAAACTgg                                                                                          <  1:1273265/149‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGaa                                                                                        <  2:264008/151‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGa                                                                                         <  1:1180106/150‑1 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GTGAAGGACAAACTGGaa                                                                                        >  1:1311289/1‑150 (MQ=255)
                                                                  gcgATGTTGCGAAGGTCAAGATTGGCCGGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCAGCC‑GTGAAGGACAAACTGGaa                                                                                        <  1:365961/149‑1 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAAc                                                                                      <  1:541855/149‑1 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACg                                                                                     <  2:703712/150‑1 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAG‑CCGGCAAAAACGCCCGGGAAGGGATCGCCGCC‑CGGCAGGGCAAAATGGaaa                                                                                       >  2:263064/1‑148 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCt                                                                                   <  1:1295419/151‑1 (MQ=255)
                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGAGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGGGACCGGCCT‑CCGGGCAAAACGACCACGAAGGGATGGCCCCC‑GCGAAGGAc                                                                                                 >  1:60118/1‑136 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTg                                                                                  <  2:354677/150‑1 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGaa                                                                                <  2:847653/150‑1 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCCAAAAAGCCCGAGAAGGGATCGCCGCC‑GGGAAGGACAAACTGGAAACGCCGaaaa                                                                              >  2:121105/1‑150 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAACCCCGGGAAGGGAGCGCCGCC‑GGGAAGGGCAAACTGGAAACGCCGaaaa                                                                              >  1:1405350/1‑150 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGAA‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCC‑GCGGAGGGCAAACTGGAAACCCTGaaaa                                                                              >  2:1013379/1‑148 (MQ=255)
                                                                              aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                          <  2:751080/151‑1 (MQ=255)
                                                                                     gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GGGAAGGACAAACTGGAAAAGCTGAAAAGAAGTCTGcc                                                                    >  1:997488/1‑150 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                  >  1:1323727/1‑148 (MQ=255)
                                                                                            cgggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                              <  2:139984/147‑1 (MQ=255)
                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                            <  1:1264067/150‑1 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTTCGGT‑CCCGCAAAAACGCCCCGGAAGGGGTCGGCGCCGGTGAAGG‑CAAAATGGGAAAGCTGAAAAAGATTCTGGCGGGAgg                                                              >  2:402782/1‑147 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCGAGAAGGGATCGCCGCC‑CGGGAGGAAAAACAGGAAAAGCTGAAAAAAAGGCTGCCGGAAgggggg                                                          >  2:905282/1‑147 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                     <  1:1015360/150‑1 (MQ=255)
                                                                                                     cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                                   <  1:762468/151‑1 (MQ=255)
                                                                                                     cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAACGCCCCAGAAGGGATCGCCGCC‑CTGAAGGACAAAATGGAAACGGTGAAAAGGAGTCTGCCGGGAGGCgggggg                                                       >  2:365755/1‑145 (MQ=255)
                                                                                                     cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCCGGGAT‑CCGGGAAAAACGCCCGAGAAGGGGGAGGCGCC‑GGGAAGGAACAAATGGGAACGCTGAAAAAGAGTCTGCCGGGAGGCgggggg                                                       >  2:926878/1‑145 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                                   <  2:1015927/150‑1 (MQ=255)
                                                                                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGAT‑TCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGAAAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGAGATAGtt                                                 >  1:1057271/1‑149 (MQ=255)
                                                                                                          gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGAT‑CTGGCAAAAACGCCCGCGAAGTGAGCGCCCCC‑GGGAAGGACAAACCGGAAACGCCGGAAAGTAGTCTGCCGGCAGGCGTGGAGATAGtt                                                 >  1:156586/1‑148 (MQ=255)
                                                                                                           ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                             <  2:760013/151‑1 (MQ=255)
                                                                                                           ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTaca                                              >  2:926369/1‑150 (MQ=255)
                                                                                                             cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGAGCGCCGCC‑GGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTacaa                                             >  2:1286026/1‑149 (MQ=255)
                                                                                                               ttGCCGAACTAAGCGGCGAAGGCGAAGGGGCGGGCGGGGTGGGTACTCTGCTAT‑CCGGAAAAAACGCCCGAGAAGGGATCGCAGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                              <  2:299020/146‑1 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGGTCCTCCGAT‑CCCGCACAAACACCCGGGAAGCGATTGGCGCC‑CCGAAGGAGAAACCAGAAACGCCGGAAAATAATCTGCCGGAAGGCGCGGAGATaagtttcaaaaa                                         >  1:125224/1‑139 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                         <  2:561921/150‑1 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                         <  2:667568/150‑1 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑CCGGCAAAAACGCCCGAGAAGAGAGCGCCGCC‑GTGAAGGACAAAATGGGAAAGCTGGAAAGAAGTCTGCCGGGAGGGGGGGGGAGAGTTAAAACATa                                         >  2:1065438/1‑150 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAct                                       <  1:619790/151‑2 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                        <  1:251385/150‑1 (MQ=255)
                                                                                                                  ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                       <  1:433311/150‑1 (MQ=255)
                                                                                                                  ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCCGCGAT‑CCCGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GTGGAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGTGGGGATAGTTTCAACATACg                                       >  1:579494/1‑150 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                     <  1:1456538/151‑1 (MQ=255)
                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGGCAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGGTAGTTACAACATACGATc                                    >  1:1392588/1‑151 (MQ=255)
                                                                                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                   <  2:1405242/151‑1 (MQ=255)
                                                                                                                       cTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGCCTGCGAA‑TCGGCAAAAACGCCCGGGAAGGGATCGCCGCC‑GGGAAAGACAAAATGGAAAAGCTGAAAAGTAGACTGCCGGGAGGGGGGGAGATAGTTacaacaaaaga                                      >  2:601965/1‑141 (MQ=255)
                                                                                                                           aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcc                              >  1:1043688/1‑150 (MQ=255)
                                                                                                                             ggcgaaggcgaaTTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccagc                           <  2:971931/151‑1 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                         <  1:476964/150‑1 (MQ=255)
                                                                                                                                   ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                     <  1:475410/151‑1 (MQ=255)
                                                                                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                  <  2:934376/151‑1 (MQ=255)
                                                                                                                                         gTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg               <  2:68326/151‑1 (MQ=255)
                                                                                                                                             cgggcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt           >  2:1187402/1‑151 (MQ=255)
                                                                                                                                                gcgggGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGa        <  1:1105632/151‑1 (MQ=255)
                                                                                                                                                   gggTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa     <  1:1146702/151‑1 (MQ=255)
                                                                                                                                                   gggTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa     <  2:685507/151‑1 (MQ=255)
                                                                                                                                                      tggtgATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  1:1260156/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCC‑GTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599426‑599726

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: ‑