Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,325,054 G→A W103* (TGG→TGA)  yciQ → DUF2207 domain‑containing protein YciQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,325,0540GA88.0% 51.9 / 1.2 25W103* (TGG→TGA) yciQDUF2207 domain‑containing protein YciQ
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (10/12);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

ATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTT  >  NC_000913/1324827‑1325282
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGtt                                                                                                                                                                                                               <  1:358143/251‑1 (MQ=255)
             aagaagGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACTAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTAACGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATTGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATCGTc                                                                                                                                                                                                  >  1:88815/1‑251 (MQ=255)
              agaagGGAGCTCGCTCTTTGAGGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGATATCAAAGAAAATAGTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTGTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGGTGGCGTCACGCGTGATGTTATTCCTGAATCAAGGCATCTTGACGGGCAAGTCGTGTTAATTAGTATTGTCg                                                                                                                                                                                                 <  2:245640/251‑1 (MQ=255)
                                    gtggAAGCCACTCTGCAAACAGATGGGGTGCGATACAGTAAAGAAAATATTCATTTTCAGGCGCGAAAGCAGTAGATTAAGCTCTGCTTTTATCGTGATTTACCACGACTATGGATGCATGCTGATGGTTACGCTGCACTGCTGAACTATCATATTGGTTGCGTCACCCGTGATGGTATTCCTGATCCCTGACAGCTTGTCTTGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCt                                                                                                                                                                            <  2:221860/247‑1 (MQ=255)
                                    gTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCtt                                                                                                                                                                           <  1:93683/251‑1 (MQ=255)
                                     ttGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTg                                                                                                                                                                                                    <  2:97205/225‑1 (MQ=255)
                                     ttGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTg                                                                                                                                                                                                    >  1:97205/1‑225 (MQ=255)
                                                       cAGATGGGGTGCTTGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATTGGGACGCTGCTCGGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAATGCGACTATCAt                                                                                                                                                        <  2:56325/251‑1 (MQ=255)
                                                            ggggTGCTTGACATCAAAGCAAATATTCATTTTCAGGCGCGAAATCAGCAGAATAAGCACGGCTTTTAACGTGATTGACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTAGCATATTGTTGGCGTCACCCGTGATTGTAGTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATCGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCa                                                                                                                                                   <  2:88815/251‑1 (MQ=255)
                                                            ggggTGCTCGACAACAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATc                                                                                                                                                    <  1:288263/250‑1 (MQ=255)
                                                                             aGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGTGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGtt                                                                                                                                   >  1:262666/1‑250 (MQ=255)
                                                                                    atTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAACATGc                                                                                                                           >  1:358125/1‑251 (MQ=255)
                                                                                              aGGCGCGAAAACAGCAGATTAAGCACGGCTTTTATCGTGATTTACCTCGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGGCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGt                                                                                                                 <  1:239177/251‑1 (MQ=255)
                                                                                                               aTTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCt                                                                                                >  1:85822/1‑251 (MQ=255)
                                                                                                                  aaGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCt                                                                                             >  1:202853/1‑251 (MQ=255)
                                                                                                                       cGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGAt                                                                                                                                                                                         <  1:114290/154‑1 (MQ=255)
                                                                                                                       cGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGAt                                                                                                                                                                                         >  2:114290/1‑154 (MQ=255)
                                                                                                                       cGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAAtt                                                                                                                                               <  1:313680/196‑1 (MQ=255)
                                                                                                                             ttATCGTGAGTTACCACGACTATGGTCGCAGCCTGATGGGGACGCTGCTCTGCTGAACTACCATATTGTTGGCGTCGCCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGTCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTTCTAATCTGGAAc                                                                                   <  2:260258/250‑1 (MQ=255)
                                                                                                                                             cGACTATGGATGCAGCCTGTTGGGGACGCTGAACTGCTGAAATATCATATTGTTGGCGTCCCCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTAAATGCGTCTTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGAATATCATTATCACATTCATACTCAGGTTAATAATGCTTTCCTGAGTGAGGGGGATTCGGATCTGCTAATCTGGACCGTGACGGGTACCCACTg                                                                  >  2:229461/1‑251 (MQ=255)
                                                                                                                                                              tGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGCGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTtat                                                  >  1:66322/1‑250 (MQ=255)
                                                                                                                                                                                                        cACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAAt        <  1:236117/250‑1 (MQ=255)
                                                                                                                                                                                                          cccGTGATGGTATTCCTGAACCCTGACATCTTGCCTGGCATATCGGGTTAATGCGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGt                                                             >  2:364649/1‑195 (MQ=255)
                                                                                                                                                                                                          cccGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGt                                                             <  1:364649/195‑1 (MQ=255)
                                                                                                                                                                                                             gTGATGGAATTCCTGAACACTGACATATTGACTGGCATATAGGGTTAATGAGTATTGTCGTGGGCGATCAACAACGTTTATTGCCTCCAGGCGACAATCATTATACAATTAATTATCAGGGTAAACATGCTTTCATACGTGAGAGGGATTCTTATCTGCTAATCTGGAACGTGACCGGTAACCACTTGCCGTTTGAAATTTATTAGACCCGTATTTATCTCCAGTTCTCTAATATTGCGGGTATTGCAttt  >  2:320045/1‑251 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
ATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTT  >  NC_000913/1324827‑1325282

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