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

CTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599431‑599722
                                                                                                                                                 |                                                                                                                                                   
cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                 >  2:994233/1‑148 (MQ=255)
cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                 <  1:866447/150‑1 (MQ=255)
  cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGAt                                                                                                                                               <  1:79706/150‑1 (MQ=255)
    acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg                                                                                                                                                 >  1:376182/1‑144 (MQ=255)
    acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                            <  2:362981/151‑1 (MQ=255)
    acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                             >  1:940417/1‑150 (MQ=255)
       acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                             >  2:488599/1‑147 (MQ=255)
             aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg                                                                                                                                   <  1:1022421/151‑1 (MQ=255)
               cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGACCgg                                                                                                                                   >  2:340024/1‑149 (MQ=255)
                  aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGG‑Caaa                                                                                                                              >  2:330841/1‑150 (MQ=255)
                    cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAAc                                                                                                                           <  2:548143/151‑1 (MQ=255)
                       tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAACACTTGATCCCG‑CAAAAACGcc                                                                                                                        >  1:394093/1‑151 (MQ=255)
                        ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGACCGGGATCCGG‑Caaaaa                                                                                                                            >  1:1162031/1‑146 (MQ=255)
                         tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAACCGG‑CCAAAAAGccc                                                                                                                       >  2:725424/1‑150 (MQ=255)
                          aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCCGG‑CCAAAAAGccc                                                                                                                       >  2:747411/1‑149 (MQ=255)
                               ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGt                                                                                                                <  2:835284/151‑1 (MQ=255)
                               ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGG‑Caaaaaaacc                                                                                                                        >  1:578129/1‑139 (MQ=255)
                                 aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGa                                                                                                              <  1:82686/151‑1 (MQ=255)
                                     aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGG‑GGGGATCCCGCGATCCGG‑CAAAAAAGCCCGAGAAGGGATcgcc                                                                                                         >  2:666003/1‑151 (MQ=255)
                                            cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCCCGAAGTGATCGCCGCCGTg                                                                                                   >  2:617220/1‑151 (MQ=255)
                                               gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                   >  2:795450/1‑148 (MQ=255)
                                                cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                               <  1:867354/151‑1 (MQ=255)
                                                   gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                            <  2:229707/151‑1 (MQ=255)
                                                      tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                          <  1:765496/150‑1 (MQ=255)
                                                      tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                          <  2:221028/150‑1 (MQ=255)
                                                       aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAAAACCCGAGAAGTGATCGCCGCCGGGAAGGACaaa                                                                                          >  2:603906/1‑149 (MQ=255)
                                                        tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                       >  2:307325/1‑151 (MQ=255)
                                                         cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCCAGAAGAGATCGCCGCCGTGGAGGACAAACTgg                                                                                      >  2:32687/1‑151 (MQ=255)
                                                         cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGG‑CAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTgg                                                                                      >  1:1024058/1‑151 (MQ=255)
                                                          tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCCTCGATCCGG‑CAAAAAAGCCCGCGAAGTGGTCGCCGCCGTGGAGGGCaaaacgg                                                                                      >  1:159378/1‑146 (MQ=255)
                                                           gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                        <  2:598502/147‑1 (MQ=255)
                                                           gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                        <  1:941794/147‑1 (MQ=255)
                                                           gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                      <  1:1037218/149‑1 (MQ=255)
                                                           gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                       <  2:242415/148‑1 (MQ=255)
                                                              cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGG‑CCAAAAACCCCGGGAAGGGATCGCCGCCGTGCAGGGCAAACTTGaaa                                                                                   >  2:80333/1‑149 (MQ=255)
                                                                aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGGTCCGG‑CAAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAACTGGaaa                                                                                   >  2:469679/1‑147 (MQ=255)
                                                                 tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGG‑CCAAAACGCCCGAGAAGAGATCGCCGCCGGGGAAGACAAACTCGAAACGCt                                                                               >  1:708897/1‑150 (MQ=255)
                                                                  gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGG‑CAAAAAAGCCCGGGGAGGGGTAGCCGGCGGGGAGGGCGAACTGGGAACGCTg                                                                              >  1:790988/1‑150 (MQ=255)
                                                                      cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGCGGTCCCG‑CAAAAACGCCCGAGAAGTGATCGCCGCCCCGGAGGGCAAACTGGAAACACCGaaaa                                                                          >  1:295508/1‑150 (MQ=255)
                                                                       gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                        <  1:1081958/151‑1 (MQ=255)
                                                                        aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGTGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                       <  2:788553/151‑1 (MQ=255)
                                                                        aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCCATCCGGCCAAAAACG‑CCGGGAAGTGATCGCCGCCGTGGGGGGCAAACTGGGAACGCTGaaa                                                                           >  1:406119/1‑147 (MQ=255)
                                                                         aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                      <  2:14205/151‑1 (MQ=255)
                                                                           gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGG‑GAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGGAAaggagac                                                                    >  1:1079274/1‑149 (MQ=255)
                                                                             ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGCTCCCG‑CAAAAAACCCCGAGAAGTGGTTGCCGCCGCGGAGGACAAACTGGAAACGCTGAAAaatagattg                                                                  >  1:416805/1‑147 (MQ=255)
                                                                                gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCTTGCATCCCG‑CGAAAACGACCGGCGAGGGCACGGCGCCGGGTAGGACCAACTGGGAACGTTggcaagaaggtggccg                                                               >  1:617657/1‑136 (MQ=255)
                                                                                 aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                              <  2:389361/151‑1 (MQ=255)
                                                                                 aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                               <  1:273725/150‑1 (MQ=255)
                                                                                     gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  1:674221/151‑1 (MQ=255)
                                                                                      cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                         <  1:934250/151‑1 (MQ=255)
                                                                                        cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                       <  2:596567/151‑1 (MQ=255)
                                                                                         ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  1:801875/151‑1 (MQ=255)
                                                                                         ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  1:722575/151‑1 (MQ=255)
                                                                                          gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                     <  1:787501/151‑1 (MQ=255)
                                                                                           agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGCGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                    <  2:693747/151‑1 (MQ=255)
                                                                                              tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCCGGAAAGTAGTCTGCCGGAAGGCgggggg                                                   >  2:305791/1‑147 (MQ=255)
                                                                                                 gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAAAGCCCGAGAAGTGATCGCCGCCGCGAAGGACAAACTGGAAACGCTGGAAAGGAGTCTGCCGGGAGGCGTGGGGATa                                                >  1:991232/1‑149 (MQ=255)
                                                                                                    ggggCATGGCCGAACTAAGCGGGGAAGGCGAGTTGGGGGGCGGGGGGGTGATCCTGCGGTCCGG‑AAAAACCCCCCGAGAAGTGAGCACCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                            <  1:557966/150‑1 (MQ=255)
                                                                                                           tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                     <  2:666848/150‑1 (MQ=255)
                                                                                                              cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                 <  1:1087981/151‑1 (MQ=255)
                                                                                                               gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  2:828597/151‑1 (MQ=255)
                                                                                                                    aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                             <  1:100903/149‑1 (MQ=255)
                                                                                                                     aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGATCCGG‑CAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGAAAGACAAAATGGAAACGCTGAAAAGGAGTCTGGCGGGAGGCGGGGGGATAGTTACAACATACGATCGca                             >  2:207800/1‑148 (MQ=255)
                                                                                                                        ggcgaaggcgaaGTGGCGGGCGGGGGGGTGATCCTGCGATCCGG‑CAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTAAAACATACGATCGcagccag                        >  2:363570/1‑150 (MQ=255)
                                                                                                                             aggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                  <  2:955747/151‑1 (MQ=255)
                                                                                                                              ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                 <  2:632650/151‑1 (MQ=255)
                                                                                                                                 gaagTGGCGGGCGGGGGGGTGGTCCTGCGCTCCCG‑CAAAAAAGCCCCGGGAGTGGTTGCCGCCGTGGAGGGCAAAATGGGAAAGCTGAAAAGTAGGTTGCCGGGAGGGGGGGGGGTAGTTACAACATACGATCGCAGCCCCCTCATTGTcc              >  1:114662/1‑151 (MQ=255)
                                                                                                                                         gggcgggGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc      <  2:489288/151‑1 (MQ=255)
                                                                                                                                           gcgggGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGa    <  1:523071/151‑1 (MQ=255)
                                                                                                                                             ggggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   <  2:209405/150‑1 (MQ=255)
                                                                                                                                             ggggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   >  1:525714/1‑150 (MQ=255)
                                                                                                                                             ggggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  <  1:756697/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599431‑599722

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