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

AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599427‑599721
                                                                                                                                                      |                                                                                                                                                   
aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGggtg                                                                                                                                                    <  1:386148/151‑1 (MQ=255)
   gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggggg                                                                                                                                                 >  2:395718/1‑149 (MQ=255)
    cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGTTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggg‑gggggggg                                                                                                                                                  >  2:404791/1‑148 (MQ=255)
       gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggg‑gggggggggggcc                                                                                                                                             >  1:628325/1‑145 (MQ=255)
          gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggggggtac                                                                                                                                             >  2:321696/1‑142 (MQ=255)
           acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCg                                                                                                                                         <  1:724377/151‑1 (MQ=255)
           acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCg                                                                                                                                         <  1:955680/151‑1 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggggggagatcgg                                                                                                                                         >  1:699491/1‑139 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGa                                                                                                                                                <  2:699468/142‑1 (MQ=255)
             tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCg                                                                                                                                         >  1:1876/1‑149 (MQ=255)
                     caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGGACCGCG‑ATcc                                                                                                                                    >  2:1188618/1‑145 (MQ=255)
                      aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑Ggggggggggttcgtcggt                                                                                                                                       >  2:1257891/1‑130 (MQ=255)
                      aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCTGCG‑GTCCCGCaaa                                                                                                                              >  2:1051492/1‑150 (MQ=255)
                           tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGCTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTTCTGCG‑GTCCCGCGAAAAAGcc                                                                                                                        >  2:1125281/1‑151 (MQ=255)
                            ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTCCG‑GCCCGGCAAAAAc                                                                                                                           >  2:21466/1‑147 (MQ=255)
                                 aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAACGCCCg                                                                                                                      >  1:397589/1‑147 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGt                                                                                                                <  1:1091354/151‑1 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGt                                                                                                                <  2:760721/151‑1 (MQ=255)
                                    cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGTGGCGGTCCTGCG‑GACTCGCCAAAAAGCCCGAGGAGTg                                                                                                               >  1:494166/1‑151 (MQ=255)
                                     aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAggg                                                                                                               >  1:496351/1‑148 (MQ=255)
                                      gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                             <  2:222616/151‑1 (MQ=255)
                                        gAAAATGGCGTGCCCGGTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                            <  2:106039/150‑1 (MQ=255)
                                          aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTTCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGATcgc                                                                                                          >  1:620731/1‑150 (MQ=255)
                                            aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGAGCCCCCG‑GGCCCGCAAAAAAGACCGCGaaagggtcccc                                                                                                         >  2:1234014/1‑140 (MQ=255)
                                             tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTTCTCCG‑GTCCCGGAAAAACGCCCGGGGAGGGATTgccgcc                                                                                                      >  2:470396/1‑151 (MQ=255)
                                             tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGATCC‑GCG‑GTCCGGCAAAAACGCCCGAGAAGGGAGCgccgcc                                                                                                      >  2:865804/1‑151 (MQ=255)
                                                 gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑GTGCCGCAAAAACACCCGAGAAGAGGTcgccg                                                                                                        >  2:688396/1‑145 (MQ=255)
                                                   gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                 <  2:416664/150‑1 (MQ=255)
                                                    cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                 <  1:423071/149‑1 (MQ=255)
                                                    cccGTTTATCGGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                               <  2:772753/151‑1 (MQ=255)
                                                      cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTGCTTCGCAGCGGGCAAAAAC‑CCCGGGAAATGCCCCTCGCCGGGAAgg                                                                                               >  2:895933/1‑149 (MQ=255)
                                                          tATCTGCGCGATGTCGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGGCCCTCG‑TGCCGGCAAAAAAGGCCCCGGAGGGATCgccgccccggagtgcaaa                                                                                          >  2:81725/1‑138 (MQ=255)
                                                           aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                         <  2:804859/150‑1 (MQ=255)
                                                                cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCG‑ATCCCGCAAAAACGCCCGAGAAGTGATTGCCGCCGTGGAGGGCAAAA‑Tgg                                                                                     >  1:601198/1‑148 (MQ=255)
                                                                  cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCCGCG‑ATCCGGCCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAAc                                                                                 >  1:1177137/1‑150 (MQ=255)
                                                                   gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCACTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAGTGGTCGCCGCCGGGAAGGACCAAA‑TGGAAAAGc                                                                               >  1:2483/1‑151 (MQ=255)
                                                                     tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCGCG‑ATCCGGGCAAAACGCCCGAGAAGTGAGTGCCGCCGGGAAGGACAAAA‑TGGgaaacctg                                                                             >  2:440457/1‑146 (MQ=255)
                                                                        tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGaa                                                                           >  1:919200/1‑150 (MQ=255)
                                                                          cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGTCTGCG‑CACCGGCAAAAACCGCCGCGGAGTGATCGGCCACGGGAAGGGCAAAC‑TCGAAAAGCCGcaaa                                                                         >  1:687562/1‑150 (MQ=255)
                                                                            aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑GACCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGAGAAAC‑TGGaaa                                                                                  >  2:428353/1‑139 (MQ=255)
                                                                             aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAagtag                                                                     <  2:1251972/151‑1 (MQ=255)
                                                                               gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAagtag                                                                     >  1:160052/1‑149 (MQ=255)
                                                                                 ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCCGCG‑GACCCGCCAAAACGCCCGAGAAAGGATCGCCGCCGTGGAGGGCCAACATGGAAAcccgaaaaataatc                                                                    >  2:385991/1‑136 (MQ=255)
                                                                                  cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCGGc                                                                <  2:723300/151‑1 (MQ=255)
                                                                                      ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCgg                                                             >  1:172987/1‑150 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGa                                                            <  2:566758/150‑1 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCCGCAAAAACGCCCGAGAAGGGATCGGCGCCGGGGAGGACAAAA‑TGGAAAAGCTGGAAAGGAGTCTGCCgg                                                             >  1:551999/1‑149 (MQ=255)
                                                                                                gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                  <  2:1038947/151‑1 (MQ=255)
                                                                                                  tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                 <  1:948819/150‑1 (MQ=255)
                                                                                                        ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCCCG‑ATCCCGCCAAAAACCCCCAGAAGGGGTCGCCCCCCTGGAGGGCAAAC‑TGGAAAAGCTGAAAAGTAGTCTGCCCGGAGGCGTGGAGATAGtt                                            >  2:250471/1‑149 (MQ=255)
                                                                                                           gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                       <  1:1209215/151‑1 (MQ=255)
                                                                                                           gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                       <  2:208740/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                     <  2:142309/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                     <  1:331456/151‑1 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCCTGGAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTacaaca                                      >  1:970545/1‑150 (MQ=255)
                                                                                                                gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                  <  2:445364/151‑1 (MQ=255)
                                                                                                                gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑GTCCGGCAAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAAA‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGGTAGTTACAACATa                                    >  1:129954/1‑149 (MQ=255)
                                                                                                                   gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                <  1:1028587/150‑1 (MQ=255)
                                                                                                                   gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                               <  2:778444/151‑1 (MQ=255)
                                                                                                                    aaCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATc                               >  1:33991/1‑150 (MQ=255)
                                                                                                                         aaCGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATACTGCG‑ATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATAGTTACAACATACCATCGcagcc                         >  1:924395/1‑151 (MQ=255)
                                                                                                                          aCGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGGAGTGATTGCCGCCGGGAAGGACAAAC‑TGGAAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGGTAGTTACAAAAAAAGATCGcagcc                         >  2:1243017/1‑150 (MQ=255)
                                                                                                                           cGGCGAAGGTGAAGTGGCGGGC‑GGGGGGGTGGTCCTGCG‑ATCCCGCAAAAAAGCCCGCGAAGGGATCGCCGCCGGGGAGGACAAAC‑CGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGAAAGTTACAACATACGATCGcagcc                         >  2:575685/1‑149 (MQ=255)
                                                                                                                           cGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                       <  1:381187/151‑1 (MQ=255)
                                                                                                                                aaggcgaagTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                  <  2:405200/151‑1 (MQ=255)
                                                                                                                                 aggcgaagTGGCGGGC‑GGGGGGGGGATGCTTCG‑ATCCGGCAAAAAAGCCCGAGGAGTGAGCGCCGCCGTGTAGGACAAAC‑TGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTGGGGAAAGTTACAACATAAGAACGCAGCCAGCCCAtt                 >  2:1084514/1‑151 (MQ=255)
                                                                                                                                     gaagTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc             <  2:361518/151‑1 (MQ=255)
                                                                                                                                             gggc‑gggGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc     >  2:265715/1‑151 (MQ=255)
                                                                                                                                                  ggggTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc  >  2:531447/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  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: ‑