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

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