Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,571,280 C→T E99K (GAG→AAG)  yjiR ← fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,2800CT94.7% 114.1 / ‑5.0 38E99K (GAG→AAG) yjiRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/27);  minor base G (0/2);  total (9/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTT  >  NC_000913/4571035‑4571525
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ggggTTATGGCAGCCGCTGGTGATGAGGATGTCATCGGGGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCCCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCa                                                                                                                                                                                                                                                   <  1:91763/248‑1 (MQ=255)
 gAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGAACACCGAGCCGCTGTCGAGCATCAAGCGGGGCATCTGCTCTCGCAAAACTCGCTGACCGGCCAACAAGTCATAACCGAGAAACGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCACCCCAGAGGGTTTAGGCTGGGCGCTTCGGCATACGGGGTGCTTTTGCTTAACGGGACAATGGGACTGTCGCTATGCGCCTTCAGCat                                                                                                                                                                                                                                                  >  1:86959/1‑251 (MQ=255)
    ttATGCCAGCCGCTGGTGATGATGATTTCATCGGCGGTGACCACCGAGCCGCTGTCGAGTATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGTTTTGCAGTTTATGCTGCACCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCatat                                                                                                                                                                                                                                                <  1:41433/250‑1 (MQ=255)
                 ggggtagtatgagttcagcggcggTGACCACCGAGCCGCTGTCGAGCTTCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATACCCGAGAACGTTTTGCAGATTATGCGGCACCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGGACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATca                                                                                                                                                                                                                                  <  2:108636/233‑2 (MQ=255)
                  ggggatgatgatGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATc                                                                                                                                                                                                                                   <  1:126979/246‑1 (MQ=255)
                                       ggTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgc                                                                                                                                                                                                            >  2:215238/1‑251 (MQ=255)
                                                       cTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGt                                                                                                                                                                                            >  1:265139/1‑251 (MQ=255)
                                                            gAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcg                                                                                                                                                                                       >  2:114841/1‑251 (MQ=255)
                                                                    ggCGGGCAACTTGCCCCCGCAATACTCGCTGCCCGGCTAAAAAGTCATAACCGAGAACGGTTTGCAGATTATGCGGCACCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACt                                                                                                                                                                               <  2:80102/251‑1 (MQ=255)
                                                                      cGGGCAATCTGCTCGCGCAATACTCGCTGACCGGCTAAAAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGGCGCTATGCGCCTGCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTgg                                                                                                                                                                             <  1:157650/251‑1 (MQ=255)
                                                                         gCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTgg                                                                                                                                                                                                                            >  1:257884/1‑201 (MQ=255)
                                                                              ccgctctccaaatactcgcttcatggCTAACAAGGCAAACCCGAGACCAGTTGGGAGAGTATGCT‑CCCCACCCCG‑GCTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCTGGACTGGCGGCTggg                                                                                                                                                                     <  2:279994/227‑1 (MQ=255)
                                                                                     cGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACAggcgggcgggggtttgg                                                                                                                                                              >  2:314573/1‑239 (MQ=255)
                                                                                      gCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGATGGGGttttagt                                                                                                                                                             >  1:258130/1‑247 (MQ=255)
                                                                                               gCTGACGGGCAAACAAGTGAGACCCCAGAGCGGTTTGCGGATTATGCGGC‑CCCCCC‑GCGCTTACCCCCCGCCGCAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGCTTTACGTTGTGCG‑ACa                                                                                                                                                   <  2:69438/251‑1 (MQ=255)
                                                                                                     cGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACt                                                                                                                                                                                                                              >  1:139598/1‑171 (MQ=255)
                                                                                                     cGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACt                                                                                                                                                                                                                              <  2:139598/171‑1 (MQ=255)
                                                                                                            cAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAAc                                                                                                                                      <  1:215238/251‑1 (MQ=255)
                                                                                                                    aaCCGAGAACGGTTTGCAGATTATGCGGCACCGCCC‑G‑GCTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTGCTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGc                                                                                                                              <  1:378592/251‑1 (MQ=255)
                                                                                                                       cGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGCCTTTCGTTGTGGGGACAAAAAAACCC‑CACGCGGCTGCggg                                                                                                                           >  1:229676/1‑250 (MQ=255)
                                                                                                                                            gctgcACCCCCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTcc                                                                                                      <  1:5997/251‑1 (MQ=255)
                                                                                                                                                               aGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTGCGACATCCGGTGTGCTTTTGCTTAACGGAACAATTGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTTGGTAATTGCCACCGGGCGCTGCACCGGACGCGTCAGCGGCGGTACGGGCGGCTGGGCTTTACGTTGTTCG‑ACAAAATAACCCGAACGCGTCTGCGGAGTGATGAGCTGCATCGTCTCCAGCGTCTGATACGCct                                                                                      <  1:382196/250‑1 (MQ=255)
                                                                                                                                                                 cccccGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCCCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgct                                                                                   <  1:372074/249‑1 (MQ=255)
                                                                                                                                                                     cGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCAc                                                                               <  1:314573/251‑1 (MQ=255)
                                                                                                                                                                          caGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGc                                                                          <  2:206949/251‑1 (MQ=255)
                                                                                                                                                                          caGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGc                                                                          <  1:322383/251‑1 (MQ=255)
                                                                                                                                                                           aGCGGTTTCATGCTGGACGCTTTGACATCCGGTGTGCTTATGCTTAACCGAACAATGGAACTGACACTATGCGCCTTCAGCATATCCAGCTCCTGATCCCACTGGGTAATTTCCACCGTGCGGTGCACCGGATGCGTCATCGGCGGTACTGGCGGCTGGGCTTTGCGTGGGGCG‑ACAAAATAACCCGAACGCGGCTGCTGAGTCATGAGCTGCATCATCTCCAGCGTCTGAGACACCTGCTGCTCGGTGCt                                                                         <  2:48333/251‑1 (MQ=255)
                                                                                                                                                                                         gggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCCCCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGcc                                                             <  2:130842/249‑1 (MQ=255)
                                                                                                                                                                                                tCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCTGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTcc                                                      <  1:190976/249‑1 (MQ=255)
                                                                                                                                                                                                tCGACAGCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGGCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACGGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACGGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCGCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCt                                                     <  2:242883/250‑1 (MQ=255)
                                                                                                                                                                                                                           ggAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCAGCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAttt                         <  1:201671/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAATGGAACTGTCGCTATGCGCCTGCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACTCGGCTGCGGAGTTATGAGCTTCATCGTCTCCa                                                                                                     <  1:96382/171‑1 (MQ=255)
                                                                                                                                                                                                                                aaTGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAAcg                                                                                                                                     <  2:219471/138‑1 (MQ=255)
                                                                                                                                                                                                                                            gCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                                                <  2:204943/195‑1 (MQ=255)
                                                                                                                                                                                                                                            gCTATGCGCCTTCAGCAGATCCAGCACCTGAGCCCACTGGGTAATTTCCTCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACa        <  2:157700/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                               <  2:385089/226‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgcCTTCAGCATATCCAGCCCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTGAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCtt  <  1:54636/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgcCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCtt  <  1:206687/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCC‑G‑GCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCG‑ACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTT  >  NC_000913/4571035‑4571525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑