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

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725
                                                                                                                                                   |                                                                                                                                                      
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                      >  2:429153/1‑150 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                    >  2:825092/1‑149 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                    <  2:527438/151‑1 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                    <  1:977783/151‑1 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                    <  2:382520/151‑1 (MQ=255)
  cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                    >  1:1114210/1‑148 (MQ=255)
     gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                                                                                     >  1:220460/1‑145 (MQ=255)
       cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAAcc                                                                                                                                                >  2:132461/1‑149 (MQ=255)
        gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggggtc                                                                                                                                                >  1:227668/1‑142 (MQ=255)
           tttAATCACATCGTTTTAAAAGCCAGCGAAAATGGCGTGCACGTTTATCTGCGCGATGTTGCAAAGGTCCAGATTGGCCCGGAGATGGGCCGGGGCATTGCCGGACTAAACGGCGAAGGCGCAGTGGCGGGGGGGggggtg                                                                                                                                                    >  1:739228/1‑141 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCggt                                                                                                                                          >  2:934176/1‑148 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggggatccgcggttc                                                                                                                                        >  1:2562/1‑137 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATc                                                                                                                                         <  1:8643/150‑1 (MQ=255)
               ctCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGGACTTAACGGCGAAGGCGGAGTGGCGGGCGGGGGGGGGGGCCTGCGGTCC‑gg                                                                                                                                     >  1:410621/2‑151 (MQ=255)
               aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTcccg                                                                                                                                      >  2:735675/1‑149 (MQ=255)
                           tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGCTCC‑GGGAAAAAAGccc                                                                                                                          >  2:874259/1‑150 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCga                                                                                                                        <  2:241971/151‑1 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCg                                                                                                                         <  1:230700/150‑1 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGAGGTGGTCCTTCGCTCC‑CGGCAAAACGccc                                                                                                                          >  2:349628/1‑149 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCC‑GGGAAAAACGccc                                                                                                                          >  1:190468/1‑147 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGGTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGaa                                                                                                                     <  1:3215/151‑1 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGGAAAAACGccc                                                                                                                          >  1:201542/1‑146 (MQ=255)
                                gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCC‑GGCAAAAACGccc                                                                                                                          >  2:998169/1‑145 (MQ=255)
                                 ccAGGGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCGGCGATCC‑GGCAAAAACGCCCGAGAAgg                                                                                                                   <  1:346483/151‑2 (MQ=255)
                                  cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTg                                                                                                                  <  1:852284/151‑1 (MQ=255)
                                     tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGGGCTCCGGGCAAAAAcccc                                                                                                                           >  1:604976/1‑140 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATc                                                                                                               <  2:993038/150‑1 (MQ=255)
                                                tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                        >  2:1068357/1‑147 (MQ=255)
                                                 gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                    <  2:104193/150‑1 (MQ=255)
                                                 gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                    <  2:129963/150‑1 (MQ=255)
                                                      tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCGGGGAGTGATCGCCGCCGGGGAgg                                                                                                  >  1:298206/1‑147 (MQ=255)
                                                        tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             <  1:396375/150‑1 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                        <  2:369604/150‑1 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  2:679377/151‑1 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  2:502127/151‑1 (MQ=255)
                                                                 gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATTC‑GGCAAAAAACCCCGGGAAGTGATCGCCGCCGCGAAGGACAAAATGGAGACg                                                                                    >  2:786398/1‑150 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 <  2:57978/150‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 <  1:956131/150‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 >  2:917136/1‑150 (MQ=255)
                                                                       gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                             <  1:241918/151‑1 (MQ=255)
                                                                         gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  1:776174/151‑1 (MQ=255)
                                                                         gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:478653/151‑1 (MQ=255)
                                                                         gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                             >  2:840243/1‑149 (MQ=255)
                                                                             gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                         <  1:243215/149‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      <  2:627953/150‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      <  1:747499/150‑1 (MQ=255)
                                                                                  gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGGAACGCTGAAAAGTAGTCTGcc                                                                   >  2:903444/1‑150 (MQ=255)
                                                                                     tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                <  1:286883/150‑1 (MQ=255)
                                                                                       gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCTGCGATCC‑GGCAAAAAAACCCGAGAAATGAACGCCGCCGGGAAAGAAAAACTGGAAAAGCTGAAAAAGAGTCTGCCCGCAgg                                                             >  2:816905/1‑151 (MQ=255)
                                                                                           ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                          <  1:11120/150‑1 (MQ=255)
                                                                                           ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                         <  2:437591/151‑1 (MQ=255)
                                                                                           ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAAAGCCCGAGAAGGGATTGCCGCCGGGAAGGACAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCggg                                                         >  2:929198/1‑149 (MQ=255)
                                                                                            gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAACC‑GGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                        >  2:571838/1‑151 (MQ=255)
                                                                                              gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                      <  2:222339/151‑1 (MQ=255)
                                                                                              gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  2:46575/150‑1 (MQ=255)
                                                                                                 gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  1:721306/150‑1 (MQ=255)
                                                                                                  cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATTC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGCGAAGGACAAACTGGGAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGgag                                                      >  1:406927/1‑147 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAAAGCCCGAGAAGTGATCGGCGCCGTGGAGGACAAACTGGGAACGCCGGAAAGGAGGCTGCCCGGAGGCGTGGGGAt                                                    >  2:797157/1‑148 (MQ=255)
                                                                                                      ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                              <  2:881857/151‑1 (MQ=255)
                                                                                                        ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                            <  1:602706/151‑1 (MQ=255)
                                                                                                         gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTacaa                                            >  2:549705/1‑150 (MQ=255)
                                                                                                               ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                      <  2:196911/150‑1 (MQ=255)
                                                                                                               ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑CGCAAAAACGCCCGCGAAGGGATCGCCGCCGTGAAGGGCAAACTGGGAAAGCTGGAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTACAACAt                                         >  1:1008523/1‑147 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  2:912159/151‑1 (MQ=255)
                                                                                                                 gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAACCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAAATGGAAAAGCTGAAAAGAAGTCTGCCGGAAGGGGGGGGGGAAGTTacaacaaa                                        >  1:736242/1‑144 (MQ=255)
                                                                                                                             gaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGCAGCCAGCTc                        >  2:37072/1‑150 (MQ=255)
                                                                                                                              aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                      <  1:1001924/151‑1 (MQ=255)
                                                                                                                                   gaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                 <  2:220500/151‑1 (MQ=255)
                                                                                                                                      gTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGGTAGTTACAACAAACGATCGCCGCCAGCTCATTTAcc                 >  2:364735/1‑148 (MQ=255)
                                                                                                                                         gcgggcggGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc            >  1:362152/1‑150 (MQ=255)
                                                                                                                                         gcgggcggGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa           <  1:653205/151‑1 (MQ=255)
                                                                                                                                             gcgggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg        >  2:11586/1‑150 (MQ=255)
                                                                                                                                               ggggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc      >  2:413958/1‑150 (MQ=255)
                                                                                                                                               ggggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc      >  2:921859/1‑150 (MQ=255)
                                                                                                                                                gggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     >  2:43943/1‑150 (MQ=255)
                                                                                                                                                  gtggtgATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  2:822202/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725

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