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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C93.8% 72.0 / ‑0.8 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (12/18);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_0009134,296,3802.G93.5% 70.4 / ‑0.9 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (12/17);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296305‑4296511
                                                                              ||                                                                                                                                     
cgcgTTTTATCAGTCCTACGCCAGGGAGCGCACTGGGCTCATCTAGTGCGATTATGTAGGTCGCCTAAGGGGTTT‑‑T‑‑CGTCGCATCTGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                    <  2:571946/143‑1 (MQ=2)
cgcgTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                    <  1:315501/143‑1 (MQ=255)
    tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                <  1:681577/143‑1 (MQ=255)
        aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                            >  2:362734/1‑143 (MQ=255)
          cggggcTACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                          <  2:708205/138‑1 (MQ=14)
          cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                          <  2:657029/143‑1 (MQ=255)
              ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCcaa                                                      >  1:257563/1‑143 (MQ=255)
              ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                      <  1:591247/143‑1 (MQ=255)
                 aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACCTCAATGCCTGATGCGACGATTGACGCGAATTATCAGACCAATCTTAACCGTTGGTTAATTTTTCCAActc                                                   >  1:11025/1‑143 (MQ=14)
                 aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTAACCGATGGTCAATTTTTCacaatc                                                   >  1:427296/1‑138 (MQ=14)
                 aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTGGGTTAATCGTTCAAActc                                                   >  2:765972/1‑143 (MQ=14)
                   gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTACCAGCCCTATCTTAACCGTTGGTCAATTTTCCAAACCCtt                                                 >  2:803498/1‑143 (MQ=14)
                    ccAGACAGCGCAATAGCCCGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                <  1:270465/143‑1 (MQ=255)
                     cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                               <  1:643757/143‑1 (MQ=255)
                       gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATCCGACGCTTGCCGCGCCTTCTCCGGCCTATCTTAACCGCTTGTTAATTATTCAAACTCATcaca                                             >  1:220289/1‑143 (MQ=14)
                         cAGCGCAATAGCCTGATTTAGCGTGCTTTGGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCCGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                           <  2:252648/143‑1 (MQ=17)
                             gcgatggCCTTATTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                       <  2:805245/137‑1 (MQ=18)
                              cAATGGCCTGATTGACGGTGATTTTGGAGGCCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                      <  2:55639/143‑1 (MQ=18)
                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATGTTTCAAACTCTTCACACGCGGTATc                                    >  2:225986/1‑143 (MQ=17)
                                   gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCGG‑‑C‑‑CGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                 <  1:144959/143‑1 (MQ=18)
                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGTTTCCTACATCAATGCCTAATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                <  1:714666/143‑1 (MQ=18)
                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCCCACTCTCCACACCCGGTATCCAAc                                >  2:524996/1‑143 (MQ=17)
                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                <  1:475333/143‑1 (MQ=255)
                                         tctAGCGTGATTTCGTAGGTCGGATAAGGCGTTT‑‑ACTC‑‑CGCATCCGACATCAATGCCTGATGCGCTAAATATCTCTCCTTGTCAGGCTTACCCGAATCGTTTGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                           <  1:364670/141‑1 (MQ=12)
                                           tAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                         <  1:246640/143‑1 (MQ=255)
                                              cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTCTTCAAACCCTTCCCACCCGGCCTCACACCCTTCCAtcc                      >  2:185593/1‑141 (MQ=17)
                                              cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCTTTACATAc                      >  2:175475/1‑143 (MQ=21)
                                                   ttttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                 >  2:487362/1‑143 (MQ=255)
                                                               gATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg     <  1:451652/143‑1 (MQ=255)
                                                                 tAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg   <  1:118255/143‑1 (MQ=255)
                                                                 tAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg   <  1:362734/143‑1 (MQ=255)
                                                                  agggCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCGGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc  <  2:156480/141‑1 (MQ=37)
                                                                              ||                                                                                                                                     
CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296305‑4296511

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: