Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,749,174 C→T P200S (CCG→TCG)  sgbH → 3‑keto‑L‑gulonate‑6‑phosphate decarboxylase SgbH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,1740CT93.9% 103.5 / ‑5.7 33P200S (CCG→TCG) sgbH3‑keto‑L‑gulonate‑6‑phosphate decarboxylase SgbH
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (10/21);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACGCTGGACGACGCCCGCGACTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATCCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCG  >  NC_000913/3748927‑3749414
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
gacgcgggacgCCGCCCGCGACTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCgg                                                                                                                                                                                                                                                <  1:331978/251‑1 (MQ=255)
       gacgacGCCCGCGACTGGCATCGTGTTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTTACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGt                                                                                                                                                                                                                                         <  2:284903/251‑1 (MQ=255)
                 gcgACTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑At                                                                                                                                                                                                                              <  2:326885/251‑1 (MQ=255)
                   gACTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑Attt                                                                                                                                                                                                                            <  1:228922/251‑1 (MQ=255)
                     cTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCAAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGGTTAACCCCG‑CTGAGCCGCCGCTGTTTTAATATAACCCCGTGGAAGCGGTTTTTGTCGGGCGGGGCCCGGCACGGCGGGGCAAATC‑GGCGCAAAtggggggagggttct                                                                                                                                                                                                                           >  2:205738/1‑239 (MQ=255)
                          aTCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATgcg                                                                                                                                                                                                                     <  1:141988/251‑1 (MQ=255)
                              tATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTAACCCTGACT‑ACCTGCCGCCGTTTAAAGATATACGCGTGGAAGCGTTTTTTTCCGGGCGGGAACTGGCA‑GGCGCGGCAAATCTGGGGCAAGGGGGTGGCGGATTTCC‑TGCGCac                                                                                                                                                                                                                  >  1:218923/1‑249 (MQ=255)
                                           cagggcattgagaatcGCGGTCGTGAGGGACAGGCGAGCGGGCAAAAGGGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGGGTTTATTGTCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTg                                                                                                                                                                                                    <  1:44610/238‑1 (MQ=255)
                                                             ggTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCCG‑CTGAACTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGGGGGCACTGGCA‑GGCGGGGGAAAATCGGCGCAAGTGGCGGGCG‑ATTTTCATGCGCAAATCGAGACTATTTGGGGTGGGGcgcgggggg                                                                                                                                                                                   >  1:307420/1‑245 (MQ=255)
                                                                cGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAat                                                                                                                                                                                <  1:150697/250‑1 (MQ=255)
                                                                      gCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGt                                                                                                                                                                         <  2:309404/251‑1 (MQ=255)
                                                                         caGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑AGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAg                                                                                                                                                                      <  1:371335/251‑1 (MQ=255)
                                                                               aGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGAttt                                                                                                                                                                <  2:310695/251‑1 (MQ=255)
                                                                                 cGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAAGCATCAGTTAGGGATTTAt                                                                                                                                                              <  1:108307/251‑1 (MQ=255)
                                                                                     cAACTGTGGGGCGACGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGAGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGaaa                                                                                                                                                          <  2:5174/251‑1 (MQ=255)
                                                                                               gCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGc                                                                                                                                                <  2:335801/251‑1 (MQ=255)
                                                                                                                                                      cATTACTGGCGGGAGTTCCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTGTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAATATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGAt                                                                                         <  1:377758/251‑1 (MQ=255)
                                                                                                                                                           cTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGGTTTGATTTTGt                                                                                    >  1:179443/1‑251 (MQ=255)
                                                                                                                                                             ggCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑GTTTCCATGCGCAAATCGACGCTATCTGGGGAGGGGCGCGTGCGTAATCCTCAGATAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCAGGGGCGGCTGGTGATGGCaaaaaaacgagggtttgattttgggg                                                                                  >  2:109282/1‑230 (MQ=255)
                                                                                                                                                                gggATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGaaa                                                                               <  1:164393/251‑1 (MQ=255)
                                                                                                                                                                               tgACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGGTGa                                                                  >  1:373115/1‑250 (MQ=255)
                                                                                                                                                                                        gCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACcgacga                                                         >  1:110707/1‑250 (MQ=255)
                                                                                                                                                                                           gTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAAc                                                       <  1:351319/249‑1 (MQ=255)
                                                                                                                                                                                                      tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCt                                          >  1:393561/1‑251 (MQ=255)
                                                                                                                                                                                                        cgcgTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTg                                        <  1:378639/251‑1 (MQ=255)
                                                                                                                                                                                                          cgTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGAc                                      >  1:330699/1‑251 (MQ=255)
                                                                                                                                                                                                           gTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACt                                     <  1:205738/251‑1 (MQ=255)
                                                                                                                                                                                                             gAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGTCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGTGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGGGAAAGGTCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTGTATTTTGTGGAAAGGTCAGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTgg                                   <  1:396496/251‑1 (MQ=255)
                                                                                                                                                                                                                     ttATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCa                                                                                                           <  1:380550/171‑1 (MQ=255)
                                                                                                                                                                                                                     ttATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCa                                                                                                           >  2:380550/1‑171 (MQ=255)
                                                                                                                                                                                                                           ccGGGCGGGCACTGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTGGATTTTGTCGAAATGTCGGTTGATGaa                                                                 <  1:303053/207‑1 (MQ=255)
                                                                                                                                                                                                                                       tGGCA‑GGCGCGGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTggg         >  2:269396/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gcgGCAAATTCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCcgcg  >  1:55427/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
GACGCTGGACGACGCCCGCGACTGGCATCGTATTGGCGTGCGGCAGGCCATTTATCATCGCGGTCGTGATGCACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTG‑CTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCA‑GGCGCGGCAAATCCGGCGCAAGTGGCTGGCG‑ATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCG  >  NC_000913/3748927‑3749414

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