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.C96.9% 77.6 / ‑3.7 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (19/12);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_0009134,296,3802.G93.5% 70.4 / ‑4.5 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (18/11);  minor base T (1/1);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCA  >  NC_000913/4296300‑4296512
                                                                                   ||                                                                                                                                    
cTTGACGCGGCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                          <  1:528225/143‑1 (MQ=1)
             aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTT‑‑ACGCCGCCTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGCTGGTTACttc                                                             >  2:45452/1‑142 (MQ=11)
             aTCAGGCCTACGCCAGACAGCGCAATAGCCTAATTTAGCGTGACTTTGTAGCTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCCACGCTTGCCGCGTCTTACCCCGCCTCTCTTAACCGTTGCTTAAtcc                                                             >  2:818429/1‑141 (MQ=11)
                aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTGTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTc                                                          >  1:712315/1‑143 (MQ=14)
                 cgCCTGCGCCAGACTGCGCAAGAGCCTGAGTGTGCGGGACTGTGTAGGCCCGATAAGGCGTTT‑‑ACGCCACCTCCGACATCAAGGCCTGATGCGACGCTTGGCACGTCTTAGCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                         <  2:65263/142‑1 (MQ=255)
                    cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAACAGTAGGTGAATGATTCAAAc                                                      >  2:71414/1‑143 (MQ=11)
                       cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAGCAGGCCTATCTTAACCGTTGGTTAATATCTCAAActct                                                   >  2:245076/1‑143 (MQ=14)
                        gCCAGACAGCGCAACAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                  <  1:546788/143‑1 (MQ=255)
                          cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                >  1:140168/1‑143 (MQ=255)
                             aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                             <  1:226684/143‑1 (MQ=255)
                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGATTAATTTTTCAAACTCTTCACACCCga                                         >  1:804782/1‑142 (MQ=17)
                                   cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTa                                       >  1:655997/1‑143 (MQ=255)
                                    aaTAGCCTGATTTAGCGTGATTTTGTAGTTCGGTTAAGGCGTTT‑‑ACTCCGCATCCGCCAGCAATGCCTTATGCGCCGCTTGCCGCGTCTGATCATGCTCATCTTCACCGTAGGTTAATTTTGCCAACTCTTCACACCCggcag                                      >  2:204830/1‑140 (MQ=9)
                                    aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCAGTTCGTTAATTTTTCAAACTCTTCACACCCCGtct                                      >  2:789645/1‑141 (MQ=17)
                                    aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                      <  2:811387/143‑1 (MQ=255)
                                        gCCTGATTCAGCGTGATTTTGCAGTTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                  <  2:367518/143‑1 (MQ=18)
                                           tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                               >  1:867041/1‑143 (MQ=255)
                                                  gCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                        <  1:760355/143‑1 (MQ=255)
                                                   cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTccttgc                       >  1:564503/1‑139 (MQ=255)
                                                     tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                     <  1:64073/143‑1 (MQ=255)
                                                     tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                     >  2:169695/1‑143 (MQ=255)
                                                     tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAACCGTTGCTTAATTTCTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAg                     >  2:591248/1‑143 (MQ=21)
                                                      gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGc                    >  2:524063/1‑143 (MQ=21)
                                                        ttttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGCTTACGTGTTCAATGTCTTAAAAAAAGTAATCAAACCCTTCCATACAGCtg                  >  1:548223/1‑142 (MQ=14)
                                                         tttGTATGTCGGAGAAGGCGTTT‑‑ACTCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                 <  2:250564/143‑1 (MQ=25)
                                                            gTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGAGTGTTAAATGGGCTGGCACTACGGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt              <  1:879156/143‑1 (MQ=18)
                                                                 tCGGATAAGGCGTTTGCG‑‑CCGCATCCGGCATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTTGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc         <  2:136055/143‑1 (MQ=25)
                                                                 tCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTCATTTTTCAAACTCTTCACACCCGGTATTAAACCCTTCCATACAGATCATGTTCAAc         >  1:435391/1‑143 (MQ=21)
                                                                   ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa       >  1:515296/1‑143 (MQ=255)
                                                                    gtgaaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg      <  2:283545/140‑1 (MQ=255)
                                                                    gATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATATTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCACACCCTTCCCTCCAGCTCGGGTTCGTCCAg      >  1:467048/1‑143 (MQ=21)
                                                                        aGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa  <  1:91005/143‑1 (MQ=255)
                                                                                   ||                                                                                                                                    
CTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCA  >  NC_000913/4296300‑4296512

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: