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

AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599427‑599721
                                                                                                                                                      |                                                                                                                                                   
aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                    <  1:386148/151‑1 (MQ=255)
   gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                 >  2:395718/1‑149 (MQ=255)
    cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGTTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg                                                                                                                                                  >  2:404791/1‑148 (MQ=255)
       gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggggcc                                                                                                                                             >  1:628325/1‑145 (MQ=255)
          gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtac                                                                                                                                             >  2:321696/1‑142 (MQ=255)
           acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg                                                                                                                                         <  1:724377/151‑1 (MQ=255)
           acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg                                                                                                                                         <  1:955680/151‑1 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggagatcgg                                                                                                                                         >  1:699491/1‑139 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                <  2:699468/142‑1 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCg                                                                                                                                         >  1:1876/1‑149 (MQ=255)
                     caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGGACCGCGATcc                                                                                                                                    >  2:1188618/1‑145 (MQ=255)
                      aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggttcgtcggt                                                                                                                                       >  2:1257891/1‑130 (MQ=255)
                      aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGTCCCGCaaa                                                                                                                              >  2:1051492/1‑150 (MQ=255)
                           tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGCTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGGTCCCGCGAAAAAGcc                                                                                                                        >  2:1125281/1‑151 (MQ=255)
                            ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTCCGGCCCGGCAAAAAc                                                                                                                           >  2:21466/1‑147 (MQ=255)
                                 aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCg                                                                                                                      >  1:397589/1‑147 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                <  1:1091354/151‑1 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                <  2:760721/151‑1 (MQ=255)
                                    cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGCGGTCCTGCGGACTCGCCAAAAAGCCCGAGGAGTg                                                                                                               >  1:494166/1‑151 (MQ=255)
                                     aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAggg                                                                                                               >  1:496351/1‑148 (MQ=255)
                                      gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                             <  2:222616/151‑1 (MQ=255)
                                        gAAAATGGCGTGCCCGGTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                            <  2:106039/150‑1 (MQ=255)
                                          aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGATCCGGCAAAAACGCCCGGGAAGTGATcgc                                                                                                          >  1:620731/1‑150 (MQ=255)
                                            aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGCCCCCGGGCCCGCAAAAAAGACCGCGaaagggtcccc                                                                                                         >  2:1234014/1‑140 (MQ=255)
                                             tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTCCGGTCCCGGAAAAACGCCCGGGGAGGGATTgccgcc                                                                                                      >  2:470396/1‑151 (MQ=255)
                                             tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGATCCGCGGTCCGGCAAAAACGCCCGAGAAGGGAGCgccgcc                                                                                                      >  2:865804/1‑151 (MQ=255)
                                                 gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTGCCGCAAAAACACCCGAGAAGAGGTcgccg                                                                                                        >  2:688396/1‑145 (MQ=255)
                                                   gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                 <  2:416664/150‑1 (MQ=255)
                                                    cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                 <  1:423071/149‑1 (MQ=255)
                                                    cccGTTTATCGGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                               <  2:772753/151‑1 (MQ=255)
                                                      cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCTTCGCAGCGGGCAAAAACCCCGGGAAATGCCCCTCGCCGGGAAgg                                                                                               >  2:895933/1‑149 (MQ=255)
                                                          tATCTGCGCGATGTCGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCCTCGTGCCGGCAAAAAAGGCCCCGGAGGGATCgccgccccggagtgcaaa                                                                                          >  2:81725/1‑138 (MQ=255)
                                                           aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                         <  2:804859/150‑1 (MQ=255)
                                                                cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCCGCAAAAACGCCCGAGAAGTGATTGCCGCCGTGGAGGGCAAAATgg                                                                                     >  1:601198/1‑148 (MQ=255)
                                                                  cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCCGCGATCCGGCCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                 >  1:1177137/1‑150 (MQ=255)
                                                                   gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCACTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGGTCGCCGCCGGGAAGGACCAAATGGAAAAGc                                                                               >  1:2483/1‑151 (MQ=255)
                                                                     tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCGGGCAAAACGCCCGAGAAGTGAGTGCCGCCGGGAAGGACAAAATGGgaaacctg                                                                             >  2:440457/1‑146 (MQ=255)
                                                                        tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                           >  1:919200/1‑150 (MQ=255)
                                                                          cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGTCTGCGCACCGGCAAAAACCGCCGCGGAGTGATCGGCCACGGGAAGGGCAAACTCGAAAAGCCGcaaa                                                                         >  1:687562/1‑150 (MQ=255)
                                                                            aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGACCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGAGAAACTGGaaa                                                                                  >  2:428353/1‑139 (MQ=255)
                                                                             aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                     <  2:1251972/151‑1 (MQ=255)
                                                                               gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                     >  1:160052/1‑149 (MQ=255)
                                                                                 ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGCGGACCCGCCAAAACGCCCGAGAAAGGATCGCCGCCGTGGAGGGCCAACATGGAAAcccgaaaaataatc                                                                    >  2:385991/1‑136 (MQ=255)
                                                                                  cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCGGc                                                                <  2:723300/151‑1 (MQ=255)
                                                                                      ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                             >  1:172987/1‑150 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                            <  2:566758/150‑1 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAACGCCCGAGAAGGGATCGGCGCCGGGGAGGACAAAATGGAAAAGCTGGAAAGGAGTCTGCCgg                                                             >  1:551999/1‑149 (MQ=255)
                                                                                                gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                  <  2:1038947/151‑1 (MQ=255)
                                                                                                  tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                 <  1:948819/150‑1 (MQ=255)
                                                                                                        ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCCCGATCCCGCCAAAAACCCCCAGAAGGGGTCGCCCCCCTGGAGGGCAAACTGGAAAAGCTGAAAAGTAGTCTGCCCGGAGGCGTGGAGATAGtt                                            >  2:250471/1‑149 (MQ=255)
                                                                                                           gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                       <  1:1209215/151‑1 (MQ=255)
                                                                                                           gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                       <  2:208740/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                     <  2:142309/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                     <  1:331456/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGATCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCCTGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTacaaca                                      >  1:970545/1‑150 (MQ=255)
                                                                                                                gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                  <  2:445364/151‑1 (MQ=255)
                                                                                                                gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGGTAGTTACAACATa                                    >  1:129954/1‑149 (MQ=255)
                                                                                                                   gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                <  1:1028587/150‑1 (MQ=255)
                                                                                                                   gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                               <  2:778444/151‑1 (MQ=255)
                                                                                                                    aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATc                               >  1:33991/1‑150 (MQ=255)
                                                                                                                         aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATAGTTACAACATACCATCGcagcc                         >  1:924395/1‑151 (MQ=255)
                                                                                                                          aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGTGATTGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGGTAGTTACAAAAAAAGATCGcagcc                         >  2:1243017/1‑150 (MQ=255)
                                                                                                                           cGGCGAAGGTGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCCGCAAAAAAGCCCGCGAAGGGATCGCCGCCGGGGAGGACAAACCGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGAAAGTTACAACATACGATCGcagcc                         >  2:575685/1‑149 (MQ=255)
                                                                                                                           cGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                       <  1:381187/151‑1 (MQ=255)
                                                                                                                                aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                  <  2:405200/151‑1 (MQ=255)
                                                                                                                                 aggcgaagTGGCGGGCGGGGGGGGGATGCTTCGATCCGGCAAAAAAGCCCGAGGAGTGAGCGCCGCCGTGTAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTGGGGAAAGTTACAACATAAGAACGCAGCCAGCCCAtt                 >  2:1084514/1‑151 (MQ=255)
                                                                                                                                     gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc             <  2:361518/151‑1 (MQ=255)
                                                                                                                                             gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc     >  2:265715/1‑151 (MQ=255)
                                                                                                                                                  ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc  >  2:531447/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599427‑599721

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