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

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599429‑599724
                                                                                                                                                    |                                                                                                                                                       
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                                                                                      <  1:292810/151‑1 (MQ=255)
    cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                    <  2:976262/149‑1 (MQ=255)
    cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGAGCCCGTTTATCTGCGCGATGATGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGCGTGGGAGTGAt                                                                                                                                                   >  2:1050812/1‑150 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggatctgcgggc                                                                                                                                          >  1:776699/1‑138 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTGCTGCGGTc                                                                                                                                         >  2:29236/1‑151 (MQ=255)
              aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGGATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggggccttgcagcc                                                                                                                                         >  1:1159682/1‑136 (MQ=255)
              aaTCACATCGTTTAAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc                                                                                                                                        <  1:105737/150‑1 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCgg                                                                                                                                      >  1:157980/1‑149 (MQ=255)
                  acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTcccg                                                                                                                                      >  1:165921/1‑146 (MQ=255)
                   caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCaa                                                                                                                                   >  2:1038546/1‑150 (MQ=255)
                           tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                           <  1:809044/150‑1 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGACCCGGCAAAACCCCCCgg                                                                                                                         >  2:719323/1‑149 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                      <  2:1168770/151‑1 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCaaaaaacgcc                                                                                                                           >  1:897317/1‑141 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                    <  1:1087557/150‑1 (MQ=255)
                                  cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                                  <  2:967158/151‑1 (MQ=255)
                                   aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                 <  2:619012/151‑1 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  2:251851/151‑1 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  2:604762/151‑1 (MQ=255)
                                        aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGCCAAAACGCCCGAGAAGGGGTcgcc                                                                                                            >  1:358824/1‑151 (MQ=255)
                                           tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGTTCCGCGAAAAAGCCCCGGAAGGGGGcgccggc                                                                                                         >  1:400403/1‑149 (MQ=255)
                                                  cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCGGTTCGGCGAAAACGCCCGAGAAGTGATcgccgc                                                                                                          >  2:238753/1‑143 (MQ=255)
                                                   ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCCTCGGGACGGCCAAAACGCCCGGGAAGCGATGGCCGCCGGGAAg                                                                                                   >  2:1027429/1‑149 (MQ=255)
                                                   ccGTTTATCGGCGCGAGGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                                 <  1:755971/151‑1 (MQ=255)
                                                    cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGACCGGCAAAAACGCCCGAGAAGAGGTCGccgccgcggag                                                                                                   >  1:315045/1‑143 (MQ=255)
                                                       ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                             <  1:409589/150‑1 (MQ=255)
                                                       ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                             <  1:853747/150‑1 (MQ=255)
                                                       ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACaa                                                                                             >  2:1147452/1‑150 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCCAAAACGCCCCGGAAGGGGTCGCCCCCGTGAAGGACaaa                                                                                            >  2:154456/1‑149 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                           <  1:872878/150‑1 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAACCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACaa                                                                                             >  2:1063985/1‑148 (MQ=255)
                                                           cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGGGAAGGACAAAATgg                                                                                        >  2:252922/1‑151 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                          <  2:314266/147‑1 (MQ=255)
                                                                 ggtGTTGCGAAGGGACAGAGTGGCCCGGGGATGCGCCGGGGGGTGGCCGAGCGAAACGGCGAAGGCGAAGCGGCGGGCGGGGTGGTGATCCTGCGATCCGGAAAAACCGCCCGAGAGGGGATCGCCGCCGTGAAGGACAAACTGGAAACGc                                                                                  <  1:1196948/149‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                <  2:328031/150‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATTCGGCAAAAACACCCGGGGAGTGAGCGCCGCCGGAAAGGGGAAAcacgaagcacgcg                                                                               >  1:1196683/1‑139 (MQ=255)
                                                                      tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAGCCGGCGAAAACGCCCCAGAAGAGATCGCCGCCGGGAAGGGCAAACTGGAAACGCGGaa                                                                              >  1:500340/1‑150 (MQ=255)
                                                                        cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGATCCCCGATCCGGCAAAAAACCCCCAGAAGGGAGTGCCGCCGGCAAAGAAAAACAAGaaaacccggaa                                                                             >  1:702348/1‑141 (MQ=255)
                                                                           aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                        <  2:590894/151‑1 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGCGAAGGAAAAACTGGAAAAGCCGAAAAGTAGTCt                                                                     >  1:804281/1‑149 (MQ=255)
                                                                                 aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAACCTGCGGTCCGGCAAAAACGCCCGAGAAGGGAACGCCGCCGTGAAAGACAAACTGGAAACGCTGGAAagtagt                                                                       >  2:71619/1‑146 (MQ=255)
                                                                                  gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGTCTGcc                                                                  >  2:57812/1‑150 (MQ=255)
                                                                                   aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATTGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGAAATCTGCCgg                                                                >  2:482031/1‑151 (MQ=255)
                                                                                   aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCAAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAAAACCCCCGAAGTGAACGCCGCCCGGAAGGACAAACTGGAAAAGCTGAAAAGCAGGCTGcc                                                                  >  2:1052830/1‑149 (MQ=255)
                                                                                     tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                               <  2:181286/150‑1 (MQ=255)
                                                                                       gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                            <  2:48729/151‑1 (MQ=255)
                                                                                        cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGACCTGCGATCCCGCCAAAACCCCCGCGGAAGGAACGCCGCCGGGGAGGGCCAACAGGAAACTCTGAAAAGTAGTcacc                                                                   >  1:1161035/1‑140 (MQ=255)
                                                                                        cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGGAGTGATCGCCGCCGGGGAGGACAAAATGGAAACGCTGAAAAAGAGGCTGCCGGGAgg                                                            >  1:851529/1‑150 (MQ=255)
                                                                                            gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                       <  2:321264/151‑1 (MQ=255)
                                                                                            gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGAGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAGCGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGGCTGCCGGAAg                                                             <  1:384589/145‑1 (MQ=255)
                                                                                              gaTGCGCCGGGGCACTGCCGAACTAAACGGCGAAGGAGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                     <  1:406669/151‑1 (MQ=255)
                                                                                               aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCCAAAAAGCCCGAGAAGGGATTGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCGGGAGGCgggg                                                       >  2:1184600/1‑148 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGGCTGCGATCCCGCCAAAAAGCCCCGGAAGGGGAGCCCCGCCGGAAGGGCAAAATTGGAAAGGCGGCAAGGAGGCTGCCGGAAGGCgggggg                                                     >  1:309295/1‑144 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                                <  2:135737/151‑1 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                                <  2:759230/151‑1 (MQ=255)
                                                                                                    ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                               <  1:1065460/151‑1 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                        <  2:930537/151‑1 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                        <  1:842402/151‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                      <  1:143059/150‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                      <  1:1129339/150‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                     <  2:711419/151‑1 (MQ=255)
                                                                                                                              aaggcgaagTGGCGGGCGGGGTGGGGACCCTGCGATCCCGCAAAAACGCCCGAGAAGGGATCGCCGCCGCGGAGGGCAAACTGGAAACGCCGAAAAGAAGTCTGGCCGGAGGCGGGGGGATAGATTCAAAAAACGacagca                               >  1:1110951/1‑136 (MQ=255)
                                                                                                                               aggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAAAGCCCGAGAAGGGATCGCCGGCGGGAAGGAAAAACTGGAAACGCTGAAAAGTAGGCTGCCGGAAGGCGTGGAGGTAGTTATAACATAAGATCGCAGCCAGCTCaat                    >  1:771692/1‑149 (MQ=255)
                                                                                                                                   gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAca                <  1:29232/151‑2 (MQ=255)
                                                                                                                                   gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                 <  1:489377/150‑1 (MQ=255)
                                                                                                                                    aagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg               <  1:1050748/151‑1 (MQ=255)
                                                                                                                                         gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa          <  1:416083/151‑1 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGACGTGATCGCCGCCGCGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc        >  2:91152/1‑151 (MQ=255)
                                                                                                                                               cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc     <  2:1132781/151‑1 (MQ=255)
                                                                                                                                               cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc     >  1:955892/1‑151 (MQ=255)
                                                                                                                                                ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc     >  1:113280/1‑150 (MQ=255)
                                                                                                                                                 gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa   <  2:1090594/151‑1 (MQ=255)
                                                                                                                                                   gtgttgATCCTGCGATCCGGCAAAAACGACCGAGAAGTGATCGCCGCCGTGCAGGACAATCTGGAAAGGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGCTCGCAGCCAGCTCATTGACCGCGCCATCGACAAc  >  1:1034081/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599429‑599724

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: