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

AGGCGTAATCACTTTTGCAACGTTCTCATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCA  >  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)
                          cATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTCTTCTGTTATGCCAGTaaa                                                                                                                                                                                                                                                                >  2:490487/1‑289 (MQ=255)
                               ccTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCg                                                                                                                                                                                                                                                            <  1:289263/288‑1 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCTGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCCTACAGCCGTTGCAGCAGGATTTTGATTGCACCTTCGTCGTCATAACAGACAGAACCAAAGCCTTTTGCGTCACGCGCCAGAAGGACCACCGTTGACTGCCAGTGGCTAAACTATCATTCTGTTCTGCCAGTcaaccggacct                                                                                                                                                                                                                                                        >  1:230455/1‑280 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTAt                                                                                                                                                                                                                                                                         >  2:402451/1‑271 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTAt                                                                                                                                                                                                                                                                         <  1:402451/271‑1 (MQ=255)
                                   gCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACa                                                                                                                                                                                                                                                        >  2:338032/1‑288 (MQ=255)
                                              aGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGt                                                                                                                                                                                                                                            <  1:490487/289‑1 (MQ=255)
                                                    ggCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACCACCCTAAGCCGTCGatat                                                                                                                                                                                                                                      >  1:44465/1‑289 (MQ=255)
                                                          gcggcAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGc                                                                                                                                                                                                                                 <  1:338032/288‑1 (MQ=255)
                                                                    gATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  2:44465/288‑1 (MQ=255)
                                                                                  gcTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGaa                                                                                                                                                                                                                                                          >  1:162593/1‑239 (MQ=255)
                                                                                  gcTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGaa                                                                                                                                                                                                                                                          <  2:162593/239‑1 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  2:129239/271‑1 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       >  1:129239/1‑271 (MQ=255)
                                                                                     ttGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTTGATATTACGCCGTTTCAGCACCCCCAAATTTTCGGCa                                                                                                                                                                                                      >  1:151614/1‑288 (MQ=255)
                                                                                                cAGGTAAGCTGCGTGACGTCGCTTACCGGGTGCCACGTCACTGTGCGGCACGTCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCACTTCGTCGTCATAACAGACCGAAGCAAAGCCATTTGCGTAACGCGCCAGCAGAACCAGCGATGAATGACAATGAGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCgg                                                                                                                                                                                           <  2:230455/288‑1 (MQ=255)
                                                                                                       gCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAAc                                                                                                                                                                                    <  1:8451/288‑1 (MQ=255)
                                                                                                         ttCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCCGTAAAACCGAACAGCACTACGCCGTCGATATTACCCCGTTTCAGCACTCCCAAATGTTCGGCCACTAATTGCGGGGAAAACTg                                                                                                                                                                                  >  1:153971/1‑288 (MQ=255)
                                                                                                                 gTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         <  1:330631/289‑1 (MQ=255)
                                                                                                                       tACCGGTTGTCACGTCCCTGTGCGGCACGCCGATATAACTGCTATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   <  1:134510/289‑1 (MQ=255)
                                                                                                                              gtgcacgtaccttgcgggCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCg                                                                                                                                                            <  1:228729/273‑1 (MQ=255)
                                                                                                                                   cgtcacTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAAcc                                                                                                                                                       <  1:456470/289‑1 (MQ=255)
                                                                                                                                      cacTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTCGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCtt                                                                                                                                                     <  2:151614/288‑1 (MQ=255)
                                                                                                                                                  cGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGtt                                                                                                                                                                                                           <  1:465569/222‑1 (MQ=255)
                                                                                                                                                  cGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGtt                                                                                                                                                                                                           >  2:465569/1‑222 (MQ=255)
                                                                                                                                                                   ttACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCg                                                                                                                                                                                                         >  1:135831/1‑207 (MQ=255)
                                                                                                                                                                   ttACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCg                                                                                                                                                                                                         <  2:135831/207‑1 (MQ=255)
                                                                                                                                                                        aTGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         <  2:271792/234‑1 (MQ=255)
                                                                                                                                                                        aTGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCcat                                                                                                                                                                         >  1:271792/1‑234 (MQ=255)
                                                                                                                                                                          gCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCgag                                                                                                                >  1:308291/1‑289 (MQ=255)
                                                                                                                                                                          gCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCgag                                                                                                                >  1:274676/1‑289 (MQ=255)
                                                                                                                                                                                ggTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGAttt                                                                                                            >  1:47669/1‑287 (MQ=255)
                                                                                                                                                                                      aCAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   >  2:13202/1‑226 (MQ=255)
                                                                                                                                                                                      aCAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGAt                                                                                                                                                                   <  1:13202/226‑1 (MQ=255)
                                                                                                                                                                                          ccGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAACCTGCCTTTCCATCATGCTTGGTTCGTACCCTTGTTCATAGAACGCTGCCAGCATGGTTTGAACGGCGAGATTTTTTGCCAACGa                                                                                                 >  1:172311/1‑288 (MQ=255)
                                                                                                                                                                                          ccGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGaa                                                                                                >  2:321075/1‑289 (MQ=255)
                                                                                                                                                                                           cGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGCTTGGGTCGTAACCTTTTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGAtttttcgaacacgca                                                                                                >  1:20859/1‑277 (MQ=255)
                                                                                                                                                                                            gTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAAt                                                                                               >  2:51747/1‑288 (MQ=255)
                                                                                                                                                                                                       aTTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGt                                                                                     <  2:308291/287‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACaa                                                                                                    >  1:293672/1‑263 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACaa                                                                                                    <  2:293672/263‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAg                                                                                           >  2:433058/1‑272 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAg                                                                                           <  1:433058/272‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAatgatg                                                                           >  2:126219/1‑288 (MQ=255)
                                                                                                                                                                                                                 ccccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       <  2:64487/147‑1 (MQ=255)
                                                                                                                                                                                                                 ccccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCa                                                                                                                                                                                                                       >  1:64487/1‑147 (MQ=255)
                                                                                                                                                                                                                  cccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCt                                                                                                                                       <  2:47233/226‑1 (MQ=255)
                                                                                                                                                                                                                  cccTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCt                                                                                                                                       >  1:47233/1‑226 (MQ=255)
                                                                                                                                                                                                                    tttcttcGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGAGGACTGCCAGTGTGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGAc                                                                      <  2:427718/288‑1 (MQ=255)
                                                                                                                                                                                                                                 aCAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCtt                                                         <  1:51747/289‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGccc                                                    >  1:99036/1‑288 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTAt                                                              >  2:199580/1‑256 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTAt                                                              <  1:199580/256‑1 (MQ=255)
                                                                                                                                                                                                                                                                      agAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgggg                       <  2:47669/286‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGc            <  2:172311/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        cAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATtcatca  <  2:274676/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                             
AGGCGTAATCACTTTTGCAACGTTCTCATAGCCTTGCTTCATAGCAAGCCCTGGCAGGGCGGCAACGGGATGCAGTTTATGCGCTTTGCAGAACGCCAGGTAGGCTTCGTGACGTCGCTTACCGGTTGTCACGTCACTGTGCGGCACGCCGAGATAACTGATATTACGATGCCCCTGGTCATACAGCCGTTGCATCAGGATTTTGATTGCCCCTTCGTCGTCATAACAGACCGAAGCAAAGCCTTTTGCGTCACGCGCCAGCAGAACCAGCGATGACTGCCAGTGGGCTAACATTTCTTCTGTTATGCCAGTAAAACCGAACAGCACTACGCCGTCGATATTACGCCGTTTCAGCACTCCCAAATGTTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCA  >  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: