Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,119,307 C→T *47* (TAG→TAA yceO ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,119,3070CT95.0% 112.9 / ‑1.5 40*47* (TAG→TAAyceOuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (21/17);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

TCTTCATAAGTTTAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAACTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTAAAAACGGCC  >  NC_000913/1119174‑1119443
                                                                                                                                     |                                                                                                                                        
tCTTCATAAGTTTAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTaaaa                                                                                                                                     >  2:825196/1‑139 (MQ=255)
     aTAAGTTTAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATaaa                                                                                                                                <  2:703664/139‑1 (MQ=255)
          tttAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAAcaca                                                                                                                           <  1:825196/139‑1 (MQ=255)
            tAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATa                                                                                                                         >  1:282627/1‑139 (MQ=255)
                    tCAGCATCAGAGATGTCTGGTGAGTGCTGTAGTAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagat                                                                                                                      <  2:859358/134‑1 (MQ=255)
                    tCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGTGAGTAAATCTTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagatcagat                                                                                                                 >  2:964118/1‑139 (MQ=255)
                    tCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagat                                                                                                                      >  1:859358/1‑134 (MQ=255)
                             gagaTGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTg                                                                                                        >  2:470210/1‑139 (MQ=255)
                              agaTGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGc                                                                                                       <  1:470210/139‑1 (MQ=255)
                                   tCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCcaaca                                                                                                  >  1:210566/1‑139 (MQ=255)
                                          aGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCa                                                                                           >  2:868345/1‑139 (MQ=255)
                                            tGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATaaa                                                                                                                                <  2:27636/100‑1 (MQ=255)
                                            tGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATaaa                                                                                                                                >  1:27636/1‑100 (MQ=255)
                                            tGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagatcag                                                                                                                   >  1:660129/1‑113 (MQ=255)
                                            tGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagatcag                                                                                                                   <  2:660129/113‑1 (MQ=255)
                                                 tAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGaaaa                                                                                    >  2:1232023/1‑139 (MQ=255)
                                                           aataatCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCa                                                                          >  1:25291/1‑139 (MQ=255)
                                                                    aTTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagat                                                                                                                      >  1:71343/1‑86 (MQ=255)
                                                                    aTTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATagat                                                                                                                      <  2:71343/86‑1 (MQ=255)
                                                                       cAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTcgc                                                              <  1:1232023/139‑1 (MQ=255)
                                                                            tttCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAata                                                         <  1:868345/139‑1 (MQ=255)
                                                                                aCGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCACGAAtattat                                                     <  2:282627/139‑1 (MQ=255)
                                                                                                agaTTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCaa                                     >  1:187987/1‑139 (MQ=255)
                                                                                                  aTTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATa                                                                               >  1:345501/1‑95 (MQ=255)
                                                                                                  aTTAAATATTTGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATa                                                                               <  2:345501/95‑1 (MQ=255)
                                                                                                           ttGTCAGTTTATGTCAGTGCGGAAAATAAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCa                          >  1:869799/1‑139 (MQ=255)
                                                                                                            tGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  >  2:1145784/1‑114 (MQ=255)
                                                                                                            tGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  >  2:371367/1‑114 (MQ=255)
                                                                                                            tGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  <  1:1145784/114‑1 (MQ=255)
                                                                                                            tGTCAGTTTATGTCAGTGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  <  1:371367/114‑1 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  >  2:1244889/1‑98 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATtga                                                  <  1:1244889/98‑1 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGa                              <  1:792213/118‑1 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGa                              >  2:792213/1‑118 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGAc                             <  2:1205781/119‑1 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGAc                             >  1:1205781/1‑119 (MQ=255)
                                                                                                                            tGCGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTaaa         >  1:823727/1‑139 (MQ=255)
                                                                                                                            tGCGGAAAACTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTaaa         >  1:640161/1‑139 (MQ=255)
                                                                                                                              cGGAAAATTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTaaaaa       <  1:964118/139‑1 (MQ=255)
                                                                                                                                   aaCTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTAAAAACGGcc  <  2:640161/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
TCTTCATAAGTTTAGTGGGTTCAGCATCAGAGATGTCTGGTGAGTGCTGTAGAAGCAGGAATAATCTCATTCAATTTTTCACGCTAAATGAAAGGGAGATTAAATATTTGTCAGTTTATGTCAGTGCGGAAAACTAAAAATAAAACACATAGATCAGATCCATAATTGCCAACAATCCCCAGAGAAAAAGATAAAGCATCAGATGCTCGCGAATATTATTGATGAGATAGTGCAACAGCCGACGCATTAAATACTCCTGTAAAAACGGCC  >  NC_000913/1119174‑1119443

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