Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 907,626 C→A L44I (CTA→ATA)  SAUSA300_RS04465 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793907,6260CA100.0% 40.9 / NA 17L44I (CTA→ATA) SAUSA300_RS04465membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/1);  total (16/1)

TTGTCAGCGATTATAGGGTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTACTACATAATATTACAACACAAGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACTTTCTTCATCACTTGGATATTTA  >  NC_007793/907536‑907712
                                                                                          |                                                                                      
ttGTCAGCGATTATAGGGTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACataata                                                                               >  1:1715570/1‑100 (MQ=255)
    cAGCGATTATAGGGTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTac                                                                           >  1:2269034/1‑100 (MQ=255)
    cAGCGATTATAGGGTCTATCGTAGGCATTGGAGGCGGTATTATTAACGTTCCAACAATGGTTTACCTCGGAGTTGAACATGGATTAATACATAATATTac                                                                           >  1:2471633/1‑100 (MQ=255)
               gggTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAg                                                                >  1:753192/1‑100 (MQ=255)
               gggTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAg                                                                <  1:5508240/100‑1 (MQ=255)
                     aTCGTAGGCATTGGAGGCGGTATTAATATCGTTCCAACAATGGTTTACCTTGGTGTTGAACATGGAATAATTAATAAAATTAAAAAACAAGTAGcagatc                                                          >  1:2224098/1‑95 (MQ=255)
                        ctAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGa                                                       >  1:532054/2‑100 (MQ=255)
                         tAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGCTAGGGAc                                                      >  1:5272644/1‑100 (MQ=255)
                         tAGGCATTGGAGGCGGTATTATTATCGGTCCAACAATGGTTTACCTCGGGGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGAc                                                      >  1:3937050/1‑100 (MQ=255)
                          aGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACg                                                     >  1:3004884/1‑100 (MQ=255)
                          aGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACg                                                     >  1:2907709/1‑100 (MQ=255)
                          aGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACg                                                     >  1:66227/1‑100 (MQ=255)
                          aGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGGTTATCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACg                                                     >  1:3640655/1‑100 (MQ=255)
                                                        aacaaTGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACtt                       >  1:5738834/1‑100 (MQ=255)
                                                          caaTGGTTTACCTCGGTGTTGAACATGGATTAATACATAATATTACAACACAAGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACTttc                     >  1:4728539/1‑100 (MQ=255)
                                                                gtgcccTCGGTGTTGAACATGGATTAATACATAATATTACAACACATGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACTTTCTtcatc               >  1:1903324/5‑100 (MQ=255)
                                                                             tGAACATGGATTAATACATAATATGACAACACAAGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACTTTCTTCATCACTTGGATATTTa  >  1:5400722/1‑100 (MQ=255)
                                                                                          |                                                                                      
TTGTCAGCGATTATAGGGTCTATCGTAGGCATTGGAGGCGGTATTATTATCGTTCCAACAATGGTTTACCTCGGTGTTGAACATGGATTACTACATAATATTACAACACAAGTAGCGATAGGGACGTCTTCAGTCATTCTAATTGTGACAGGACTTTCTTCATCACTTGGATATTTA  >  NC_007793/907536‑907712

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