Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG43.3% 90.2 / 29.5 61V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (0/34);  new base G (26/0);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.43e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599439‑599726
                                                                                                                                         |                                                                                                                                                       
cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggggttctgggg                                                                                                                                            >  2:966090/1‑140 (MQ=255)
 tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                                 >  2:877757/1‑145 (MQ=255)
  ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGa                                                                                                                                            <  2:579225/149‑1 (MQ=255)
   tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGAT‑cc                                                                                                                                        >  2:209084/1‑151 (MQ=255)
   tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑c                                                                                                                                         <  1:295109/150‑1 (MQ=255)
    aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTTCGAT‑cc                                                                                                                                        >  2:651182/1‑150 (MQ=255)
    aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑cc                                                                                                                                        <  1:1027385/150‑1 (MQ=255)
    aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCg                                                                                                                                       <  2:125357/151‑1 (MQ=255)
     aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCgg                                                                                                                                      <  1:931759/151‑1 (MQ=255)
      tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGT‑CCgg                                                                                                                                      >  2:478256/1‑150 (MQ=255)
         caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCaa                                                                                                                                   <  1:252611/150‑1 (MQ=255)
             gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTTCCGTGAG‑CCGGCaaaa                                                                                                                                 >  1:725509/1‑148 (MQ=255)
             gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCTTGCGAT‑CCGGCaaaaa                                                                                                                                >  2:897029/1‑149 (MQ=255)
             gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACg                                                                                                                              <  1:217112/151‑1 (MQ=255)
              ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGAG‑CCGGCcaaaa                                                                                                                                >  2:897653/1‑148 (MQ=255)
                 tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCCGCGATCCCGGCaaaacaccc                                                                                                                            >  2:116864/1‑145 (MQ=255)
                  aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCga                                                                                                                         <  1:722385/151‑1 (MQ=255)
                  aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCga                                                                                                                         <  1:966324/151‑1 (MQ=255)
                    aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCgaga                                                                                                                       <  2:262461/151‑1 (MQ=255)
                    aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGCCTGCGAT‑CCGGCAAAAAAGCCCggg                                                                                                                        >  1:295313/1‑148 (MQ=255)
                         aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCCAAAAAGCCCGAGAAGTg                                                                                                                   >  1:995845/1‑150 (MQ=255)
                          gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                 <  2:771268/151‑1 (MQ=255)
                            gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAA‑CCGGCGAAAAAGCCCCAGAAgaggtc                                                                                                                >  1:868647/1‑145 (MQ=255)
                             aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGATCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCGCGAAGGGATc                                                                                                                >  2:397619/1‑149 (MQ=255)
                                 tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGTGGGGGGGGGGTCCTTTCAC‑ACGGCAAAAAAGCCCGGGAAgggggcaccccc                                                                                                          >  2:575274/1‑140 (MQ=255)
                                  ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                          <  2:50634/150‑1 (MQ=255)
                                  ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                         <  2:375148/151‑1 (MQ=255)
                                       gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                    <  2:705387/151‑1 (MQ=255)
                                        cccGTTTACCTGCGCGATGTTGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAACCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                     <  2:961059/149‑1 (MQ=255)
                                         ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCCGCAAAAACGCCCGAGAAGTGAGCGCCGCCGGGAAgg                                                                                                   >  2:684951/1‑150 (MQ=255)
                                           gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCCGCGAT‑CCGGCCAAAACGCCCCCGAAGGGATTGGCGCCGGGGAgg                                                                                                   >  2:210371/1‑148 (MQ=255)
                                           gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAAAGGCCGAGAAGGGATAGCCGCCGTGGAgg                                                                                                   >  1:619805/1‑148 (MQ=255)
                                              tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                              <  1:917993/150‑1 (MQ=255)
                                                   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                          <  1:227060/149‑1 (MQ=255)
                                                   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                           <  1:575158/148‑1 (MQ=255)
                                                    cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGAT‑CCGGCCAAAACGCCCCAGGAGGGATCGCCGCCGGGAAGGACAAAATggga                                                                                        >  2:1023373/1‑148 (MQ=255)
                                                            tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                                <  1:215749/150‑1 (MQ=255)
                                                               gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                            <  1:689038/151‑1 (MQ=255)
                                                                aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGGAAAGCTGaaaa                                                                              >  2:693695/1‑148 (MQ=255)
                                                                  ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCGGGGAGGGATCGCCGCCGTGAAGGGAAAAATGGAAAAGCTGAAAaata                                                                           >  2:315637/1‑149 (MQ=255)
                                                                   gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTGCTGCGAA‑CCCGCAAAAAAGCCCGGGGAGAGATCCGCGCCGGGGAGGGAAAACTGGAAAAGCCGGAaagaaggc                                                                        >  2:761737/1‑149 (MQ=255)
                                                                    tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAACCCGCGAAGGGAGAGCCGCCGGGAAGGGAAAACTGGAAAAGCTGAAAAGGAGTCt                                                                       >  2:628117/1‑151 (MQ=255)
                                                                     ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                       <  2:995733/150‑1 (MQ=255)
                                                                     ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                       <  1:252407/150‑1 (MQ=255)
                                                                      cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCCGGGGGGGGGTCCCGCGCA‑TCGGCAAAAACGCCCGCGGAGGGATCCCCGGCGCGGAGGGAAAACCAGGAAAGCTTAAAAGTAGTCt                                                                       >  2:981474/1‑149 (MQ=255)
                                                                         aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGTAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                  <  2:95128/151‑1 (MQ=255)
                                                                              cccGGAGATGCGCCGGGGCATTGCCGAACTAAAGGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAATCTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                             <  1:880667/151‑1 (MQ=255)
                                                                                      tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                      <  2:250828/150‑1 (MQ=255)
                                                                                        cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                    <  2:314005/150‑1 (MQ=255)
                                                                                          ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCCAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAggt                                                 >  1:93145/1‑149 (MQ=255)
                                                                                               gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                            <  1:376516/151‑1 (MQ=255)
                                                                                                cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGCGGGGATCCTTTGGT‑CCGGGCAAAAAGGCCGAGAAGAGGTCGCCGCCGCGAAGGACGAACTGGGAACCGTGAGAAGCAGACTGGCGGACGGCGTGgggaggggtaaaa                                             >  1:1035342/1‑140 (MQ=255)
                                                                                                 aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTacaaca                                           >  1:810203/1‑150 (MQ=255)
                                                                                                   tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGGGAAGGACAAACTTGAAAAGCCGAAAAGGAGTCTGCCGGAAGGGGTGGGGGTAGATtcaaca                                           >  2:506948/1‑148 (MQ=255)
                                                                                                    gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                        <  2:822221/150‑1 (MQ=255)
                                                                                                       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                    <  2:819597/151‑1 (MQ=255)
                                                                                                       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                    <  1:645639/151‑1 (MQ=255)
                                                                                                            aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                                <  2:428306/150‑1 (MQ=255)
                                                                                                                   gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                         <  2:295241/150‑1 (MQ=255)
                                                                                                                       gcgaagTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGAGAAACTGGGAACGCTGGAAAGTAGGCTGCCGGAAGGGGGGGGGATAGTTACAACATACGAACCCAGCCAGCTCATTg                     >  2:333479/1‑150 (MQ=255)
                                                                                                                          aagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg                 <  1:730295/151‑1 (MQ=255)
                                                                                                                                     ggggTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa      <  2:179162/151‑1 (MQ=255)
                                                                                                                                       ggtggtGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc    >  2:509320/1‑151 (MQ=255)
                                                                                                                                         tggtgATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  2:407995/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                       
CTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599439‑599726

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: ‑