Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG45.2% 89.1 / 25.9 63V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (4/30);  new base G (28/0);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599429‑599722
                                                                                                                                                    |                                                                                                                                                  
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Gggggggg                                                                                                                                                 >  1:729755/1‑150 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggggg                                                                                                                                                >  1:676399/1‑149 (MQ=255)
    cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggg‑ggggggggg                                                                                                                                                >  2:430515/1‑146 (MQ=255)
    cgacgaCTTTAATCACATCGTTTTAAAAATCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGGTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGCACTAAACGGCGAAGGCGAAGTGGCAggg‑ggggggggg                                                                                                                                                >  2:929284/1‑146 (MQ=255)
          cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTg                                                                                                                                          >  2:105855/1‑148 (MQ=255)
          cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCTGc                                                                                                                                         >  2:687953/1‑149 (MQ=255)
           tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTcc                                                                                                                                            >  2:36710/1‑145 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCCGCGATc                                                                                                                                     >  1:1103196/1‑151 (MQ=255)
                      cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCaaaaa                                                                                                                            <  1:899085/150‑1 (MQ=255)
                       gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                          <  1:438934/151‑1 (MQ=255)
                       gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGGTCCCGCaaaaa                                                                                                                            >  2:304789/1‑149 (MQ=255)
                        ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCaaaaa                                                                                                                            >  1:316514/1‑148 (MQ=255)
                          ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                        <  1:642608/150‑1 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGccc                                                                                                                       >  1:937943/1‑146 (MQ=255)
                                gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                 <  1:1057427/151‑1 (MQ=255)
                                  cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCCTCCGA‑CCCGCAAAAAACCCCgag                                                                                                                    >  2:1056468/1‑146 (MQ=255)
                                        aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑AGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGGGAAGGGGTcgcc                                                                                                         >  1:1053597/1‑151 (MQ=255)
                                          aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                       <  1:639969/151‑1 (MQ=255)
                                          aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGCGAAGGGATcgccgc                                                                                                       >  1:537930/1‑151 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                      <  1:151228/150‑1 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                     <  2:771491/151‑1 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTTCGGTCCGGCAAAAAAGCCCGAGAAATGATCgccgcc                                                                                                      >  1:605974/1‑150 (MQ=255)
                                                           cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGGAAAAAAGCCCGAGGAGGGAACGCCGCCCTGAAGGACaaaaagg                                                                                      >  2:44681/1‑147 (MQ=255)
                                                           cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                      >  1:1110015/1‑151 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                        <  1:983470/147‑1 (MQ=255)
                                                              cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGGCCCGCGGACCGGGCAAAAAGCCCCAGGAGGGGTCGCCGCCGGGGAGGACaaaacggga                                                                                    >  1:63489/1‑144 (MQ=255)
                                                              cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                   <  1:976929/151‑1 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGaaa                                                                                   >  1:253467/1‑150 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGCGGGCCTGCGATACGGCCAAAACCCCCGAAAAGAGATGTCCGCCGCGaaagcaaaaccagaa                                                                                    >  2:1023871/1‑136 (MQ=255)
                                                                  aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                               <  2:191925/151‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                             <  1:466753/151‑1 (MQ=255)
                                                                     ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                             >  1:99101/1‑150 (MQ=255)
                                                                          aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                       <  1:647750/151‑1 (MQ=255)
                                                                            ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCCGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGAAGGACAAACTGGAAACGCTGGAAagaa                                                                       >  2:115173/1‑147 (MQ=255)
                                                                            ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGTGGTCGCCGCCGGGAAGGACAAACTGGAAAAGCTGcaaa                                                                          >  2:211348/1‑146 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                   <  1:646862/150‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGGCAAACTGGAAACGCTGCAAAGTAGTc                                                                    >  2:360458/1‑149 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                 <  1:143259/151‑1 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCGCGATCCCGCAAAAAAGCCCGAGAAGGGATTGCCGCCGGGAAGGACAAACTGGGAACGCTGAAAAGTAGGCTg                                                                  >  2:582814/1‑150 (MQ=255)
                                                                                  gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                               <  2:1041326/151‑1 (MQ=255)
                                                                                   tttagttcGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                <  1:881202/110‑1 (MQ=255)
                                                                                     tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGAAAAGTAGTCTGCCgg                                                              >  1:235618/1‑149 (MQ=255)
                                                                                       gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  2:274585/151‑1 (MQ=255)
                                                                                       gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGAAgg                                                          >  1:353814/1‑151 (MQ=255)
                                                                                       gCCCGGAGATGCGCCGGAGCATTGCCGAGCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACCCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  2:162452/151‑1 (MQ=255)
                                                                                         ccGGAGATGCGCCGGGGCATTGCCGAACTAAAAGGCGAAGGCGAAGTGGCGGCC‑GGGGTGGCGAACCTTCCATCCGGCCAAAACGCCCGCGAAGGGCATGCCCCCCGGCAGGCCAAAATGGTGAAGCTGAAAaggagggcgcaccg                                                             >  2:723085/1‑137 (MQ=255)
                                                                                          cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                       <  1:728221/151‑1 (MQ=255)
                                                                                             agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                    <  1:38922/151‑1 (MQ=255)
                                                                                               aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTTCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGAAACGCTGGAAAGTAGTCTGCCGGGAGGCgggg                                                     >  2:816167/1‑148 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                              <  2:597788/151‑1 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                              <  2:414632/151‑1 (MQ=255)
                                                                                                    ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                             <  2:939901/151‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                    <  1:990575/150‑1 (MQ=255)
                                                                                                               ccGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                  <  2:178160/151‑1 (MQ=255)
                                                                                                               ccGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGGTACAACATa                                     >  1:1106621/1‑148 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:382529/151‑1 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  1:402990/151‑1 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATTCTGCGATACCGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGGAAGCGGGGGGATAGTTACAACATACGATc                                >  2:281916/1‑151 (MQ=255)
                                                                                                                       aaCGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGGCTGCCGGAAAGCGTGGAGATAGGTACAAAATACGAACGcagcc                          >  1:439907/1‑151 (MQ=255)
                                                                                                                        aCGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGcagcc                          >  1:550638/1‑150 (MQ=255)
                                                                                                                         cGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCCCGTAGGGCCAACCGGAAACCCCGCAAAGGAGTCTGGCGGCAGGCGGGGAGATAGGTGCAACATACGAACGcaggc                          >  2:293908/1‑147 (MQ=255)
                                                                                                                                     agTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc            <  1:168382/151‑1 (MQ=255)
                                                                                                                                        ggcgggc‑gGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc         <  2:934344/151‑1 (MQ=255)
                                                                                                                                           gggc‑gggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc      >  1:1028494/1‑151 (MQ=255)
                                                                                                                                             gc‑gggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg     >  1:264726/1‑150 (MQ=255)
                                                                                                                                              c‑gggGAGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   >  2:271117/1‑151 (MQ=255)
                                                                                                                                                ggggTGGTGAGCCTTCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  <  2:1081987/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599429‑599722

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