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

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725
                                                                                                                                                   |                                                                                                                                                       
cgcgcgACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                                                                                      <  2:76617/151‑1 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggg                                                                                                                                                       >  1:175949/1‑149 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                     >  1:146047/1‑149 (MQ=255)
   tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                       >  1:710767/1‑147 (MQ=255)
     gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggggcc                                                                                                                                                 >  2:186022/1‑145 (MQ=255)
     gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATAGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGTAGTGGCGGTCGGGGGGGCGGTcc                                                                                                                                                 >  2:473344/1‑151 (MQ=255)
      acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                                 >  2:754821/1‑150 (MQ=255)
           tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAACCTGC‑g                                                                                                                                            >  1:746054/1‑149 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑Ga                                                                                                                                           <  1:166213/149‑1 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCTGC‑GAttc                                                                                                                                        >  2:1151380/1‑149 (MQ=255)
              aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTTC‑GGT‑cc                                                                                                                                       >  1:17615/1‑150 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGc                                                                                                                                    <  2:448397/150‑1 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCa                                                                                                                                   <  2:767784/151‑1 (MQ=255)
                      cGTTTTATAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGCTGTTGCGAAGGCCTAGATTGGCCCGGAGATGCCCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGCCGGGGGGGGGGCCTGG‑CAT‑CCCGCcaaa                                                                                                                                >  2:1194911/1‑149 (MQ=255)
                      cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAAc                                                                                                                              <  1:969828/151‑1 (MQ=255)
                      cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GAT‑CCGGCaaaaa                                                                                                                               >  2:719599/1‑150 (MQ=255)
                       gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACg                                                                                                                             <  1:471193/151‑1 (MQ=255)
                       gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAATCTCC‑GAT‑CCGGGgaaaa                                                                                                                               >  1:441914/1‑149 (MQ=255)
                           tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGC‑GAT‑CCGGCCAAAAAGccc                                                                                                                          >  1:243535/1‑150 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCg                                                                                                                         <  2:614246/150‑1 (MQ=255)
                             aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGACCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGGTCCTGC‑GCA‑ACGGCcaaaaaac                                                                                                                            >  1:985008/1‑143 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACT‑AACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGC‑GAT‑CCCGCAAAAACGCCCggg                                                                                                                       >  1:428317/1‑147 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGG‑CAA‑CCCGCAAAAACCACCGCGaa                                                                                                                     >  1:644390/1‑149 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGt                                                                                                                   <  1:499012/151‑1 (MQ=255)
                                  cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTT‑CGGGCAAAAAC‑CCCGCGAAggg                                                                                                                  >  1:311260/1‑149 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  2:1160333/151‑1 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAT‑CCCGCAAAAAAGCCCGGGAAggg                                                                                                                  >  1:825831/1‑147 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  1:815327/151‑1 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               <  1:113450/150‑1 (MQ=255)
                                        aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGC‑GGT‑CCGGCAAAAAAGCCCGGGAAGTGATcgcc                                                                                                            >  1:245863/1‑151 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                        <  2:1084499/151‑1 (MQ=255)
                                                tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAT‑CCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTg                                                                                                      >  1:148087/1‑149 (MQ=255)
                                                 gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAT‑CCGGCCAAAACGCCCGAGAAGTGATCGccgccggg                                                                                                      >  2:841360/1‑146 (MQ=255)
                                                        tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             <  1:1218851/150‑1 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                            <  2:1171040/150‑1 (MQ=255)
                                                                  aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                  <  1:1253400/151‑1 (MQ=255)
                                                                            ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGC‑GAT‑CCCGCAAAAACCCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAAAGGAAACGCTGaaaaaaag                                                                         >  2:777301/1‑146 (MQ=255)
                                                                             gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                       <  1:831252/151‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      <  1:90683/150‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                     <  1:849784/151‑1 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGG‑GGT‑CCGGCAAAAACGCCCGAGAAGAGGTCGCCGCCGTGAAGGGCAAAAATGAAAAGCTGgaaa                                                                             >  2:541705/1‑142 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGC‑GGT‑CCGGCCAAAACACCCGAGACGTGGTGGCCGCCGGGAAGGAGCAACTGGGAACCCTCaa                                                                               >  2:68643/1‑140 (MQ=255)
                                                                                 aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                   <  1:877165/151‑1 (MQ=255)
                                                                                     tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                <  1:1005083/150‑1 (MQ=255)
                                                                                        cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGC‑GATCCCGGAAAAAACCCCCG‑GAAATGATCGCCGCCGTGAAGGACAAACTGGAAAACCTGAAAAGAAGACTGCCGGAAgg                                                             >  1:777164/1‑150 (MQ=255)
                                                                                         ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGAGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  1:943790/150‑1 (MQ=255)
                                                                                            gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                        <  1:553282/151‑1 (MQ=255)
                                                                                            gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGACCCTGC‑GAT‑CCGGCAAAAAAGCCCGCGAAGTGGGCGCCCGCCCGAAGGGGAAAATTGGAACGGTGAAAAGTAGTGTGCCGGCAggag                                                           >  2:561265/1‑146 (MQ=255)
                                                                                                 gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  1:988726/150‑1 (MQ=255)
                                                                                                  cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGC‑GAT‑CCCGCAAAAACGCCCCCGAAGGGATCGCCGCCGCGGAGGGCAAACTGGAAAAGCTGAAAAGGAGTCTGCCCGGAGGCgggggg                                                      >  2:527003/1‑145 (MQ=255)
                                                                                                    ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                                 <  2:235800/150‑1 (MQ=255)
                                                                                                       gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                             <  1:141980/151‑1 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                         <  1:76273/151‑1 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                         <  2:311247/151‑1 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGAGGGGGTGATGCCGT‑GAT‑ATGGCGAAAAAGGCCGGGAAATGAGCGCCGCCGAGAAGGGGAAACAGATAACACCGCAGAGGATTTTGCCGGAAGGCGGGGGGGTAGGTaccaca                                          >  2:990934/1‑150 (MQ=255)
                                                                                                             tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAT‑CCGGGAAAAAAGCCCGAGAAGTGCTGGCCGCCGTGAAAGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTACAACATa                                        >  1:918007/1‑150 (MQ=255)
                                                                                                                cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAT‑CCCGCAAAAACGCCCGAGAAGGGGTCGCCGCCGTGGAGGAAAAAACGGAAACGCTGAAAAGAAGGCTGCCGGAAGGCGGGGGGATAGTTACAAAATACGa                                     >  2:162615/1‑150 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  2:429631/151‑1 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  1:860854/151‑1 (MQ=255)
                                                                                                                    cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGC‑GAT‑CCGGCAAAAACGCCCGAGAAGGGAACGCCGCCGTGAAGGACAAACCGGAAAAGCTGAAAAGTAGGCTGCCGGAAAGGGGGGAGAAAGTTACAACATACGATCGc                                 >  1:809790/1‑150 (MQ=255)
                                                                                                                      aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                <  2:1010212/149‑1 (MQ=255)
                                                                                                                        aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GAT‑CCGGCCAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGcagcc                             >  1:1222866/1‑150 (MQ=255)
                                                                                                                           gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCt                         <  2:746026/151‑1 (MQ=255)
                                                                                                                              aaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                      <  1:792000/151‑1 (MQ=255)
                                                                                                                               aggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGGAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  2:995184/151‑1 (MQ=255)
                                                                                                                               aggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  2:907325/151‑1 (MQ=255)
                                                                                                                                ggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  2:17614/150‑1 (MQ=255)
                                                                                                                                    aagTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg                <  1:562027/151‑1 (MQ=255)
                                                                                                                                           gggcgggGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         <  2:428301/151‑1 (MQ=255)
                                                                                                                                             gcgggGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg        >  1:484239/1‑150 (MQ=255)
                                                                                                                                                gggTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa    <  2:710739/151‑1 (MQ=255)
                                                                                                                                                  gtggtgATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  1:187127/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725

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