Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG44.1% 96.5 / 31.0 69V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (1/37);  new base G (30/0);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.73e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599426‑599721
                                                                                                                                                       |                                                                                                                                                  
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGggt                                                                                                                                                    <  1:117693/151‑1 (MQ=255)
    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGTTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGgggggggg                                                                                                                                                >  1:606917/1‑149 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGa                                                                                                                                               <  1:120552/149‑1 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGa                                                                                                                                               <  2:216971/149‑1 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGTGATcc                                                                                                                                            >  2:1088189/1‑150 (MQ=255)
          cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTg                                                                                                                                          <  1:521621/151‑1 (MQ=255)
              tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGgggggggg                                                                                                                                                >  2:783165/1‑139 (MQ=255)
              tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCt                                                                                                                                           >  2:181568/1‑146 (MQ=255)
                 aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑Tcc                                                                                                                                   <  2:75733/150‑1 (MQ=255)
                 aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCTGCGA‑Tc                                                                                                                                    >  1:408201/1‑149 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCttgggacccg                                                                                                                                  >  2:1169438/1‑142 (MQ=255)
                      caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGA‑TCCGGCaa                                                                                                                              >  2:1107352/1‑150 (MQ=255)
                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAAc                                                                                                                          <  2:223744/151‑1 (MQ=255)
                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGcc                                                                                                                       <  1:926373/150‑1 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGccc                                                                                                                      <  2:964091/150‑1 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGccc                                                                                                                      <  1:841246/150‑1 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGccc                                                                                                                      <  2:453667/150‑1 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCCGGGGGGGTG‑TCCTTCGG‑CCCGGCCAAAAAGccc                                                                                                                      >  2:818849/1‑149 (MQ=255)
                                aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCCGCGA‑TCCGGCCAAAAAGccc                                                                                                                      >  1:700161/1‑148 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGGTACCGGGA‑TCCAGCAAAAACCCCCg                                                                                                                     >  1:156402/1‑147 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGaa                                                                                                                 <  2:1198319/151‑1 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGG‑TCCGGCAAAAACGccc                                                                                                                      >  2:971244/1‑146 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGATCCTGCGG‑TCCGGCAAAAAAGCCCGAGAAGTg                                                                                                              >  2:905531/1‑147 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                          <  1:548119/151‑1 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                           <  2:134321/150‑1 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGGGCCCGCGA‑TCCGGCCAAAAACCCCgcggagg                                                                                                               >  1:846465/1‑138 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCAGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGTGATCCTGCGG‑ACCGGCAAAAACGCCCGCGAAGAGAGcggc                                                                                                        >  2:155294/1‑149 (MQ=255)
                                                  gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                 <  2:824064/151‑1 (MQ=255)
                                                  gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGA‑TCCGGCAAAAAAGCCCGGGAAGAGATCGCCGCCGTg                                                                                                  >  2:595754/1‑150 (MQ=255)
                                                     cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                              <  2:730708/151‑1 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGA‑TCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCGGGAAgg                                                                                              >  2:220338/1‑150 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATCGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                             <  2:965392/151‑1 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGGGGTCCTGCGC‑TCCGGCAAAAACGCCCGGGGAGGGGTCGCCGCCGTGaa                                                                                                >  2:448654/1‑146 (MQ=255)
                                                         tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCCGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGCGGAgg                                                                                              >  2:369244/1‑147 (MQ=255)
                                                          ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                          <  2:459457/150‑1 (MQ=255)
                                                           tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATTCTGCGG‑TCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGGGAAGGaaaa                                                                                          >  1:928145/1‑149 (MQ=255)
                                                             tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGA‑TCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAccg                                                                                      >  2:139101/1‑149 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGGCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAACt                                                                                       <  1:658680/147‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                     <  2:108539/149‑1 (MQ=255)
                                                                  gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                 <  1:1171673/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGGTCCTGCGA‑TCCGGCGAAAAAGCCCGGGAAGGGATCGCCGCCCCGAAAGACAAAATGGAAAAGcagggaa                                                                         >  1:907779/1‑143 (MQ=255)
                                                                                 tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGTGGTCCTGCGG‑TCCGGCAAAAAAGCCCGAGAAGGGATTGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAagta                                                                      >  2:178518/1‑147 (MQ=255)
                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATAGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                 <  2:663821/151‑1 (MQ=255)
                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGGTCCTGCGA‑TCCGGCAAAAAAGCCCGAGAAAGGATCCCCGCCGGGAAGGGAAAAACAGAAAAGCTGaaaaaaagttcg                                                                 >  1:851439/1‑141 (MQ=255)
                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                               <  1:1192422/151‑1 (MQ=255)
                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGTGATTCTGCGA‑AACGGCCAAAAAGCCCCGGAAGAGAGCGCCGCCGTGACGGACAAAACGGGAACGCGGAAAAGTAGTTTGCCggga                                                           >  1:650823/1‑144 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTCCTGCGA‑TCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGGAGGGCAAACTGGAAACGCTGAAAAAGAGGCTGCCGGGAGGCggg                                                     >  1:818414/1‑149 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGTACTGGGA‑TGCGGCAAAAAAGCCCGGGAAGTGATCGCCGCCGGGAAGGACAAAATGGAAAAGCTGAAAAGCAGTTTGCCGGGAg                                                          >  1:26142/1‑146 (MQ=255)
                                                                                               gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGTGATCCTGCCATTCCGGC‑AAAACGCCCGAGGAGGGATCGCCGGCGGGAAGGACAAACTGGGAAAGCCGAAAAGTGGTCCGCCGGGCGGCgggg                                                    >  1:124221/1‑151 (MQ=255)
                                                                                               gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGTGATCCTGCGA‑TCCGGCAAAAAAGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGGAAAGTAGTCTGCCGGGAGGCGTgg                                                    >  2:848566/1‑151 (MQ=255)
                                                                                                 gaTGCGCCGAGGCATTGCCGAACTAACCGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGT‑TCCGGCAAAAACACCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                    <  1:140812/149‑1 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                <  1:386725/150‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                             <  2:976476/151‑1 (MQ=255)
                                                                                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGG‑TCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGAAAAAACGGGAAAGCCGAAAAATAGTCTGCCGGGAgggggg                                                     >  1:935948/1‑138 (MQ=255)
                                                                                                                 gCCTAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATGGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                     <  2:660087/148‑1 (MQ=255)
                                                                                                                  ccGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGATCCTGCGA‑TCCGGCAAAAAAGCCCGGGGAGTGATCGGCCCCCGGAAGGACAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGGGGGGAGGTAGTTacaaacta                                    >  2:933294/1‑144 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                <  2:79265/151‑1 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                <  2:990523/151‑1 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                <  2:505680/151‑1 (MQ=255)
                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                               <  1:129687/151‑1 (MQ=255)
                                                                                                                     aaCTAAACGGCGAAGGCGAGGTGGCGGG‑CGGGGTGGTGATCCTGCGC‑TCCGGCAAAACCGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                               <  2:604265/150‑1 (MQ=255)
                                                                                                                       cTAAACGGCGAAGGCGAAGTGGCGGG‑GGGGGGGGGGGATCTGCGA‑TCCCGCACAAACGCCCGAGAAGGGGACGCCGCCGTGGAGGACAAAACGGAAAAGCTGGAAAAGAGTCTGCCGGGAGGGGGGGGGGGAGTTAAAAAATACGATCg                              >  1:823864/1‑149 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagc                          <  2:1092617/151‑1 (MQ=255)
                                                                                                                              gcgaaggcgaagTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCt                     <  2:112699/151‑1 (MQ=255)
                                                                                                                              gcgaaggcgaagTGGCGGG‑CGGGGGGGTGATTCTGCGA‑TCCCGCAAAAAACCCCCAGAAGTGATCGCCGCCGAGAAAGAAAAACTGGAAACGCTTAAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAACATACGATCGcagc                          >  2:1190049/1‑146 (MQ=255)
                                                                                                                               cgaaggcgaagTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                    <  2:405421/151‑1 (MQ=255)
                                                                                                                                  aggcgaagTGGCGGG‑CGGGGGGGGGATCCTGCGC‑TCCCGCGAAAAAGCCCGAGGAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGGTTGGCGGGAGGCGTGGAGATAGTTATAAAATACGATCGCAGGCAGCCCAtt                 >  1:50122/1‑151 (MQ=255)
                                                                                                                                      gaagTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc             <  2:287865/151‑1 (MQ=255)
                                                                                                                                       aagTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc             <  2:1184994/150‑1 (MQ=255)
                                                                                                                                        agTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc           <  1:617645/151‑1 (MQ=255)
                                                                                                                                              ggg‑cgggGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc     <  2:103951/151‑1 (MQ=255)
                                                                                                                                                   ggggAGGTGAACCTGCGA‑TCCGGCAAAAACGCCCGAGAAGCGATCGCCGCCGCGAAGGACCAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGCGGTGATAGCTACAACATACGATCGCAGCCAGCTCATTGACCGCGTCATCggc  >  1:456661/1‑148 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGTGGTGATCCTGCGA‑TCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599426‑599721

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