Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5790TG40.9% 104.7 / 15.3 67V289G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (5/34);  new base G (27/0);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 8.22e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599433‑599727
                                                                                                                                                  |                                                                                                                                                    
cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                    <  1:566924/149‑1 (MQ=255)
 gacgtcTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGggtggt                                                                                                                                                      >  1:343820/1‑146 (MQ=255)
      cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                              <  1:154259/149‑1 (MQ=255)
      cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCt                                                                                                                                                >  1:254194/1‑147 (MQ=255)
          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc                                                                                                                                         <  2:479730/150‑1 (MQ=255)
          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCg                                                                                                                                        <  2:1062199/151‑1 (MQ=255)
              acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaa                                                                                                                                    <  1:326581/151‑1 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGTCCCGCaaaaa                                                                                                                                 >  1:292107/1‑150 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCCGCGGTCCGGCaaaaa                                                                                                                                 >  2:121924/1‑150 (MQ=255)
                      ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGAGCCCGCAAAAAAGccc                                                                                                                            >  1:379732/1‑151 (MQ=255)
                      ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                             <  1:350362/150‑1 (MQ=255)
                      ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGCGATCCCGCgaaaaaaccc                                                                                                                            >  2:669350/1‑146 (MQ=255)
                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGccc                                                                                                                            >  1:424832/1‑144 (MQ=255)
                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTGCCGCGGTCCGGGAAAAAAGGCCGGGAAggg                                                                                                                    >  1:86874/1‑148 (MQ=255)
                                 tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTCCGGTCCCGCAAAAACGCCCgcg                                                                                                                         >  2:1066257/1‑141 (MQ=255)
                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCTGCGATCCCGGAAAAAAGCCCGCGAAGTGAAcg                                                                                                                >  2:59007/1‑151 (MQ=255)
                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCGAGGAGGGATcg                                                                                                                >  1:1173366/1‑151 (MQ=255)
                                   aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                                                                               <  2:986004/151‑1 (MQ=255)
                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGATCCGGCCAAAAAGCCCGAGAAGGGATCgccgcc                                                                                                           >  2:195686/1‑150 (MQ=255)
                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  2:1035472/150‑1 (MQ=255)
                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGCTCCGGCAAAAACGCCCGCGGAGGGATCgccgcc                                                                                                           >  1:715395/1‑150 (MQ=255)
                                            tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTCCGGTCCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGt                                                                                                         >  1:47495/1‑148 (MQ=255)
                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCgccgcc                                                                                                           >  1:744280/1‑145 (MQ=255)
                                               ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                                   <  1:339372/151‑1 (MQ=255)
                                                 gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATTCGGCAAAAAACCCCGAGAAGTGATCGCCGCCGGGCAGGAc                                                                                                  >  1:653538/1‑150 (MQ=255)
                                                        tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTTCGGCCAAAAAGCCCGAGAAGGGATCGCCGCCGCGGAGGGGAAACTgg                                                                                           >  2:1012921/1‑150 (MQ=255)
                                                         gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                          <  2:128973/150‑1 (MQ=255)
                                                            cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACg                                                                                      <  1:792168/151‑1 (MQ=255)
                                                              aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                    <  1:1105432/151‑1 (MQ=255)
                                                                gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGc                                                                                     >  1:490670/1‑148 (MQ=255)
                                                                   gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                               <  1:904605/151‑1 (MQ=255)
                                                                    cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGAGCGCCGCCGGGGAGGGAAAACTGGAAAAGCCGgaaa                                                                               >  2:211926/1‑150 (MQ=255)
                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  1:1201020/151‑1 (MQ=255)
                                                                         gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                         <  1:781410/151‑1 (MQ=255)
                                                                         gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCTCGGTCCGGCAAAAAAGCCCCGGGAGTGATCGGCGGCGTGAAGGACAAACTGGAAAAGCGGAAAAATATTc                                                                         >  2:581662/1‑151 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCCAGAAGGGAGCGGCGCCGGGAAGGGAAAAATGGGAAAGCCGAAAAGGAGTCTGcc                                                                     >  2:175541/1‑151 (MQ=255)
                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCTAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                   <  1:663298/151‑1 (MQ=255)
                                                                                 tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                  <  2:1117487/150‑1 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  2:325612/151‑1 (MQ=255)
                                                                                    cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                              <  1:13070/151‑1 (MQ=255)
                                                                                     ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                             <  1:1127641/151‑1 (MQ=255)
                                                                                     ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGATCCGGCAAAAACCCCCGAGAAAGGATCGCCGCCGTGAAAGGCAAACAGGAAAACCTGAAAAGGAATCTGCCGGaa                                                                 >  1:745036/1‑147 (MQ=255)
                                                                                          gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                        <  2:955922/151‑1 (MQ=255)
                                                                                          gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                        <  1:1136414/151‑1 (MQ=255)
                                                                                           aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGCTCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAATAGTCTGCCGGAAGGCGTgg                                                          >  1:378009/1‑148 (MQ=255)
                                                                                            tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                       <  2:671209/150‑1 (MQ=255)
                                                                                             gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                      <  1:801933/150‑1 (MQ=255)
                                                                                             gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                     <  1:388431/151‑1 (MQ=255)
                                                                                              cgcCGGGGCAGTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGGCCTACGAACCCGCAAAAAGAGCCGGGAAAGGAGTGCCGCCGCGAAGGGCAAACCGGAAAAGCCGGAAAGTAGTATGCCgggaagag                                                             >  2:162457/1‑136 (MQ=255)
                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATTGCCGCCGGGAAGGAGAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGGGGGGAGAAAGtt                                                  >  2:956159/1‑150 (MQ=255)
                                                                                                  ggggCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGGGGGGGACCCGCCATCCCGCAAAAACGCCCGCGGGGCGAAGCCCCCCCTGCAGGACAAACAGGAAACGCTGAAAAGTAGGCCGGCGGAAGGGGTGGGGAGAGTTac                                                >  2:753703/1‑151 (MQ=255)
                                                                                                     gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCCGGGAGGCGACCCTGCGGTTGGGGAAAAACTCCCGGGAGGAGGTACCCGCCGCGGAGGAGAAAAAGGAAAAGAAGAAAAGAAGGGTGCCGGTAGGCGTGGAGGTAgat                                                  >  1:173925/1‑144 (MQ=255)
                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                           <  1:669383/151‑1 (MQ=255)
                                                                                                         tGCCGACCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACTGCGAACGCGCCAAAACGCCCGAGGCGTGAGCGCCCCCGGGCAGGGAAAAACGGAAAAGATGACAAAAAGGGTGCCGGCAGGGGCGGGGATATTTACAacacac                                         >  2:965549/1‑151 (MQ=255)
                                                                                                          gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                        <  2:947982/151‑1 (MQ=255)
                                                                                                             gagcTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                          <  2:303385/143‑1 (MQ=255)
                                                                                                             gAACTAAACGGTGAAGGGGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                      <  2:876068/150‑1 (MQ=255)
                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGGAACGCCGAAAAGTAGTCTGCCGGGAGGGGGGGGGGTAGGTTCAACATa                                          >  2:694344/1‑145 (MQ=255)
                                                                                                                cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATACTCCGATCCCGCCAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGGTGAAAAGCAGTCTCCCGGAAGGCGGGGGGGTAGGTATAACATaaggacgcc                                  >  2:879194/1‑143 (MQ=255)
                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGATCCGGCAAAAAAACCCGGGGAGGGAATGCCGCCCTGAAGGACAAACCGGGAACGGTGAAAAATAGTCTGCCGGGAGGGGGGGGGGTAGTTACAACATACGATCGcag                                 >  1:7734/1‑150 (MQ=255)
                                                                                                                   aaCGGCGAAGGCGAAGTGGCGGGCGGGGAGGTGATCCTGCGGTCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACCAACTGGAAACGCTGCAAAGTAGGCTGCCGGGAGGCGTGGAGAAAGTTAAAACATa                                          >  1:712965/1‑140 (MQ=255)
                                                                                                                       gcgaaggcgaagTGGCGGGGCGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGGCAAACTGGAAAAGCTGGAAAGAAGTCTGCCGGAAGGCGGGGGGATAGGTAAAACATACGATCGcagcc                               >  2:255507/1‑147 (MQ=255)
                                                                                                                           aggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                       <  1:455309/151‑1 (MQ=255)
                                                                                                                            ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                      <  1:943189/151‑1 (MQ=255)
                                                                                                                                       gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc           <  2:361552/151‑1 (MQ=255)
                                                                                                                                        ggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg          <  2:729080/151‑1 (MQ=255)
                                                                                                                                        ggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg          <  2:903217/151‑1 (MQ=255)
                                                                                                                                           ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc        >  2:789243/1‑150 (MQ=255)
                                                                                                                                                ggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  <  1:119458/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599433‑599727

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: