Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,466,839 G→C 6.9% 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,8390GC 6.9% 197.4 / 7.5 59H137D (CAC→GAC) treRtrehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base C (2/2);  total (27/32)
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.98e-01

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

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: ‑