Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C95.0% 104.2 / ‑1.6 41intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (18/20);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G95.1% 105.7 / ‑1.6 41intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (19/20);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAG  >  NC_000913/4296284‑4296508
                                                                                                    ||                                                                                                                                  
aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTGGGCTAGTAACGCAACT‑‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                      >  1:1308545/1‑143 (MQ=1)
                 ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGAAATAAATGCCTAATACGACGCTTGCCGCGTCTTATCAGGCATATCTTAAcc                                                                     >  2:1034357/1‑143 (MQ=11)
                  tAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                    >  1:846434/1‑143 (MQ=255)
                    aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                  <  1:1313242/143‑1 (MQ=255)
                    aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                  <  1:167308/143‑1 (MQ=255)
                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGGGTCTTATCATACCTATCTTAACCTTTGGttc                                                             >  1:326142/1‑142 (MQ=11)
                             aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTATCTTACACGTTGGTTAAtct                                                         >  2:730192/1‑141 (MQ=11)
                              tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                        <  1:424457/143‑1 (MQ=255)
                              tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGGCGG‑‑A‑TAAGGCGTTTA‑‑CGCGGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGACGTATCTTAACCGTTGGTTCGtttt                                                        >  1:840697/1‑143 (MQ=11)
                              tCAGGCCTACGCCAGACAGCGCAATAGCCGTATTTCGCGTGATTTTGTATGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                        <  2:1177329/143‑1 (MQ=14)
                               cAGGCCTACGCCAGACAGCGCAATCGCCTGATTTCGCGTGATTCTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                       <  2:459197/143‑1 (MQ=14)
                               cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTAtttttt                                                       >  1:232992/1‑143 (MQ=14)
                               cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                       >  2:409633/1‑143 (MQ=255)
                               cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                       <  1:494407/143‑1 (MQ=255)
                                aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                      <  1:1306817/143‑1 (MQ=255)
                                  gCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                    >  1:604092/1‑143 (MQ=255)
                                  gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                    <  1:811219/143‑1 (MQ=255)
                                  gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACACCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAACGTTGGTTAATTTTTCaa                                                    >  2:554916/1‑143 (MQ=14)
                                        ggcAGACCGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTGTACGC‑C‑GGATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                              <  1:660830/141‑1 (MQ=17)
                                           aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTGAAACTGTTcac                                           >  1:1151386/1‑143 (MQ=14)
                                            gACAGCGCAATAGCCTGATTTAGCGTGCTTTGGTAGGTCGG‑‑A‑TAAGGCGTTTATGC‑T‑GGAACAGACACAGCTGGATGATGCGACGCGTGCCGCGTCGTATCAGGCCTAGCTTAACCGGTGGTTAATGTTTCAAACTCCTcgca                                          >  2:279027/1‑143 (MQ=9)
                                               aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                       >  1:1230279/1‑143 (MQ=255)
                                                        gCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                              <  1:279514/143‑1 (MQ=255)
                                                        gCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                              <  1:277777/143‑1 (MQ=255)
                                                              tttAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                        <  1:860752/143‑1 (MQ=255)
                                                                 aGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTCTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACAACACACACCCGGTATCAAGCCCTTCCAt                     >  2:1207022/1‑143 (MQ=18)
                                                                  gCGTGAGTTTGTAGGCCGG‑‑A‑CAAGGCGTCTACGC‑C‑GCATCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                    <  1:464150/143‑1 (MQ=21)
                                                                     tGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATACCCCCCTtccatccac                 >  1:902672/1‑142 (MQ=18)
                                                                        ttttGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGTCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc              >  2:984555/1‑143 (MQ=255)
                                                                        ttttGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCCAACCCTTCCATACAGCTc              >  1:517454/1‑143 (MQ=21)
                                                                          ttGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTCCTTTTTCAAACTCTTCACACCCGGTATCAACCACTTCCATACAGcacaa            >  1:139699/1‑142 (MQ=18)
                                                                          ttGTAGGTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg            <  2:111667/143‑1 (MQ=255)
                                                                           ctgtggTCGG‑‑T‑TATGGCGTTTACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg           <  2:775940/139‑1 (MQ=25)
                                                                              aGGTCGG‑‑A‑TAAGGCGTTGACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc        <  1:1318562/143‑1 (MQ=255)
                                                                               ggTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa       >  2:133921/1‑143 (MQ=255)
                                                                               ggTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa       <  2:604092/143‑1 (MQ=255)
                                                                                gTCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCGGGTATCAAACCCTTCCGCACAGAACAGGATCac      >  2:442079/1‑142 (MQ=21)
                                                                                 tCGG‑‑A‑TATGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc     >  2:1178011/1‑143 (MQ=255)
                                                                                 tCGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc     <  1:828564/143‑1 (MQ=255)
                                                                                 gcGG‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc     <  2:509900/142‑1 (MQ=255)
                                                                                    g‑‑A‑TAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTAAACCGTGTGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg  <  1:185304/143‑1 (MQ=25)
                                                                                                    ||                                                                                                                                  
AACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGG‑‑A‑TAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAG  >  NC_000913/4296284‑4296508

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