Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 798,688 C→T P35S (CCG→TCG)  pgl → 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913798,6880CT97.9% 168.1 / ‑3.9 47P35S (CCG→TCG) pgl6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/20);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TTTTAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGCCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCCGGGTAG  >  NC_000913/798544‑798833
                                                                                                                                                |                                                                                                                                                 
ttttAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGc                                                                                                                                             >  1:565654/1‑151 (MQ=255)
    aGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAgg                                                                                                                                          >  2:551841/1‑150 (MQ=255)
     gCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTg                                                                                                                                        <  2:51238/151‑1 (MQ=255)
      cTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGc                                                                                                                                       >  2:407548/1‑151 (MQ=255)
          tACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGcc                                                                                                                                   >  2:331216/1‑151 (MQ=255)
             agtcagtTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGAt                                                                                                                                <  1:499937/151‑1 (MQ=255)
             agtcagtTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGAt                                                                                                                                >  2:573321/1‑151 (MQ=255)
                   ttGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGCGCAGCCGATggtggt                                                                                                                          >  1:393249/1‑151 (MQ=255)
                         aaTGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGccc                                                                                                                    <  2:81304/151‑1 (MQ=255)
                              aaaGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGGTGGTGGTCAGCCCGGAc                                                                                                                >  2:779785/1‑150 (MQ=255)
                                      aTTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTAt                                                                                                       <  2:725546/151‑1 (MQ=255)
                                       ttCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTAt                                                                                                       >  1:758606/1‑150 (MQ=255)
                                       ttCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGGGCAGCCGATGGTGGTCAGCCCGGACAAAAGTTAtc                                                                                                      >  1:156931/1‑151 (MQ=255)
                                        tCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGCCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGGCGCCGATGGTGGTCAGCCCGGACAAACGTTAtct                                                                                                     >  1:51692/1‑151 (MQ=255)
                                                cAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGt                                                                                               <  1:452224/149‑1 (MQ=255)
                                                    cAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCCGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGATgggg                                                                                          >  1:590194/1‑148 (MQ=255)
                                                    cAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGgtgt                                                                                         <  1:19415/151‑1 (MQ=255)
                                                      gTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTc                                                                                       <  1:114721/151‑1 (MQ=255)
                                                       tttATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGGAGCCGATGGGGGTCAGCCCGGACAAACGt                                                                                                          >  1:544091/1‑131 (MQ=255)
                                                       tttATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGt                                                                                                          <  2:544128/131‑1 (MQ=255)
                                                             tCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTg                                                                                 >  2:471306/1‑150 (MQ=255)
                                                             tCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACACTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTg                                                                                 >  2:616008/1‑150 (MQ=255)
                                                              cGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTg                                                                                 >  2:124069/1‑149 (MQ=255)
                                                                      cTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGt                                                                       >  2:280279/1‑151 (MQ=255)
                                                                         agagCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTcc                                                                     >  2:848269/1‑150 (MQ=255)
                                                                          gagCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTcc                                                                     >  1:749350/1‑149 (MQ=255)
                                                                          gagCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTg                                                                   <  1:551804/151‑1 (MQ=255)
                                                                           agCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGGGTTCGCCCTGAGTTTCGCGTCCTgg                                                                  >  1:741503/1‑151 (MQ=255)
                                                                            gCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGc                                                                 <  1:554621/151‑1 (MQ=255)
                                                                             ccAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGcg                                                                <  2:653023/151‑1 (MQ=255)
                                                                             ccAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGcg                                                                <  1:115049/151‑1 (MQ=255)
                                                                                    aTTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGcgtatcgta                                                         <  2:409181/151‑1 (MQ=255)
                                                                                              ggAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCgg                                                <  2:212586/150‑1 (MQ=255)
                                                                                                       aTCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCCGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGCGTTTCGCGTCCTGGCGGATCGTATCGCCCCGGACGATGgcgc                                      >  1:501641/1‑151 (MQ=255)
                                                                                                       aTCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGgcg                                       >  2:312426/1‑150 (MQ=255)
                                                                                                         cATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCAc                                    <  2:434226/151‑1 (MQ=255)
                                                                                                          aggAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGGGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTCTGTTGGGGTTCTCCCTGAGGTTCGCCCCCTGGCGTATCGTATCGCCCCGGACGATgg                                         >  1:757373/3‑145 (MQ=255)
                                                                                                          aTGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACt                                   <  2:395998/151‑1 (MQ=255)
                                                                                                             aaGGCGCACTGACGCTGACACAGGTTGTCGATGTGCCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGAc                                >  2:824755/1‑151 (MQ=255)
                                                                                                              aGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGAcc                               >  1:853827/1‑151 (MQ=255)
                                                                                                                 cgcACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCttt                            >  1:270391/1‑151 (MQ=255)
                                                                                                                        aCGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTTAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGc                       >  2:264855/1‑149 (MQ=255)
                                                                                                                             gACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTc                 >  1:337678/1‑150 (MQ=255)
                                                                                                                                  aGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCt           <  2:445100/151‑1 (MQ=255)
                                                                                                                                  aGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCt           <  1:739532/151‑1 (MQ=255)
                                                                                                                                         cGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCCGGGt    <  1:147788/151‑1 (MQ=255)
                                                                                                                                           aTGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCCGGGTAg  <  1:529116/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTTTAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGCCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCCGGGTAG  >  NC_000913/798544‑798833

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