Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,841,612 G→T P41P (CCG→CCT nudG → 5‑hydroxy‑CTP diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,841,6120GT100.0% 115.7 / NA 34P41P (CCG→CCTnudG5‑hydroxy‑CTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/18);  total (16/18)

TGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCGGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTC  >  NC_000913/1841480‑1841733
                                                                                                                                    |                                                                                                                         
tgtgtAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGaaa                                                                                                                    <  1:448643/140‑1 (MQ=255)
 gtgtAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAg                                                                                                                   >  2:231872/1‑140 (MQ=255)
 gtgtAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAg                                                                                                                   >  1:506052/1‑140 (MQ=255)
  tgtAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGt                                                                                                                  >  2:284292/1‑140 (MQ=255)
  tgtAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTc                                                                                                                 >  2:35411/1‑140 (MQ=255)
          aTGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCa                                                                                                          >  1:109260/1‑140 (MQ=255)
          aTGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCa                                                                                                          >  2:556450/1‑140 (MQ=255)
           tGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAg                                                                                                         >  2:44050/1‑140 (MQ=255)
               aaTGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGgcgc                                                                                                     >  1:116448/1‑140 (MQ=255)
                     tGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGc                                                                                               <  1:556450/140‑1 (MQ=255)
                          ttgttgCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAg                                                                                          <  2:282147/140‑1 (MQ=255)
                            gttgCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGtt                                                                                        >  1:244017/1‑140 (MQ=255)
                                 cgccATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTAcgcg                                                                                   >  1:135369/1‑140 (MQ=255)
                                 cgccATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTAcgcg                                                                                   <  2:343704/140‑1 (MQ=255)
                                      tcatTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCgaagaa                                                                              <  2:406091/140‑1 (MQ=255)
                                                tGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCg                                                                    <  1:573361/140‑1 (MQ=255)
                                                 gATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGa                                                                   >  2:4514/1‑140 (MQ=255)
                                                  aTGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGaa                                                                  <  2:196424/140‑1 (MQ=255)
                                                      cAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCaa                                                              <  1:44050/140‑1 (MQ=255)
                                                         aaTTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACtg                                                           <  1:456745/140‑1 (MQ=255)
                                                          aTTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACtgt                                                          <  2:116448/140‑1 (MQ=255)
                                                            tttACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTgg                                                        <  2:276375/140‑1 (MQ=255)
                                                              tACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGt                                                      <  1:413223/140‑1 (MQ=255)
                                                              tACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGt                                                      <  1:4514/140‑1 (MQ=255)
                                                                     gcAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATg                                               >  1:503564/1‑140 (MQ=255)
                                                                                      gCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCga                              >  1:427571/1‑140 (MQ=255)
                                                                                          tCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGGGAATATGTTGCCAGCCATCAGCGAGAAg                          >  2:215891/1‑140 (MQ=255)
                                                                                                  gATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCgggg                  <  1:284292/140‑1 (MQ=255)
                                                                                                   aTTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGc                 >  1:496257/1‑140 (MQ=255)
                                                                                                    ttATGGGAGTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCg                <  1:271174/140‑1 (MQ=255)
                                                                                                            gTTTGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATcc        <  2:390119/140‑1 (MQ=255)
                                                                                                               tGCCGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACCGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATc     >  1:527464/1‑140 (MQ=255)
                                                                                                                 ccGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCtt   <  1:462210/140‑1 (MQ=255)
                                                                                                                  cGGTGGTAAAGTCGAGCCTGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTc  <  1:446914/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                         
TGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCGCAACGCCCCGCCCAGAGCGATCAGGCGGGATTATGGGAGTTTGCCGGTGGTAAAGTCGAGCCGGATGAAAGTCAGCGGCAGGCGCTGGTGCGTGAGTTACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTC  >  NC_000913/1841480‑1841733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: