Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 602,835 G→A A293T (GCT→ACT)  pheP → phenylalanine:H(+) symporter PheP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913602,8350GA98.4% 178.8 / ‑4.8 65A293T (GCT→ACT) phePphenylalanine:H(+) symporter PheP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/36);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

TTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTG  >  NC_000913/602595‑603080
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
ttggctatgtggggaaTTATGGTCTCGTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCTCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGGATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTTGCTCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTTTGATGATTTTCCATGATCTCGACAGCAACGTGGTAACTCCTgcgc                                                                                                                                                                                                                                              <  2:73193/236‑1 (MQ=255)
         gCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTGCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCGGAACTTCGt                                                                                                                                                                                                                                    <  1:313915/251‑1 (MQ=255)
          cGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTc                                                                                                                                                                                                                                   <  2:204325/251‑1 (MQ=255)
                             cggcggTCTGGATCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGt                                                                                                                                                                                                                <  2:468295/251‑1 (MQ=255)
                                  gTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGGGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAgtg                                                                                                                                                                                                           >  1:86470/1‑251 (MQ=255)
                                                  gaatACTGCCGATGAAGCGCGCGCTCCGGAGAAAAGCATTCCAAAAGCGGTAAATAAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGATTTACTGGAGCTCTATCCGTGGGTGGAAGTGAAATCCAACAACAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCgggggtt                                                                                                                                                                                           >  2:412750/4‑248 (MQ=255)
                                                   aTTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAAGCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTa                                                                                                                                                                                          <  1:414369/251‑1 (MQ=255)
                                                      acgggcGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAATCGGTAAATCAGGTGGTGGATCGCAGCCTGCTGTGTTACAGCGGTTCACTGGTGGTTTTACGGGCGGTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTGTGTGATGATTTTCCATAATCTCGACAGCAACGTGGGAACATCGGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGTTTTActc                                                                                                                                                                                       <  2:276221/246‑1 (MQ=255)
                                                                     cgcgATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCt                                                                                                                                                                        <  2:112655/251‑1 (MQ=255)
                                                                     cgcgATCCGGAAAAAAGCATGCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCt                                                                                                                                                                        <  2:448330/251‑1 (MQ=255)
                                                                      gcgATCCGGAAAAAAGCATTCCAAACGCGGTAAATCAGGTGGTGTATCGCATCATGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGAGGGGTTTACTCTAACAGCCGCATGCTg                                                                                                                                                                       >  2:564342/1‑251 (MQ=255)
                                                                      gcgAACCGGAAATCAGCATTCCAAAAGCGGTAAATACGGTGGTGGATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAGGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGACTAACAGCGGGGTTTACTCTAAAAGCCGcctggtg                                                                                                                                                                       >  2:54573/1‑245 (MQ=255)
                                                                       cgATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGt                                                                                                                                                                      <  1:332945/251‑1 (MQ=255)
                                                                         aTCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGttt                                                                                                                                                                    >  1:525191/1‑251 (MQ=255)
                                                                             ggAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGTTTCACTGGTTGTTTTACTGGCGTTCTATCCGTGGGTGGAAGTGAAATCCAACAGGAGCCCGTTTGTGATGATTTTCCATAATCGCTGCAGCAACGTGGTAACTTCTGCGCGGAACGTCTGCATTCTTGTAGCAGCTCTGTCAGTGGATAACAGCGGGCTTTACCCTAACGGCCGCATGCTGTGTGGcc                                                                                                                                                                <  2:294176/251‑1 (MQ=255)
                                                                             ggAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGcc                                                                                                                                                                >  1:475311/1‑251 (MQ=255)
                                                                                  aaatcATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTATCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCTTGCTGTTTGGCCTTTCt                                                                                                                                                           <  1:203502/247‑1 (MQ=255)
                                                                                  aaaGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCt                                                                                                                                                           >  1:406641/1‑251 (MQ=255)
                                                                                  aaaGCATTCCAAAAGCGGGAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTTACTGGTTGTTTGACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTTATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCGGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGGATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCt                                                                                                                                                           <  2:402219/251‑1 (MQ=255)
                                                                                       aTTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGc                                                                                                                                                       >  1:24827/1‑250 (MQ=255)
                                                                                             aaaGCGGTAAATCAGGTGGTGGATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCTCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTGTGGCCTTTCTGTGCAGGGTaa                                                                                                                                                <  1:444006/251‑1 (MQ=255)
                                                                                              aaGCGGTAAATCAGGTGGAGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGCAATCCAACAGTAGCCCGTTTGTGATGCTTTTCCATAATCTCGACAGCCACGTGGTACCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTact                                                                                                                                               >  2:122421/1‑249 (MQ=255)
                                                                                               aGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGGTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTTTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACATCGGGGTTTACGCTAGCAGCTGCATGCTGATGGGCCTTTCTGTGCAGGGTAATg                                                                                                                                              <  2:119263/251‑1 (MQ=255)
                                                                                                  ggTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATgcgc                                                                                                                                           >  2:465581/1‑251 (MQ=255)
                                                                                                  ggTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATgcgc                                                                                                                                           <  1:134982/251‑1 (MQ=255)
                                                                                                   gTAAATCAGGTGGTGTATCGCATCCTGCTTGTTTACCGCGGTTCACGGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTCAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACGTCTGCGCTGAACTTCGTCATTCGGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTAATCTAACAGCCGCAGGCTGTTTGGCCTGTCTGTGCAGGGTAATGCGcc                                                                                                                                          <  2:356293/251‑1 (MQ=255)
                                                                                                   gTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGcc                                                                                                                                          >  1:74419/1‑251 (MQ=255)
                                                                                                   gTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGcc                                                                                                                                          <  1:142587/251‑1 (MQ=255)
                                                                                                     aaaTCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTTTCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAATGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCg                                                                                                                                         <  2:179894/250‑1 (MQ=255)
                                                                                                         cAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGt                                                                                                                                     <  1:172761/250‑1 (MQ=255)
                                                                                                               gtgtATCGCATCCTGCTGTTTTCCATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGGAGTGAAATCCAACAGTAGCACGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGTCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAGGTTTTTGAc                                                                                                                              >  2:16924/1‑251 (MQ=255)
                                                                                                                 gtttCGCATCCTGCGGTTTGACAGCGGTTCACTGGGGGTTTTACTGGCGCTGTATCCGTGGGTGGAAGTGAAATCCAGCAGTATCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACGTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACGGCGGTGTTGACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCTAAGTTGTTGACt                                                                                                                             <  2:69551/247‑1 (MQ=255)
                                                                                                                               tgTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTcgc                                                                                                              >  1:69128/1‑251 (MQ=255)
                                                                                                                                 ttttACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCgg                                                                                                            >  1:192587/1‑251 (MQ=255)
                                                                                                                                   ttACATCGGATAACTGGTGGTTTTAATGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCACGTTTGTGATGATTTTCCATAATCTCGCCAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGATGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCAGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTAGCGTAAGCCGTCGCGgtg                                                                                                          >  2:34270/1‑251 (MQ=255)
                                                                                                                                          ggTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGAt                                                                                                   >  2:155618/1‑251 (MQ=255)
                                                                                                                                             tCACTGGTGGTTTTACTGGCTCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTGTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTa                                                                                                 <  1:83323/250‑1 (MQ=255)
                                                                                                                                                                   tttATCCGTGGGTGGAAGTGAAATCCAACTGTAGCCCTTTTGTGATGAGTTTCCATAATCTCGTCAGCAACGGGGTAACTTCTGCGCTGAACGTCGTCATTCTGGTAGCATCGCTTTCATTGTATAACTGCGTGGTGTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTTCAGGGTAATGCTCCGATGGTTTGGGCTCGCGTCAGCCGGCGCTGCGTGCCGATTAACTCTCTGATGCTTTCCGTAGc                                                                           <  2:289296/248‑1 (MQ=255)
                                                                                                                                                                   tctATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAAATTCGTCATTCTGGTAGCATCGATGTCAGTGTATAAAAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGAAGGGAAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCTGTGTGCCGATTAACTCTCTGATGCTTTCCGGAg                                                                            >  1:51151/1‑249 (MQ=255)
                                                                                                                                                                    ctATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCAGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCg                                                                          >  2:480648/1‑250 (MQ=255)
                                                                                                                                                                    ctATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGGTTGTGATGATGTTCCATAATCTCGACAGCAACTTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATTCTTTCCGGAGc                                                                           <  1:530774/249‑1 (MQ=255)
                                                                                                                                                                            gggtggAAGTGAAATCCTACAGTAGCCCGTTTGTGATGATGTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGGTGCACGTCGTCATTCTGGTAGCATGGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCGGTTTGGCCTTGCTGTGCAGGATAATGCGTCGAAGTTTTTGACGTGTGTCAGCCGTCGCGGTGTGCCAGTGAACTCGCTGATGCATTCTGGAGCGATCATTTc                                                                  <  2:154276/250‑1 (MQ=255)
                                                                                                                                                                            gggtggAAGTGAAATCCAACAGTAGTCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGGTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTc                                                                  <  1:32128/250‑1 (MQ=255)
                                                                                                                                                                                   aGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCtggtgg                                                          <  1:241357/251‑1 (MQ=255)
                                                                                                                                                                                    gTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTggtggt                                                         <  1:518160/251‑1 (MQ=255)
                                                                                                                                                                                       aaaTCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGtt                                                      <  1:511119/251‑1 (MQ=255)
                                                                                                                                                                                             aaCAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATGCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCaa                                                <  1:55636/251‑1 (MQ=255)
                                                                                                                                                                                                   aGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATct                                          >  1:390111/1‑251 (MQ=255)
                                                                                                                                                                                                    gCCCGTTTGAGATGATTTTCCATAATCTCGCCAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTCGCATCGCTGTCAGTGTATAACAGCGGGGATTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTCATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctg                                         >  1:596792/1‑251 (MQ=255)
                                                                                                                                                                                                      ccGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctgct                                       <  1:564342/251‑1 (MQ=255)
                                                                                                                                                                                                        gTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCtgctgc                                     >  2:1114/1‑251 (MQ=255)
                                                                                                                                                                                                           tgtgATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCg                                   <  1:181906/250‑1 (MQ=255)
                                                                                                                                                                                                            gtgATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATGCTGGTAGCATCGCTGGCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATTCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAAGCAACGATCTGCTGCCGCa                                 <  2:252972/251‑1 (MQ=255)
                                                                                                                                                                                                                    tttCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTAATACTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCAGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAACGCGt                         >  1:47982/1‑251 (MQ=255)
                                                                                                                                                                                                                     ttCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGtt                        >  1:520017/1‑251 (MQ=255)
                                                                                                                                                                                                                          tAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTc                   >  1:94830/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCTCGACAACCACGTGGTAACTTATGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATCACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGCTTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAACAGCGTTTGGTctg                 >  2:45730/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTTCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTctg                 <  2:74419/247‑1 (MQ=255)
                                                                                                                                                                                                                                  aCAGCCACGTGGTAGCTTCTGCGCTGAAACTCGTCCTTCTGGTAGCATCGCTGTCAGTGAATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATACGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATg           >  2:321707/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATgg          <  2:69128/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGgc         <  1:128557/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCGTGGTAACTTCTGCGCGGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGGCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCt      <  1:542910/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCtggt   >  2:45055/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCtggtg  >  1:572020/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggTAACTTCTGCGCTGAACTTCGTCATTCGGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCtggtg  <  1:407360/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
TTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTG  >  NC_000913/602595‑603080

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