Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,606,034 G→A noncoding (48/87 nt) leuV ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,606,0340GA100.0% 43.7 / NA 16noncoding (48/87 nt)leuVtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/9);  total (7/9)

CTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCT  >  NC_000913_3_pae_tpiA/4605891‑4606091
                                                                                                                                               |                                                         
cTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                      <  2:370833/149‑1 (MQ=255)
     aTCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                 >  1:7529/1‑149 (MQ=255)
          gCCTAGCACTTATTGCTGGCATCATGCTGAGTATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                            <  2:198347/149‑1 (MQ=255)
               gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCACCACGTCCGTAAGAACACTAACGCCTGAAGCTAg                                       >  2:504332/1‑149 (MQ=255)
                cACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                      <  1:459846/149‑1 (MQ=255)
                          tGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGGAGTAAAGATAATTGGTGCGAAGGGGGAGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGATAGCGCGTCTAcga                            >  1:197578/1‑147 (MQ=37)
                                 atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                     <  2:334684/149‑1 (MQ=255)
                                     tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAAAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                 >  2:438504/1‑149 (MQ=255)
                                     tgatgaTAAGGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                 <  1:71844/149‑1 (MQ=255)
                                           tAATGAATTGGAATCAGGGATGAAAAATGAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc           <  2:247543/149‑2 (MQ=255)
                                               gAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt       >  2:89733/1‑149 (MQ=255)
                                                    tGAAGCAGGGATTACAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGCACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGtgcgag  <  1:480722/149‑6 (MQ=17)
                                                                             agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGtgc     >  1:141927/1‑119 (MQ=255)
                                                                             agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGtgc     <  2:141927/121‑3 (MQ=255)
                                                                                      aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                              >  1:26831/1‑87 (MQ=255)
                                                                                      aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                              <  2:26831/87‑1 (MQ=255)
                                                                                                                                               |                                                         
CTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCT  >  NC_000913_3_pae_tpiA/4605891‑4606091

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