Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,998,712 (C)5→4 coding (117/795 nt) ygeQ ← protein YgeQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,998,7080C.91.3% 79.3 / 2.1 23coding (121/795 nt)ygeQprotein YgeQ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (10/11);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AGCAATCCTTGAATTTCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGCCCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCTTTCATAAAATATT  >  NC_000913/2998578‑2998836
                                                                                                                                  |                                                                                                                                
aGCAATCCTTGAATTTCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGcct                                                                                                                                >  1:442310/1‑132 (MQ=255)
 gCAATCCTTGAATTTCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGcccc                                                                                                                               <  2:330392/133‑1 (MQ=255)
  cAATCCTTGAATTTCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGccccg                                                                                                                              <  1:32180/133‑2 (MQ=255)
             tttCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGa                                                                                                                  <  1:27465/133‑1 (MQ=255)
                  tCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGtt                                                                                                             <  2:61766/133‑1 (MQ=255)
                        gTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGCCCCCGTGGAATTAAGAATGTTAAATa                                                                                                        <  1:106235/133‑1 (MQ=255)
                                   ctTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCg                                                                                            >  1:71172/1‑133 (MQ=255)
                                     ttATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTc                                                                                          <  2:71172/133‑1 (MQ=255)
                                            aTTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTg                                                                                   <  1:432799/133‑1 (MQ=255)
                                                      ttGGCTGAGTCAATCGGAAGACGATAATAAGGTTACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGc                                                                         <  2:180207/133‑1 (MQ=255)
                                                          cTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTgg                                                                     >  1:60068/1‑133 (MQ=255)
                                                                   tCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCg                                                            >  2:455831/1‑133 (MQ=255)
                                                                   tCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCg                                                            <  2:19508/133‑1 (MQ=255)
                                                                                  tAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCa                                             <  1:9116/133‑1 (MQ=255)
                                                                                    ggATACTTCTCTGGAATTCACGCTAAATAAGTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGc                                           <  1:135476/133‑1 (MQ=255)
                                                                                           tATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCCCTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCaga                                    <  2:442310/133‑1 (MQ=255)
                                                                                                     tAACGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATg                          >  2:51496/1‑133 (MQ=255)
                                                                                                      aaCGATAAATAACTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGTAAATATAACAGCGACCAGAAGAGTTAATGc                         >  1:236798/1‑133 (MQ=255)
                                                                                                       aCGATAAATTCCTTTTATCTGTTG‑TTCCCCCGGGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCa                        <  2:236798/133‑1 (MQ=255)
                                                                                                            aaataaCTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCACACGAAAAGATAGCAGGGAGCAGAAGAGTTAATGCAGAttt                   >  2:152360/1‑133 (MQ=255)
                                                                                                              ataaCTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTcc                 >  1:437769/1‑133 (MQ=255)
                                                                                                                 aCTTTTATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCtt              <  2:60068/133‑1 (MQ=255)
                                                                                                                 aCTTTTATCTATTGCTG‑CCCCGTGTAATTAAGTATGTTAAAGACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCtt              <  2:210881/133‑1 (MQ=255)
                                                                                                                     ttATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTATTGAAGATTTACCTTTCAt          >  1:210881/1‑133 (MQ=255)
                                                                                                                     ttATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAACATACCAGCGAGCTGAAGAGTTAATGGAGATATCCCTTTCAt          >  2:140332/1‑133 (MQ=255)
                                                                                                                      tATCTGTTGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCTTTCATa         <  2:508943/133‑1 (MQ=255)
                                                                                                                            ttGCTG‑CCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCACCGAGCAGAAGAGGTAATGCAGATTTACCTTCCATAAAatat   >  2:217306/1‑133 (MQ=255)
                                                                                                                             tgctg‑cCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCTTTCATAAAATAtt  >  2:80859/1‑133 (MQ=255)
                                                                                                                                  |                                                                                                                                
AGCAATCCTTGAATTTCGTCTTTGGTATAGTGGCTCTTTATTTCATTTTGTAATTTGGCTGAGTGAATCAGAAGACGATAATTAGGATACTTATCTGGATTTAACGATAAATAACTTTTATCTGTTGCTGCCCCCGTGGAATTAAGAATGTTAAATACAGATTCACTCGTCACTTCTGCCCCTGTCGCCTGGCATGTGCCGCACGAAAATATACCAGCGAGCAGAAGAGTTAATGCAGATTTCCCTTTCATAAAATATT  >  NC_000913/2998578‑2998836

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