Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,719,429 G→A intergenic (‑385/‑49) yiaF ← / → yiaG DUF3053 domain‑containing protein YiaF/putative DNA‑binding transcriptional regulator YiaG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,4290GA100.0% 78.0 / NA 24intergenic (‑385/‑49)yiaF/yiaGDUF3053 domain‑containing protein YiaF/putative DNA‑binding transcriptional regulator YiaG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/12);  total (12/12)

GTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCC  >  NC_000913/3719190‑3719637
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
gttgttTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGAGGTAAGAGATATATCGGCATTTATCAGCTCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTAAATAAGCTTCGATGc                                                                                                                                                                                                       >  1:108272/1‑251 (MQ=255)
       cATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGTCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCGGCTGTGTCTGCTTTGCCCGACTATTCTTAAGAAGCTTCGATGCAATTCAc                                                                                                                                                                                                <  1:124386/251‑1 (MQ=255)
          tattaacaaaagacagtgactttGCATTTTTTGTAGCCAGTGATTTACACCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAATACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACATTCTGAATCGAATAATAGCCTGAGCATGCCATGACTTCACCCCGCTGTGTTTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGTAATTTACGa                                                                                                                                                                                              >  2:98681/22‑250 (MQ=255)
                 aaaTGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAg                                                                                                                                                                                      <  1:12126/251‑1 (MQ=255)
                     gATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCCGtgtg                                                                                                                                                                                  >  2:243581/1‑251 (MQ=255)
                            ccGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGGCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGCTTTg                                                                                                                                                                             >  2:168052/1‑249 (MQ=255)
                                             gCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTGCGAAATCATTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGAGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTGCGATGCAATTCACGATCCCGCAGTGTGATTTGAAGAGTTTTCAATGGAata                                                                                                                                                          <  2:265733/251‑1 (MQ=255)
                                                           ccTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGc                                                                                                                                            >  1:253141/1‑251 (MQ=255)
                                                              tCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGGTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCa                                                                                                                                           <  2:19039/249‑1 (MQ=255)
                                                                                       cTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAAc                                                                                                                    >  1:254589/1‑247 (MQ=255)
                                                                                                    gAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGAt                                                                                                   <  1:58995/251‑1 (MQ=255)
                                                                                                                  atatCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGtt                                                                                                                                                                      >  2:35035/1‑170 (MQ=255)
                                                                                                                  atatCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGtt                                                                                                                                                                      <  1:35035/170‑1 (MQ=255)
                                                                                                                   tatCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGCCAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAACAtt                                                                                    >  2:274806/1‑251 (MQ=255)
                                                                                                                      cGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTAcc                                                                                 >  2:131481/1‑251 (MQ=255)
                                                                                                                                                       cACGCGAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAAt                                                 <  2:77497/250‑1 (MQ=255)
                                                                                                                                                        aCGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTg                                               >  1:112370/1‑251 (MQ=255)
                                                                                                                                                             aaaTGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGc                                            <  1:131481/249‑1 (MQ=255)
                                                                                                                                                                  aCAGGCTGAATCGAATAATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTCATACGCTTCGATGCAATTCACGAACCCGCAGTGAGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGCTTGTTTTTAAAGATGAAACGCAGAAAATTACCCAGCCGCACAGAAACACGTGCTGTACCGAAATTGAGCAGCTAcc                                     >  2:150854/1‑250 (MQ=255)
                                                                                                                                                                              gaatcATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAAATGAGCAGTTACGAAAAGGGACAgg                         <  2:227640/251‑1 (MQ=255)
                                                                                                                                                                                  cATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATGCACGATCCCGCAGTGTGATTGGAGGAGTGTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGTAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTa                     <  1:47490/251‑1 (MQ=255)
                                                                                                                                                                                   aTAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCGTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTaa                    <  1:326055/251‑1 (MQ=255)
                                                                                                                                                                                                  ccTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGCAATTGAGCAGTTACGAAAAGGGACAGGCTTAAACATCGATGATTTCg     >  2:221751/1‑251 (MQ=255)
                                                                                                                                                                                                     gACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGccc  <  2:1300/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
GTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCC  >  NC_000913/3719190‑3719637

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