Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG27.7% 68.6 / 10.3 47V288G (GTG→GGG) cusAcopper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base T (1/33);  new base G (12/1);  total (13/34)
Fisher's exact test for biased strand distribution p-value = 3.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCG  >  NC_000913/599327‑599826
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                              
aaGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                                                                                                                               <  1:671981/251‑1 (MQ=255)
      aaaGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGGTCGTCGATCGAACGGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                                                                                                                                                                                                                                                         <  2:45045/251‑1 (MQ=255)
         gcgcgcTGGATGCGTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAAGATATGGTGCGCGCCAGCGGCTATCTGCAAACGCGCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTAGCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTg                                                                                                                                                                                                                                                      <  2:71308/251‑1 (MQ=255)
          cgcgcTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                                                                                                                                     <  1:150239/251‑1 (MQ=255)
                 gATGCGTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATGGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑gg                                                                                                                                                                                                                                             <  2:337665/251‑1 (MQ=255)
                                       gcggTTCGTCGATCGAACTGGCGGAAGCGGAAAATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATc                                                                                                                                                                                                                      <  2:641438/251‑1 (MQ=255)
                                          gTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATcgcc                                                                                                                                                                                                                   <  1:439946/251‑1 (MQ=255)
                                             cgtcgATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCTATCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                                                                                                                                <  2:92673/251‑1 (MQ=255)
                                                              gaagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑a                                                                                                                                                                                             <  1:367201/251‑1 (MQ=255)
                                                               aagcggaaTATATTGTGCGCGCCGGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCAGTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑aa                                                                                                                                                                                            <  2:350342/251‑1 (MQ=255)
                                                               aagcggaaTATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑aa                                                                                                                                                                                            <  1:45924/251‑1 (MQ=255)
                                                                ggcggaaTATAGGGTGCGCGCAAGCGGCTATCTACAAACCCTCGACGACTTTAATCACATCGTTTGAAAAGCCAGTGAAAAGGGCGTGCCCGTTTATCTGCGCGATGTGGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑aaa                                                                                                                                                                                           <  2:612834/250‑1 (MQ=255)
                                                                        atatGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGaa                                                                                                                                                                                    <  1:132244/250‑1 (MQ=255)
                                                                            ggTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTGGCGAATGTCCAGATTGGCCCGGAGATGCGCCGGGGCAGTGCCGAACTAAACGGCGAGGGCGAGGTGACGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCGCGATAAGCTATCTCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAagt                                                                                                                                                                                <  1:398174/250‑1 (MQ=255)
                                                                            ggTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGACCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTGCTGCGATCC‑GG‑CAAAAACGCCCCAGGAGTGACTGGCGCCGGGGAGGACAAACCG‑G‑GAACGCTGTAAagg                                                                                                                                                                                >  1:67659/1‑249 (MQ=255)
                                                                               gggcgccagggtcaaacTGCAAACGCTCGCCTACTTTAATCGCATCGTTTTAAAAGCCAGTAAAAAGGGCGTGCCCGTTTATCTGCGTGATGTTGCGAAGGTCCAGATTGGCCCGTAGATGCGCAGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGTGGGGCGGGGTGGTGCTCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTc                                                                                                                                                                            <  2:184326/235‑1 (MQ=255)
                                                                                                   aaaCGCTCGACGCCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                                                                                                                        <  1:521034/251‑1 (MQ=255)
                                                                                                      cGCTCGACGACTTTAATCACATCGTTTTACAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAGCC‑GG‑CCACAACGCCCGGGAAGGGATCGGCGCCCTGGAGGGGAAACTGTG‑AAA‑GCTGGAAAATAGTCTggcgcggaggggggggggaaat                                                                                                                                                     >  2:667081/1‑230 (MQ=255)
                                                                                                           gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGG‑CAAAAA‑GCCCGGGAAGGGATAGGCGCCGGGGAGGGCAAAATG‑G‑GAACGCTGAAAAGTAGGCTGCCGGGCGGCGGGGGGATAGTTtcaac                                                                                                                                                >  1:71665/1‑251 (MQ=255)
                                                                                                             cgccttgactcacttcgctgtAAAAGCCAGTGAAAATGGCGTGCCCGTTTACCGGCGCCATGGTGCGAAGGTCCAGATGGGCGCGGAGAGGCGCCGGGGTATTGCGGAAGTAAACGGGGAAGGCGAAGTCGCGGCCGTGTTGCTGATCCTACGCTCC‑TG‑CAGAAAGGCCCTAGAAGTGATCGCCGCCGTGAAGGACACGCTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAATGCGTGTAGGTAGTTCCAACAt                                                                                                                                              <  2:672507/231‑1 (MQ=255)
                                                                                                                    aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                                                                                                                        <  1:590646/250‑1 (MQ=255)
                                                                                                                     aTCACAGCGTTTTAAAAGCCAGGGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATGGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAATTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGGGTGGAGATAGTTACAGCATACGATCGc                                                                                                                                      <  2:660695/251‑1 (MQ=255)
                                                                                                                        acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GG‑CCAAAACGCCCGAGAAGTGGTCGCCGCCGCGAAGGGAAAACAG‑G‑AAAAGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGGTACAACATACGAACGccac                                                                                                                                   >  1:604277/1‑248 (MQ=255)
                                                                                                                            cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                                                                                                                                <  1:632621/250‑1 (MQ=255)
                                                                                                                                ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                                                                                                                            <  1:489507/250‑1 (MQ=255)
                                                                                                                                   aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGGACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATTCCGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGGAAAGGAATCTGGCGGGAGGCGGGGGGAGAGTTACAACATACGATCGcagacagccccttg                                                                                                                         >  1:333343/1‑244 (MQ=255)
                                                                                                                                     aGCCAGTGAAAGTGGCGTGCCCGTTTATCTGCGCGATGTGGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATGGCCGAACGAAACGGCGAAGGCGAAGTGGCGTGCGGGGTGGGGATCCTGCGATCC‑GT‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAGCGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                                                                                                                      <  2:329834/251‑1 (MQ=255)
                                                                                                                                      gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCC‑GG‑CAAAAACGCCCGCGAAGTGATCGCCGCCGGGCAGGACAAACCG‑G‑AAACGCGGAAAAGTAGTCTGCCGGAAGGCGTGGAGGTAGTTTCAAAAAAAGAACGAAGCCAGCTccattaccc                                                                                                                     >  1:187542/1‑243 (MQ=255)
                                                                                                                                       ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc                                                                                                                    <  1:323396/251‑1 (MQ=255)
                                                                                                                                          gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc                                                                                                                 <  1:106718/251‑1 (MQ=255)
                                                                                                                                              aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGCGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGGTATGCGATCC‑GG‑CCAAAACGCCCGAGAAGTGAACGCCGCCGTGAAGGACAAACCG‑G‑GAACGCTGGAAAGTATTCTGCCGGCAGGGGGGGAGATAGTTACAAAATACGATCCCCGCCCGCTCAGTGACCGCGCCATCg                                                                                                             >  1:329341/1‑251 (MQ=255)
                                                                                                                                                  ggCGTGCCCGTTTATCTGCGCTATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGTCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACTCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa                                                                                                         <  1:643297/251‑1 (MQ=255)
                                                                                                                                                                  ggggcgagggtGCGAAGGTCGCGATTGGCCCGGAGAGGCGCCGGGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACTCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCGGAAAAGTAGTCGGCCGGAAGGCGTGGAGATAGTTACAACATACTATCGCAGCCAGGTCATTGACCGCGGCATCGACAACCGCAGCGGCAAGTTg                                                                                         <  2:315800/241‑1 (MQ=255)
                                                                                                                                                                   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAg                                                                                            <  1:459876/247‑1 (MQ=255)
                                                                                                                                                                             gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGt                                                                               <  1:140169/250‑1 (MQ=255)
                                                                                                                                                                                    tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGCTCC‑GGGCAAAAA‑ACCCCCGAAGGGGTCGCCGCCGCGGAGGACAAACTG‑G‑GAAAGCCGGAAAGGAGTCTGCCGGGAGGCGGGGGGATAGTTAAAACATACGATCGCAGCCAGCTCATTGTCCGCGCCATCGCCAACCCCAGCGGCCAGCTGCTGGTAGAGGTTTTtgtg                                                                       >  2:500856/1‑251 (MQ=255)
                                                                                                                                                                                     ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGTCC‑GG‑CCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGGAGGCTGCCGGGAGGCGTGGGGATAGGTACAACATACGGTCGCTGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGGTTATTgtgg                                                                      >  1:496450/1‑251 (MQ=255)
                                                                                                                                                                                                cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGTTTAGTGTGGTGGCGGTAGTc                                                           <  2:67659/251‑1 (MQ=255)
                                                                                                                                                                                                     gtGCGCCGGGGGTGTGCCGAATGAAACGGAGAAGGGGAAGGGGCGGGAGGGGGGGTGAGCCGGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGGTGCTGGAAGAGTTTAGTGTGGTGGCGGTAGTCTGTgc                                                      <  2:218202/250‑1 (MQ=255)
                                                                                                                                                                                                           cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGCTCCTGCGCTCC‑GG‑CAAAAACGCCCGAGAAGAGATCGCCGGCGCGAAGGGCAAACTG‑GGAAAC‑CCGAAAAGTAGTCTGGCGGGAGGGGTGGGGATAGTTACAACATACGGTCCCAACCAGCTCATTGAACGCGCCACCCACAACCTCCGCGGCAAGCTGCTGGTAGAGATTATTGTGGTGGCGGGAGGCtgtggggtgtt                                                >  1:157783/1‑244 (MQ=255)
                                                                                                                                                                                                            ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCC‑GG‑CCAAAACGCCCGAGCAGTGATTGCCGCCGTGAAGGACCAACTT‑G‑GAACGCTGTAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCCTTTACCCCGCCCTCGACAACCTCAGCGGCAAGTTTCTGGGAGAGTTTATTGTGGGGGCGGGAGGCTGTGCGCTGttt                                               >  1:486049/1‑251 (MQ=255)
                                                                                                                                                                                                                     ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCAt                                      <  1:436379/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGAGATCC‑GG‑CAAAAACGCC‑GAGGAGGGATCGCCGCCGTGAAGGACAAACTG‑G‑GAAAGCTGAAAAGTAGGCTGCCGGGAGGGGGGGGGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGGAGAGTTTTTTGTGGGGGCGGGAGGCTGTGCGCTGGTTTTCTgggctgtggggtc                              >  1:649108/1‑240 (MQ=255)
                                                                                                                                                                                                                                         gaagTGTCGGGCGGGGTGGTGATCGTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCGGTGCGCTGTTTCTCTGGCATGTGCGCTCGTCGCTGGTGGc                  <  1:437496/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gggggggTGGTGATCCGGCGGTCC‑GG‑GAAAAACGCCCGAGAAGTGATCGACGCCGAGAAGGACAAACTG‑G‑AAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTc         <  2:379937/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgggGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGt       <  1:136210/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggggTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGCAGTGATCACCGCCGTGAAGGACACACTG‑G‑AAACGCTGACACGTCGTCTGCCGGAAGCCATGGAGATAGTTACAACATACGATCCCAGCTAGCCCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt      >  1:316803/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     ggggTGGGGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGtt      <  2:491730/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        gtggtgATCCTGAGATCC‑GG‑CAAAAACGCCCGAGAAGTGATAGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAACAGTAGTCTGCCGGCAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAACTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGCCTGTGCGCTGTTTCTCTGGCCTGTGCGCTCGGAGCTGGTGGCGATTATTTCGTTGcc   >  1:324821/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         tggtgCTCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGCTCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAACAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTTGCGATTATTTCGTTGccc  >  2:536874/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                              
AAGTAAAAAGCGCGCTGGATGCTTCAAACCAGGAAGCGGGCGGTTCGTCGATCGAACTGGCGGAAGCGGAATATATGGTGCGCGCCAGCGGCTATCTGCAAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTG‑G‑AAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCG  >  NC_000913/599327‑599826

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