Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG49.1% 64.6 / 24.9 59V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (5/25);  new base G (29/0);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725
                                                                                                                                                      |                                                                                                                                                     
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                       <  1:937170/151‑1 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                     >  1:1029372/1‑150 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                                                                                    <  2:941468/151‑1 (MQ=255)
      tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg                                                                                                                                                   >  1:716584/1‑147 (MQ=255)
          cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggggac                                                                                                                                                >  2:1153129/1‑143 (MQ=255)
           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                            <  1:358099/151‑1 (MQ=255)
              tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGTGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGc                                                                                                                                            >  2:851570/1‑148 (MQ=255)
               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGa                                                                                                                                          <  1:1120225/149‑1 (MQ=255)
                     acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCggtgc                                                                                                                                       >  2:195808/1‑142 (MQ=255)
                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCCCGCAAAAACAccc                                                                                                                          >  2:158264/1‑151 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                          <  2:1072238/150‑1 (MQ=255)
                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga                                                                                                                      <  2:442455/151‑1 (MQ=255)
                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga                                                                                                                      <  1:998479/151‑1 (MQ=255)
                                   gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                       >  1:381696/1‑148 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGTCTGCGGTCCCGCAAAAACAcc                                                                                                                           >  2:1036133/1‑143 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGACCCGGAAAAAAGccccgg                                                                                                                       >  1:849878/1‑143 (MQ=255)
                                       gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGTg                                                                                                                  >  1:5332/1‑149 (MQ=255)
                                        tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGATCCTGCCATCCGGCCAAAACGcc                                                                                                                           >  2:18533/1‑139 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATCCGGCGAAAAAGCCCGCGGAGGGGTcgc                                                                                                             >  2:449968/1‑150 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTACCGCGATCCGGCAAAAACGCCCGCGAAGGGATcg                                                                                                              >  2:308426/1‑149 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                        <  2:149780/151‑1 (MQ=255)
                                                gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt                                                                                                       <  1:43173/151‑1 (MQ=255)
                                                  gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                     <  1:36151/151‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                    <  2:114029/150‑1 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCGCGGAGTGGTCGCCGCCGCGGAgg                                                                                                  >  2:1035398/1‑150 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                               <  1:120866/151‑1 (MQ=255)
                                                           tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGATCTGCGATCCGGCAAAAACGCCCGAGAAGAGATCGCCCCCCTGAAGGACaaa                                                                                             >  1:1146351/1‑150 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAc                                                                                            >  1:1149430/1‑150 (MQ=255)
                                                             tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCCGTCCGGCAAAAACGCCCGGGGAGCGAGCGCCGCCGTGAAGGACAAAccg                                                                                          >  2:796873/1‑149 (MQ=255)
                                                              cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGGAGTGGTCGCCGCCGTGGAGGGCAAACTgg                                                                                         >  2:474832/1‑151 (MQ=255)
                                                               tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTgg                                                                                         >  1:1141488/1‑150 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGGTCCAGCGAGCCGGGAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAACTGGAAACGc                                                                                   >  1:425762/1‑150 (MQ=255)
                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGc                                                                                   >  1:845028/1‑149 (MQ=255)
                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGAGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 >  1:785525/1‑151 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                               <  2:518668/150‑1 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTCCGATTCGGCAAAAACGCCCCAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGgaa                                                                              >  2:1132585/1‑150 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:728599/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:231638/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                              >  1:579539/1‑148 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                       <  2:569531/151‑1 (MQ=255)
                                                                                 tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCCCCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGCGAAGGAACAACTGGAAAAGCTGGAAAATAATCt                                                                      >  1:230536/1‑151 (MQ=255)
                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      <  2:746852/150‑1 (MQ=255)
                                                                                   cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                    <  2:248234/151‑1 (MQ=255)
                                                                                   cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAACGCTGAAAAGAAGTCTg                                                                     >  2:252383/1‑150 (MQ=255)
                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGcc                                                                   >  1:545197/1‑151 (MQ=255)
                                                                                    aGATTGGCCCGCAGATGCCCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGCACGGGGGGGGATCTCCCGTTTCGGGCAAAAAGCGCAGGAAGTGATCGACGGCGTGACGGACCAACTGGAAAAGCTCaaaaaaagaattg                                                                    >  2:697882/1‑141 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                         <  2:1063012/151‑1 (MQ=255)
                                                                                                  aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                     <  1:581991/151‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGATCCTGCGATCCGGCAAAAAAAACCGAGAAGAGAACGCCGCCGTGAAGGAGCAACTGGGAACCCTGAAAAGGAGTCTGCCGGAAGGCGGGGAGATAg                                                  >  1:539234/1‑150 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGGAAAGGAGTCTGCCGGAAGGCGGGGGGATa                                                   >  2:950354/1‑149 (MQ=255)
                                                                                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGGCAAAACGCCCCAGAAGGGAACGCCGCCGGGCAAGACAAACTGGAAACGCTGGAAAGTAATCTGCCGGGAGGCGGGGGGATAGtt                                                >  1:681091/1‑151 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                    <  2:19025/151‑1 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                    <  1:574779/151‑1 (MQ=255)
                                                                                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTTCGATCCGGCAAAAAAGCCCGAGAAGGGAGCGGCGCCGTGAAGGGCAAACTGGGAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTACAACATAAGAACg                                  >  2:169167/1‑151 (MQ=255)
                                                                                                                      aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                                 <  2:1029201/151‑1 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGAGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagc                              <  1:605734/151‑1 (MQ=255)
                                                                                                                          aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGctg                               <  2:431086/149‑3 (MQ=255)
                                                                                                                          aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACTGGGAACGCTGGAAAGAAGTCTGCCGGAAGGCGTGGAGATAGGTACAACATACGATCGcagcc                             >  1:917752/1‑151 (MQ=255)
                                                                                                                                  aggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGGTAGTTACAACATACGATCGCAGCCAGCTCAtt                     >  1:565796/1‑151 (MQ=255)
                                                                                                                                  aggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGGAAAAGAGTCTGCCGGGAGGCGGGGAGGTAGTTACAACATAAGATCGCAGCCAGCTCAtt                     >  2:493664/1‑151 (MQ=255)
                                                                                                                                   ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                    <  2:853816/151‑1 (MQ=255)
                                                                                                                                     cgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                  <  1:691471/151‑1 (MQ=255)
                                                                                                                                       aagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                 >  2:371667/1‑150 (MQ=255)
                                                                                                                                          tGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             >  1:136768/1‑151 (MQ=255)
                                                                                                                                                     gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  2:906254/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725

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