Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,381:1 +GC 100% 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.C91.0% 80.2 / ‑2.5 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (11/15);  minor base T (2/1);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G96.7% 88.8 / ‑2.7 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (13/16);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAC  >  NC_000913/4296345‑4296469
                                    ||                                                                                           
aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                             <  1:731416/100‑1 (MQ=255)
 tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                               <  2:628092/97‑1 (MQ=255)
  ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGC‑C‑GCATCCGACATCAACGCCTGATGCGATGCTTAACGCGTCTTATCAGGCCTATGTTagacgg                            >  2:1364493/1‑94 (MQ=1)
   tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                          <  1:1181361/100‑1 (MQ=255)
    aGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCATGACGCGTCTTATCAGGCCTATCTTAACCGTTg                          >  2:1181357/1‑99 (MQ=11)
    aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                         <  2:518492/100‑1 (MQ=255)
     gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCACTGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTATCTTAAACGTTGGt                        >  1:1081092/1‑100 (MQ=11)
     gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTATCTTAACCGTTGGt                        <  2:1297007/100‑1 (MQ=14)
     gCGTGATTTTGTAGGTCGGATAAGGCGATTACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTACCGCGTCCTATCAGGCCTATCTTAACCGTTg                          >  1:1754255/1‑98 (MQ=11)
         gATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                    <  2:711419/100‑1 (MQ=255)
         gATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                    >  2:1617577/1‑100 (MQ=255)
          aTTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                   <  1:510092/100‑1 (MQ=255)
            tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                 <  2:1857788/100‑1 (MQ=255)
             ttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                <  1:1661712/100‑1 (MQ=255)
              tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc               >  1:1147910/1‑100 (MQ=255)
              tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc               <  2:2039014/100‑1 (MQ=255)
               gTAGGTCGTATAAGGCGTTTA‑‑CGCTGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCAACGGTTaagataggc              <  2:1023837/100‑8 (MQ=11)
                tAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa             <  1:87153/100‑1 (MQ=255)
                  ggTCTGATAAGGGGTTTATGC‑C‑GCATCCCACATCTATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc           <  2:1845207/100‑1 (MQ=18)
                  ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc           <  2:1728003/100‑1 (MQ=255)
                  ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc           >  2:257651/1‑100 (MQ=255)
                  ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc           >  1:1224683/1‑100 (MQ=255)
                   gTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct          >  2:1742855/1‑100 (MQ=255)
                    tCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc         >  2:2207420/1‑100 (MQ=255)
                    tCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc         <  1:1291526/100‑1 (MQ=255)
                     cGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct        <  2:32069/100‑1 (MQ=255)
                      ggATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa            >  1:714078/1‑95 (MQ=255)
                        aTAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca     <  1:1662176/100‑1 (MQ=255)
                         tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac    >  1:1615004/1‑100 (MQ=255)
                           aGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac  >  1:1541115/1‑100 (MQ=255)
                                    ||                                                                                           
ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAC  >  NC_000913/4296345‑4296469

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑