Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_pgi 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_pgi3,560,4551.G97.0% 97.4 / ‑3.6 33intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (10/22);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGT  >  NC_000913_3_vch_pgi/3560314‑3560591
                                                                                                                                               |                                                                                                                                        
tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGAt                                                                                                                                    >  1:157208/1‑149 (MQ=255)
 cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATc                                                                                                                                   <  2:85030/149‑1 (MQ=255)
 cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATc                                                                                                                                   >  2:35048/1‑149 (MQ=255)
  gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCAATGACAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCa                                                                                                                                  >  1:291331/1‑149 (MQ=255)
   atcaaCAGCGTTGAGCCATGGGGAATTTGCTCCGCCACTTCGCGGCCGATGCGCTCTCCGTCTTCGGTCTGGGTGGCCTTGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAg                                                                                                                                 <  2:19115/146‑1 (MQ=255)
     aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGt                                                                                                                               <  2:157208/149‑1 (MQ=255)
      aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGtt                                                                                                                              <  1:379059/149‑1 (MQ=255)
       tcAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGGAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGttt                                                                                                                             <  1:219456/148‑1 (MQ=255)
          gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTGTACCGTTCGCACTGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGc                                                                                                                          <  1:544825/149‑1 (MQ=255)
                        gggATTTGCTGCGCCACTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGGTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                            <  2:149056/149‑1 (MQ=255)
                        gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                            >  1:32747/1‑149 (MQ=255)
                        gggATTTGCACCGCCACTTGGCGGGCGATGCGCTCTTGTTCTTCGGGCTGGGTGGGCTTGCGATCGTGCCACGGCGTGTGAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGTATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                            <  2:22796/149‑1 (MQ=255)
                           aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                         <  1:31749/149‑1 (MQ=255)
                             ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                       <  1:434357/149‑1 (MQ=255)
                                   cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGCGG‑CAGGAACAGGTTTTGCCACGCCAGCTCTTTGACGTcgcgtc                                                                                                 >  2:70469/1‑147 (MQ=255)
                                        cTTTGCGGGCGATTCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                            <  1:377806/149‑1 (MQ=255)
                                            tggggCAATGCGCTCTTTTTTGTCGGTCTGCGTGCGCTTACTGTCGTGCCACGGCCTGTTAAGCGAACTGGAAGGCCGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                        <  2:12938/147‑1 (MQ=255)
                                                                    ggTCTGGGCGGCCTTGCGATCGTGCAACGGCGTCTCAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGCCTGCGGGCTGACGGAGAAATGCTCTa                                                                <  2:381895/149‑1 (MQ=255)
                                                                     gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                               >  1:337764/1‑149 (MQ=255)
                                                                     gTCGGGGTGGCCTTGGGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                               <  2:424216/149‑1 (MQ=255)
                                                                     ctttGGGTGGCCTGGCGATCGTGCGACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                               <  2:553884/146‑1 (MQ=255)
                                                                             ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                       <  2:32747/149‑1 (MQ=255)
                                                                               cctgtcGAACGGGCCATTGCGTTTTAACCGAACTGGAAGCCAGCGCTGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGCCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                     <  2:82806/144‑1 (MQ=255)
                                                                                      aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                              >  2:429244/1‑149 (MQ=255)
                                                                                            ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                    <  1:251573/89‑1 (MQ=255)
                                                                                            ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                    >  2:251573/1‑89 (MQ=255)
                                                                                                   gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGAACCGGTGTGACGTCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                 <  1:406598/149‑1 (MQ=255)
                                                                                                        aacgcaaCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGCTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                            <  2:55869/144‑1 (MQ=255)
                                                                                                              aCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                      <  2:332447/149‑1 (MQ=255)
                                                                                                                tGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGGTTTGCTCCGCCAGCTAATTCATGTGGCGGCTTATAGTCGGCGTGGTGACGGAGAAATGCTCTACCAGCTCTTCCGTACTGAAATAACACTGCTGTTTAACCAGTTCg                    >  2:505750/1‑149 (MQ=255)
                                                                                                                        aGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc            >  1:527533/1‑149 (MQ=255)
                                                                                                                              gcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCCGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCCGCTGTTTAACCCGTACCAAAATACCAttatt      >  2:296392/1‑146 (MQ=255)
                                                                                                                                  cGCCATGATG‑GCGCACGATTAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACCGAGAGATGCTCTCCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt  <  2:463394/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGT  >  NC_000913_3_vch_pgi/3560314‑3560591

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