Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 62,467 T→G N41T (AAT→ACT)  SAUSA300_RS00270 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779362,4670TG93.8% 33.0 / ‑4.2 16N41T (AAT→ACT) SAUSA300_RS00270transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (5/10);  minor base A (0/1);  total (5/11)
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.92e-01

TCAATGCTGTAATCTTACGATAACCGTAGGTATATTGTTTACCAACATGTTGTGAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAATTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTGTTGGAAATCCTGCTGCTTTCATTTCAACTGCTTTATA  >  NC_007793/62372‑62555
                                                                                               |                                                                                        
tCAATGCTGTAATCTTACGATAACCGTAGGTATATTGTTTACCAACATGTTGTAAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAAgtttt                                                                                      <  1:976821/100‑1 (MQ=255)
 cAATGCTGTAATCTTACGATAACCGTAGGTATATTGTTTACCAACATGTTGTAAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAAATTTTc                                                                                     <  1:3646151/100‑1 (MQ=255)
                   tctgccGTTGGTATATTGTTTACCAACATGTTGTAAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTa                                                                   <  1:2797804/96‑1 (MQ=38)
                                          cAACATGTTGTAAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTc                                            >  1:5192952/1‑100 (MQ=255)
                                                              ccgatctCCCCACTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCctgctg                        <  1:5706808/94‑1 (MQ=255)
                                                                tcTTTCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTgctgct                      <  1:5686312/99‑1 (MQ=255)
                                                                tcTTTCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTgctgct                      <  1:3382139/99‑1 (MQ=255)
                                                                gatctCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTgctgct                      <  1:4383812/96‑1 (MQ=255)
                                                                 cTTTCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTgctgcta                     <  1:745130/100‑2 (MQ=255)
                                                                 cTTTCCCCATTTCGATAACATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTgctgctg                     <  1:4560042/100‑2 (MQ=255)
                                                                  tttCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCttt                    <  1:4151114/100‑1 (MQ=255)
                                                                   ttCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCtttc                   >  1:5443018/1‑100 (MQ=255)
                                                                   ttCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCCTTTTTTGAAAATCCTGCTGCtttc                   >  1:5527977/1‑100 (MQ=255)
                                                                    tCCCCATTTCGATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCtttca                  >  1:3053264/1‑100 (MQ=255)
                                                                              gATACCATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCTTTCATTTCAACTgg        <  1:1366119/100‑2 (MQ=255)
                                                                                   cATCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCTTTCATTTCAACTGCTTtag   <  1:1819990/100‑2 (MQ=255)
                                                                                    aTCGCCACCAAGTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTTTTGAAAATCCTGCTGCTTTCATTTCAACTGCTTtata  >  1:2678912/1‑100 (MQ=255)
                                                                                               |                                                                                        
TCAATGCTGTAATCTTACGATAACCGTAGGTATATTGTTTACCAACATGTTGTGAAAATCTATAACTTTCCCCATTTCGATACCATCGCCACCAATTTTTCACTTGTGTTCTATTTCTAATATTTAATTCTTTCATAATTTCTTTTGTTGGAAATCCTGCTGCTTTCATTTCAACTGCTTTATA  >  NC_007793/62372‑62555

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