Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,466,838 G→C 7.0% H137D (CAC→GAC)  treR ← trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,466,8380GC 7.0% 194.0 / 7.7 58H137D (CAC→GAC) treRtrehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base C (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

AGGCGTAATCACTTTTGCAACGTTCTCATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCA  >  NC_000913/4466554‑4467122
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                             
aGGCGGAATCACGTTTGCAACGTTCTCATAGCCTTGCTTCATGGCAAGCCCGGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCt                                                                                                                                                                                                                                                                                           <  2:103955/289‑1 (MQ=255)
   cGTAATCACGTTTGCAACGTTCTCATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCt                                                                                                                                                                                                                                                                                           <  1:30475/286‑1 (MQ=255)
   cGTAATACCTTTTGCAACGTTCTCATGTCCTTGCGTCATAGCAAGCCCTGGGAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCt                                                                                                                                                                                                                                                                                           <  1:460113/286‑1 (MQ=255)
                          cATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTC‑TTCTGTTATGCCAGTaaa                                                                                                                                                                                                                                                                >  1:490487/1‑289 (MQ=255)
                               ccTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCg                                                                                                                                                                                                                                                            <  2:289263/288‑1 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCTGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCCTACAGCCGTTGCAGCAGGATTTTGATTGCACCTTCGTCGTCATAACAGACAGAACCAAAGCCTTTTGCGTCACGCGCCAGAAGGACCACCGTTGACTGCCAGT‑GGCTAAACTATCATTCTGTTCTGCCAGTcaaccggacct                                                                                                                                                                                                                                                        >  2:230455/1‑280 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTAt                                                                                                                                                                                                                                                                         >  1:402451/1‑271 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTAt                                                                                                                                                                                                                                                                         <  2:402451/271‑1 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACa                                                                                                                                                                                                                                                        >  1:338032/1‑288 (MQ=255)
                                              aGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGt                                                                                                                                                                                                                                            <  2:490487/289‑1 (MQ=255)
                                                    ggCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACCACCCTAAGCCGTCGatat                                                                                                                                                                                                                                      >  2:44465/1‑289 (MQ=255)
                                                          gcggcAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGc                                                                                                                                                                                                                                 <  2:338032/288‑1 (MQ=255)
                                                                    gATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  1:44465/288‑1 (MQ=255)
                                                                                  gcTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGaa                                                                                                                                                                                                                                                          >  2:162593/1‑239 (MQ=255)
                                                                                  gcTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGaa                                                                                                                                                                                                                                                          <  1:162593/239‑1 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  1:129239/271‑1 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       >  2:129239/1‑271 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTTGATATTACGCCGTTTCAGCACCCCCAAATTTTCGGCa                                                                                                                                                                                                      >  2:151614/1‑288 (MQ=255)
                                                                                                cAGGTAAGCTGCGTGACGTCGCTTACCGGGTGCCACGTCACTGTGCGGCACGTCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCACTTCGTCGTCATAACAGACCGAAGCAAAGCCATTTGCGTAACGCGCCAGCAGAACCAGCGATGAATGACAATGAGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCgg                                                                                                                                                                                           <  1:230455/288‑1 (MQ=255)
                                                                                                       gCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAAc                                                                                                                                                                                    <  2:8451/288‑1 (MQ=255)
                                                                                                         ttCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCCGTAAAACCGAACAGCACTACGCCGTCGATATTACCCCGTTTCAGCACTCCCAAATGTTCGGCCACTAATTGCGGGGAAAACTg                                                                                                                                                                                  >  2:153971/1‑288 (MQ=255)
                                                                                                                 gTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         <  2:330631/289‑1 (MQ=255)
                                                                                                                       tACCGGTTGTCACGTCCCTGTGCGGCACGCCGATATAACTGCTATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   <  2:134510/289‑1 (MQ=255)
                                                                                                                                   cgtcacTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAAcc                                                                                                                                                       <  2:456470/289‑1 (MQ=255)
                                                                                                                                      cacTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCtt                                                                                                                                                     <  1:151614/288‑1 (MQ=255)
                                                                                                                                                  cGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGtt                                                                                                                                                                                                           <  2:465569/222‑1 (MQ=255)
                                                                                                                                                  cGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGtt                                                                                                                                                                                                           >  1:465569/1‑222 (MQ=255)
                                                                                                                                                                   ttACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCg                                                                                                                                                                                                         <  1:135831/207‑1 (MQ=255)
                                                                                                                                                                   ttACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCg                                                                                                                                                                                                         >  2:135831/1‑207 (MQ=255)
                                                                                                                                                                        aTGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         >  2:271792/1‑234 (MQ=255)
                                                                                                                                                                        aTGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         <  1:271792/234‑1 (MQ=255)
                                                                                                                                                                          gCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCgag                                                                                                                >  2:274676/1‑289 (MQ=255)
                                                                                                                                                                          gCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCgag                                                                                                                >  2:308291/1‑289 (MQ=255)
                                                                                                                                                                                ggTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGAttt                                                                                                            >  2:47669/1‑287 (MQ=255)
                                                                                                                                                                                      aCAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   <  2:13202/226‑1 (MQ=255)
                                                                                                                                                                                      aCAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   >  1:13202/1‑226 (MQ=255)
                                                                                                                                                                                          ccGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAACCTGCCTTTCCATCATGCTTGGTTCGTACCCTTGTTCATAGAACGCTGCCAGCATGGTTTGAACGGCGAGATTTTTTGCCAACGa                                                                                                 >  2:172311/1‑288 (MQ=255)
                                                                                                                                                                                          ccGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGaa                                                                                                >  1:321075/1‑289 (MQ=255)
                                                                                                                                                                                           cGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGCTTGGGTCGTAACCTTTTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGAtttttcgaacacgca                                                                                                >  2:20859/1‑277 (MQ=255)
                                                                                                                                                                                            gTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAAt                                                                                               >  1:51747/1‑288 (MQ=255)
                                                                                                                                                                                                       aTTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGt                                                                                     <  1:308291/287‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACaa                                                                                                    <  1:293672/263‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACaa                                                                                                    >  2:293672/1‑263 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAg                                                                                           >  1:433058/1‑272 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAg                                                                                           <  2:433058/272‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAatgatg                                                                           >  1:126219/1‑288 (MQ=255)
                                                                                                                                                                                                                 ccccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  1:64487/147‑1 (MQ=255)
                                                                                                                                                                                                                 ccccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       >  2:64487/1‑147 (MQ=255)
                                                                                                                                                                                                                  cccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCt                                                                                                                                       <  1:47233/226‑1 (MQ=255)
                                                                                                                                                                                                                  cccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCt                                                                                                                                       >  2:47233/1‑226 (MQ=255)
                                                                                                                                                                                                                    tttcttcGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGAGGACTGCCAGTGTGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGAc                                                                      <  1:427718/288‑1 (MQ=255)
                                                                                                                                                                                                                                 aCAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCtt                                                         <  2:51747/289‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGccc                                                    >  2:99036/1‑288 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTAt                                                              >  1:199580/1‑256 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTAt                                                              <  2:199580/256‑1 (MQ=255)
                                                                                                                                                                                                                                                                      agAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgggg                       <  1:47669/286‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGc            <  1:172311/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        cAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATtcatca  <  1:274676/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                             
AGGCGTAATCACTTTTGCAACGTTCTCATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTC‑TTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCA  >  NC_000913/4466554‑4467122

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