Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 4,296,411 +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_000913_3_vch_tpiA4,296,4101.C90.0% 72.3 / 2.9 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (11/16);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_vch_tpiA4,296,4102.G90.6% 83.6 / 2.8 32intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (12/17);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_vch_tpiA/4296325‑4296548
                                                                                          ||                                                                                                                                          
ggACGCTTACTGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGACAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                               <  1:126121/148‑1 (MQ=11)
   cGCTGAACGCGGCTTATGAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCATT‑‑TCCT‑‑CCCCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                            <  2:32743/149‑1 (MQ=11)
    gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                           <  1:712239/149‑1 (MQ=255)
      ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                         >  1:421918/1‑149 (MQ=255)
         aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGCCGGTTGCCGCGTCTTGTCCGGCCTATCTTAACCGTTGGTTaa                                                                      >  2:288470/1‑149 (MQ=11)
              tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                 <  1:419255/149‑1 (MQ=255)
                ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                               >  2:379148/1‑149 (MQ=255)
                ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACACTTGCCGCAACTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                               >  2:118242/1‑149 (MQ=14)
                  aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                             <  2:421918/149‑1 (MQ=255)
                       gCCTACGCCAGATAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  2:223912/149‑1 (MQ=17)
                       gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  1:201970/149‑1 (MQ=255)
                         cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                      >  1:683375/1‑149 (MQ=255)
                          tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACCCTTccc                                                     >  1:437105/1‑147 (MQ=17)
                               cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                >  2:651781/1‑149 (MQ=255)
                                       gcAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                        <  1:314120/149‑1 (MQ=255)
                                         aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                      >  1:509609/1‑149 (MQ=255)
                                           tggcgtGATTTGGCGTGATTTTGTTGGTTGGAGCAGTTGTT‑‑TTCG‑‑CCGCCGCTGTTATCATTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                    <  2:292446/144‑1 (MQ=18)
                                             ggcTGCTTTAGGGTGATTTTGTAGGTCGGGTAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  2:453476/147‑1 (MQ=21)
                                             gCCTGATTTAGCGTGATTTTGTAGGTCTGATAATGCGTTTCT‑‑G‑‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTGTCAAACCCTTc                                  <  2:225056/149‑1 (MQ=18)
                                              ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTCATTTTCCCAACTCTTCACACCCGGTATCAAACCCTTcc                                 >  2:590598/1‑149 (MQ=18)
                                               cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                 <  2:289688/148‑1 (MQ=255)
                                               cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                 >  1:289688/1‑148 (MQ=255)
                                               cTGATTTACCTTGATTTTGTTGGCCTGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTGCACACCCGGTATCAAACCCTTCCa                                <  1:40555/149‑1 (MQ=18)
                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                              <  1:69896/149‑1 (MQ=255)
                                                   tttAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                            <  1:427322/149‑1 (MQ=255)
                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                        <  2:692074/149‑1 (MQ=255)
                                                            aTTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        >  1:48233/1‑112 (MQ=14)
                                                            aTTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  2:48233/112‑1 (MQ=14)
                                                            aTTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                   <  1:435538/149‑1 (MQ=255)
                                                                tGTAGGTCGGATAAGGCGTT‑‑C‑‑ACGCCCCATCCAACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc               <  1:466770/149‑1 (MQ=25)
                                                                  tAGGTGGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa             <  2:160628/149‑1 (MQ=255)
                                                                             aGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGGTTCAACCAGTGCAGCGCCt  >  2:403443/1‑149 (MQ=255)
                                                                                          ||                                                                                                                                          
CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_vch_tpiA/4296325‑4296548

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: