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.C97.2% 104.0 / ‑3.9 36intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (14/21);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
*NC_0009134,296,3802.G97.2% 106.5 / ‑3.8 36intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (14/21);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAG  >  NC_000913/4296313‑4296513
                                                                      ||                                                                                                                                     
aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAATAAAGACTTACTGC‑‑CCGTATCGTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                              >  1:1066325/1‑143 (MQ=17)
   aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGCCACCAATGCCGGAGGCTACGCTTGCCGCGCCTTACCCGGCCTATCTTACCCGTTGGCTactattgc                                                           >  2:640363/1‑136 (MQ=11)
      ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                        <  1:273930/143‑1 (MQ=255)
      ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                        <  1:58603/143‑1 (MQ=255)
      ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                        <  1:967166/143‑1 (MQ=255)
        tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                      <  2:632529/143‑1 (MQ=255)
           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                   <  1:359609/143‑1 (MQ=255)
            ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                  <  1:571121/143‑1 (MQ=255)
             gaGACAGCGCAATGGCCTGAGTTAGCGTGATTTTGTAGTTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                 <  2:755163/142‑1 (MQ=17)
             cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                 <  1:347784/143‑1 (MQ=255)
             cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                 >  2:595079/1‑143 (MQ=255)
                  aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTGCACAccg                                            >  1:691308/1‑142 (MQ=17)
                   gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                           <  1:605670/143‑1 (MQ=255)
                     gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGATTAATTTTTCAAAATCTTCACACCCggc                                         >  2:85043/1‑142 (MQ=17)
                      cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                        <  1:898421/143‑1 (MQ=255)
                           gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                   <  1:100089/143‑1 (MQ=255)
                           gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                   <  1:664748/143‑1 (MQ=255)
                             cttatgtAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                 <  1:337411/137‑1 (MQ=18)
                              tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                >  2:155207/1‑143 (MQ=255)
                               gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCAGGCCTATCTTAACCGATGGTTAATTTTTCAAACTCTTCAACCCCGGTATCAAAccat                               >  2:1061517/1‑141 (MQ=18)
                                 tttAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                             <  1:595079/143‑1 (MQ=255)
                                       gTGATTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                       >  2:988611/1‑143 (MQ=255)
                                          aTTTTGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                    >  2:435642/1‑143 (MQ=255)
                                           ttttGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTc                   >  2:756583/1‑143 (MQ=21)
                                            tttGTAGGTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTTTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                  <  1:834745/143‑1 (MQ=255)
                                                   gTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa           <  1:599383/143‑1 (MQ=255)
                                                   gTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa           >  2:1945/1‑143 (MQ=255)
                                                   gTCGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCCTACAGCTCTGGTTCaa           >  1:290978/1‑143 (MQ=25)
                                                     cGGATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCCTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACCCCCGGCATCGAACCCTTCCATACAACTCCGGTTctccc         >  2:300174/1‑139 (MQ=21)
                                                      ggATAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCCGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa        >  1:577757/1‑143 (MQ=255)
                                                      ggACAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa        <  1:37588/143‑1 (MQ=255)
                                                        aTAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                           >  2:973860/1‑106 (MQ=255)
                                                        aTAAGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                           <  1:973860/106‑1 (MQ=255)
                                                          aaGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc    <  2:366618/143‑1 (MQ=255)
                                                           aGGCGT‑‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa   <  1:771307/143‑1 (MQ=255)
                                                            ggCGT‑‑TTACGCCGCATCCCACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg  <  2:383651/143‑1 (MQ=255)
                                                                      ||                                                                                                                                     
ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑‑TTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAG  >  NC_000913/4296313‑4296513

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: