Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,014 C→T G100G (GGC→GGT nhaR → DNA‑binding transcriptional activator NhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,0140CT91.5% 156.1 / 2.2 59G100G (GGC→GGTnhaRDNA‑binding transcriptional activator NhaR
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (24/30);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGA  >  NC_000913/18779‑19259
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
ccGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACt                                                                                                                                                                                                                                         <  1:231598/250‑1 (MQ=255)
 cGTGGTGGGCGCAGCTGAGGCGCTTTATTTAACTCCACAAACCATTAACGGTCATATTCGAGCGCTGGAAGAGCGCCTGAAAGGCAAATTATTTAAACGCAATGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGAGATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACtgt                                                                                                                                                                                                                                       <  2:524266/251‑3 (MQ=255)
                  aGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGt                                                                                                                                                                                                                      <  1:209639/251‑1 (MQ=255)
                           atttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACtt                                                                                                                                                                                                             >  1:344211/1‑251 (MQ=255)
                             ttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTa                                                                                                                                                                                                            >  1:427201/1‑250 (MQ=255)
                             ttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTAGCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTa                                                                                                                                                                                                            <  2:57265/250‑1 (MQ=255)
                                  tCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAAcgccg                                                                                                                                                                                                      >  1:212593/1‑251 (MQ=255)
                                        aaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTgg                                                                                                                                                                                                <  2:339240/251‑1 (MQ=255)
                                        aaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTgg                                                                                                                                                                                                <  1:411004/251‑1 (MQ=255)
                                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGGATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGa                                                                                                                                                                                            <  2:664767/251‑1 (MQ=255)
                                            aTTACCGGACAGATTAGAGCGATGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCATGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCAGTTAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGGCGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGt                                                                                                                                                                                               >  2:174155/1‑248 (MQ=255)
                                                  ggACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGa                                                                                                                                                                                      <  2:65601/251‑1 (MQ=255)
                                                           cGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCa                                                                                                                                                                             <  2:264816/251‑1 (MQ=255)
                                                            gAGCGCTGGAAGAGAGACTGCAAGTCAAATTATTTAATCGAAAGGGACGTGGTCTCGAACCCAGCGAGCCGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTTTCGCAAAGAATCCAATTTATTGTTTGACGTGGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCAt                                                                                                                                                                            <  2:548801/251‑1 (MQ=255)
                                                                        agCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTc                                                                                                                                                                >  1:91865/1‑251 (MQ=255)
                                                                        agCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTCCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTc                                                                                                                                                                <  1:35594/251‑1 (MQ=255)
                                                                            ccTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAAt                                                                                                                                                            <  2:541618/251‑1 (MQ=255)
                                                                                cAAGGCAAATTATTTAAACGCAAGGGTCGTGGTTTCTAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCTCCTTAAGCCAGGAAATGCGGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGTCTTATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATccac                                                                                                                                                        <  2:663065/251‑1 (MQ=255)
                                                                                 aaGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGATCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGAAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGACCATTAATCTTCGCTGCTTCGAATccacc                                                                                                                                                       >  1:630658/1‑251 (MQ=255)
                                                                                    gCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGCCGTTGGTGTGGCTGATTCACTTTCCAAACGCCTGGTCAGTAGCTTACTTAACGCCGCAGTGGTAGAAGGCGATCCCATTCATCTTCGCGTCTTCGACTccacccac                                                                                                                                                    >  2:327657/1‑251 (MQ=255)
                                                                                     cAAATTATTTAAACGATATGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTAACGCTATGGCGAGAAAATGTTCACCTTAAGCCAGGTAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGAGGCACTTTCCAAACGCCGGGTCAGTAGCGTACTTAACGCCGCAGCGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACg                                                                                                                                                   <  2:51747/251‑1 (MQ=255)
                                                                                           attTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGGCAGTAGCGTACTTaa                                                                                                                                                                                                           >  1:312056/1‑189 (MQ=255)
                                                                                           attTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACGTTCCAAACGCCTGGTCAGTAGCGTACTTaa                                                                                                                                                                                                           <  2:312056/189‑1 (MQ=255)
                                                                                                   gCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAg                                                                                                                                     <  1:327657/251‑1 (MQ=255)
                                                                                                      aGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCaa                                                                                                                                  <  1:447338/251‑1 (MQ=255)
                                                                                                      aGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACGTGTTGATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCCCCCACGAAATGCTGCTGGAGCaa                                                                                                                                  <  2:273346/251‑1 (MQ=255)
                                                                                                            gTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGt                                                                                                                            >  1:333971/1‑251 (MQ=255)
                                                                                                              ggTCTCGAACCCAGCGAGCTGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAg                                                                                                                         <  2:11897/251‑1 (MQ=255)
                                                                                                                   cGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAAt                                                                                                                                               <  1:9299/225‑1 (MQ=255)
                                                                                                                        ccAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTc                                                                                                                                                                        >  1:187419/1‑195 (MQ=255)
                                                                                                                        ccAGCGAGCTGGGAGAACTGGGCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATGTATTGTTTGACGTTGGTGTGGCTGAAGCACTTGCCAATCGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTc                                                                                                                                                                        <  2:187419/195‑1 (MQ=255)
                                                                                                                              aGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACAATCGCAAAGCATCCAATTTATTGTTTGACGTTGGTGTGGCTGGTGCACTTTCCAAACGCCTGGTCAGAAGAGTACTTAACGCCGCAGGGGCAGAAGGCGAGCACATTCCTCTTCGCTGCTTCGAATCCACACACGAAAAGCTGCTGGAGCAATTAAGTCAGCATAATCTGGaaat                                                                                                           >  2:158967/1‑247 (MQ=255)
                                                                                                                                    gagaACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCAtt                                                                                                    >  1:474093/1‑251 (MQ=255)
                                                                                                                                          tGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGAc                                                                                              >  1:671879/1‑251 (MQ=255)
                                                                                                                                              cgaacGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCGGGATATTGTGAACTAGCGCAAAGAATCCAATTTAGTGTTTGACGTTGGTGTGGCTGATGCACGTTCCAAACGCCTGGTCAGTAGCGTGCTTAACGCCGCAGTGGTAGAAGGCGAGCCCAATCATCTTCGCTGCTTCGAAGCCACCCACGAAATGCTGCTGGTGCAATTAAGTCAGCATAAACTGTATATGAGCATGTCTGACTGt                                                                                           <  2:334792/246‑1 (MQ=255)
                                                                                                                                                aTCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAAACCAATTTATTGGTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTAAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAAGTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTcc                                                                                         <  2:344211/250‑1 (MQ=255)
                                                                                                                                                   gCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATa                                                                                     <  1:571947/251‑1 (MQ=255)
                                                                                                                                                         acGATAAAATGTTCACCATAAGCCAGGAAATGCTGGATATTGTGAACTATCGCACAGACTCCAATTTATTGTTTGACGTTGGTGTTGCTAATTCACTTTCCAAACGCCTGGTTCGTACCGTACTTACAGCCGCAGTGGTAGAAGCCGAGTCCTTTCAACTTCGTTGCTTCGATACCACCAACGAAATGCAGCACGAGCAAATAAGTCCACCTTAACTGGAAATGAACTTTTCTGAcatatcaatacaatcc                                                                               >  2:6127/2‑236 (MQ=255)
                                                                                                                                                            aTAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACg                                                                            >  1:445248/1‑251 (MQ=255)
                                                                                                                                                                       aCCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGGCGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAgg                                                                 <  1:174155/251‑1 (MQ=255)
                                                                                                                                                                                      aTGCTGGATATTGTGAACTATCGCAAAGAATCCATTTTATTGTTTGACGTTGTCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCATACTTAACGCCGCAGTGGTAGAATGCGAGCCCATTCATCTTCGCTGCTTCGAATCCCCCCACGAAATGCTGCTGGGGCAATTAAGTCCGCATAAACTGGATATGATCATTTCTGACTGTCCGAGAGACTCTCCGCAGCAGGAAGGCCTGTTCTCCGTg                                                    >  2:527825/1‑249 (MQ=255)
                                                                                                                                                                                                    gAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCACTTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGACGGCCTGTTCTCCGTGAGAATTGGCGAAtgtg                                    >  2:94412/1‑251 (MQ=255)
                                                                                                                                                                                                     aaCTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAATCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTgg                                   >  1:641019/1‑251 (MQ=255)
                                                                                                                                                                                                      aCTATCGCAAAGAATCCAATTTATTGTTTGTCGTTGGTGTGGCTGATGCACTTGCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGTTTAGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGc                                  <  2:14574/251‑1 (MQ=255)
                                                                                                                                                                                                         ttCGCAAAGAATCCAATTTATTTTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTTGTCAGTAGCGTACTTAACGCCTCAGTGGTAGAAGGCTAGCCCATGCATCTTCGCTGCTTCGAATCCGCCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACGCTACGCAGCAGGAAGGCCTGTTCTCCTTGAGAATTGGTGAATGTGTCGt                                <  2:427201/249‑1 (MQ=255)
                                                                                                                                                                                                           cGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgcagca                                                                       <  1:637582/209‑1 (MQ=255)
                                                                                                                                                                                                           cGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgcagca                                                                       >  2:637582/1‑209 (MQ=255)
                                                                                                                                                                                                           cGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGa                              <  2:390789/250‑1 (MQ=255)
                                                                                                                                                                                                             cAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGtt                           >  1:451673/1‑251 (MQ=255)
                                                                                                                                                                                                                     ccAATTTATTGTGTGACGTGGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCCCCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGgtgt                   <  1:224580/251‑1 (MQ=255)
                                                                                                                                                                                                                       aaTTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgcagc                                                                        <  1:586030/196‑1 (MQ=255)
                                                                                                                                                                                                                       aaTTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCAGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgcagc                                                                        >  2:586030/1‑196 (MQ=255)
                                                                                                                                                                                                                            aTTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCAttt                                                                                                   >  2:318345/1‑164 (MQ=255)
                                                                                                                                                                                                                            aTTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCAttt                                                                                                   <  1:318345/164‑1 (MQ=255)
                                                                                                                                                                                                                            aTTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATc            <  1:507859/251‑1 (MQ=255)
                                                                                                                                                                                                                                tttGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGGACAAATccacc        >  1:639336/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATccaccac      >  1:258984/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATccaccacc     >  1:372415/1‑248 (MQ=255)
                                                                                                                                                                                                                                       ttGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGGGTACAAATCCACCACCAGa  >  1:55690/1‑250 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
CCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGA  >  NC_000913/18779‑19259

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