Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,257,400 A→G intergenic (‑65/+29) psuT ← / ← psuG putative pseudouridine transporter/pseudouridine‑5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,257,4000AG90.9% 27.8 / ‑3.9 11intergenic (‑65/+29)psuT/psuGputative pseudouridine transporter/pseudouridine‑5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (7/3);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAACAGAAATGCTATTGCCAGTAACACCACCATCCCCACAACACTTCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAG  >  NC_000913/2257274‑2257535
                                                                                                                              |                                                                                                                                       
gACAACAGAAATGCTATTGCCAGTAACACCACCATCCCCACAACACTTCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCtt                                                                                                                            >  2:104865/1‑140 (MQ=255)
                            gaactgtCCCCACAACACTTCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCgaga                                                                                                <  1:227527/134‑1 (MQ=255)
                                  tCCCCACAACACTTCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGa                                                                                          >  1:318726/1‑140 (MQ=255)
                                              ttCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCa                                                                              >  2:62397/1‑140 (MQ=255)
                                               tctcATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCaa                                                                             <  1:62397/140‑1 (MQ=255)
                                                        atCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTtcgctcg                                                                    >  2:266398/1‑140 (MQ=255)
                                                            aTAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAg                                                                <  2:318726/140‑1 (MQ=255)
                                                             tAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGa                                                               >  2:134204/1‑140 (MQ=255)
                                                                                                cAGCGCCAGGTTTTGGTATGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGAttt                            <  1:266398/140‑1 (MQ=255)
                                                                                                                  tGCTTGATGAGTGCGGGCGACGGCTTTCTGCCCGTAAAAAATTACCCCGCGAGACGCTGATATTCTTTGGCATTTTCGCTCGCAAGAATGGCGTTGTTGAACACCAGCGGGATGTTGGATTTCAGGCTGTCACCGCCGGt          <  2:101906/140‑1 (MQ=255)
                                                                                                                          gaGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAg  >  1:101906/1‑140 (MQ=255)
                                                                                                                          gaGTGCGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGGTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAg  >  1:327707/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                                       
GACAACAGAAATGCTATTGCCAGTAACACCACCATCCCCACAACACTTCTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAG  >  NC_000913/2257274‑2257535

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