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.C87.9% 75.3 / 6.8 34intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (14/15);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G88.2% 79.0 / 7.2 34intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (14/16);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑A‑‑CCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTG  >  NC_000913/4296304‑4296510
                                                                                  ||                                                                                                                                    
aCGCTTGTTATCAGGTCTACGACAGACAGCGCAATACCATGATTTAGCACGATTTCGTAGGTCGGATAA‑GGCGTC‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                    <  2:613408/143‑1 (MQ=11)
aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTAACCGtt                                                                    <  2:806609/143‑1 (MQ=255)
            aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                        >  1:828306/1‑143 (MQ=255)
                     ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTCT‑GGCGAT‑‑G‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                               <  1:745936/143‑1 (MQ=17)
                     ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGAATTAAT‑‑A‑‑CCGCG‑TTAGGCATATA‑GAAGTAATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTGCCCGTTGGTTAATTTTTCAAACTCTTc                                               >  1:221544/1‑143 (MQ=12)
                      cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTGTCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTca                                              >  1:281196/1‑143 (MQ=17)
                      cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTCACCGTTGGTTAATTTTTCATCCTCTTcc                                              >  1:465233/1‑142 (MQ=14)
                      cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATCTTCCAAACTCTTca                                              >  2:372674/1‑143 (MQ=14)
                          cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ATGCCGCA‑TCCGACATCAA‑CGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                          <  1:618216/143‑1 (MQ=14)
                            gcggAATAGGCTGATCTCGCGTGTTGTTGTAGGTCGGAGAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                        <  2:455323/143‑1 (MQ=17)
                             cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATTTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                       >  1:730184/1‑143 (MQ=255)
                                     ccTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAAc                               >  2:159347/1‑143 (MQ=18)
                                     ccTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                               <  1:892479/143‑1 (MQ=255)
                                     ccTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGCTATCAAACTCTTCACACCCGGCATCAAAc                               >  2:628230/1‑143 (MQ=17)
                                     ccTGATTTAGCGCGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                               <  1:315184/143‑1 (MQ=255)
                                         aTTTGGCCTCATCTTCTAGGTCGGATAC‑GCCGTT‑‑T‑‑A‑‑CGCCACTTCGACAAGAACCTCCTGCTGCGAAGATTGGCTCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                           <  2:202059/143‑1 (MQ=255)
                                          tttAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                          <  1:463395/143‑1 (MQ=255)
                                          ttgaGCGTGATTTTGTGGGCCCGATAAGGGCGTTA‑C‑‑G‑‑CGGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                          <  2:448315/140‑1 (MQ=21)
                                           ttAGTGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCAACCTCTTCACGCCCGGTATCAAACCCTTcc                         >  2:327633/1‑143 (MQ=18)
                                           ttAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑GCCGACATCAA‑TGCCTGAAGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                         <  1:523783/143‑1 (MQ=21)
                                                 tGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTCATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAtcccg                   >  1:507422/1‑139 (MQ=21)
                                                     tttGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa               <  1:167695/143‑1 (MQ=255)
                                                     tttGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGATTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAATGTTTCAAACTCTTCAAAACCGGTATCAAACCCTTGAATACAGcata               >  2:264116/1‑140 (MQ=18)
                                                      ttGTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGATTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCCACCCGGTATCAAACCCTTCCATACAGCTCAg              >  2:627271/1‑143 (MQ=25)
                                                        gTAGGTCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt            <  1:785133/143‑1 (MQ=255)
                                                        gTAGGTCGGATAA‑CGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt            <  2:709145/143‑1 (MQ=255)
                                                             tCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTGTCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc       <  2:281196/143‑1 (MQ=255)
                                                             tCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc       >  2:443181/1‑143 (MQ=255)
                                                             tCGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCATTCCATACAGCTCCGGTTAAAc       >  2:443166/1‑143 (MQ=25)
                                                             gcGGATAA‑GGCGTT‑‑TTCG‑‑CCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc       <  2:442816/142‑1 (MQ=255)
                                                             ccGGATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc       <  1:345042/142‑1 (MQ=255)
                                                               ggATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACACCCGGTATCAAACCCTTCCATACAACTAAGAATCAAACa     >  1:371932/1‑143 (MQ=21)
                                                               ggATAA‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa     <  1:308842/143‑1 (MQ=255)
                                                                  tca‑GGCGTT‑‑T‑‑ACGCCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg  <  1:65506/141‑1 (MQ=255)
                                                                                  ||                                                                                                                                    
ACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAA‑GGCGTT‑‑T‑‑A‑‑CCGCA‑TCCGACATCAA‑TGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTG  >  NC_000913/4296304‑4296510

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: ‑