Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.96.8% 125.6 / ‑3.2 31intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (13/17);  minor base C (1/0);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_pae_pgi/4234520‑4234810
                                                                                                                                                   |                                                                                                                                                
tCTATATACGCCGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                >  2:154041/1‑147 (MQ=255)
    tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcg                                                                                                                                           <  1:154041/149‑1 (MQ=255)
      taCGTAGTAGAGATGGACCCCAGGGGCGTTTTATGCTTTCTACGGGATGTGGGCATCGCCTCTGTAAAACTGGCAGCGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                         <  2:116863/149‑1 (MQ=255)
         gtaGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTGACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                      <  2:106318/147‑1 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                  >  1:328674/1‑138 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                  <  2:328674/138‑1 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                      >  2:289400/1‑149 (MQ=255)
             tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                >  1:237021/1‑134 (MQ=255)
             tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                <  2:237021/136‑3 (MQ=255)
               gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAACCTGGCAGAGATGAGGGGCGTAAATCCCCTGGCCACCCCTAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGGTATTTAg                                                                                                                                >  2:333146/1‑149 (MQ=255)
                  aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGcc                                                                                                                             >  1:281608/1‑149 (MQ=255)
                       ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATa                                                                                                                        >  2:380249/1‑149 (MQ=255)
                           aaGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                    >  1:9804/1‑149 (MQ=255)
                                cGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                               <  2:143741/149‑1 (MQ=255)
                                    ttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCTAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCcgcc                                                                                                           <  1:102114/149‑1 (MQ=255)
                                     tATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACCCACCTAATGAGGATGCCTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTGGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGGTAAcgccg                                                                                                          >  1:98106/1‑149 (MQ=255)
                                           ttttACGAGATGTGGGCATCGCTCCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  1:12798/149‑1 (MQ=255)
                                               aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAAACCCCTCGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                <  1:95418/149‑1 (MQ=255)
                                                cGAGATGTGGGCATCGTCTTTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               <  1:352172/149‑1 (MQ=255)
                                                cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGACCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               >  2:68190/1‑149 (MQ=255)
                                                               gCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                                <  1:225529/149‑1 (MQ=255)
                                                                                agggaTTAGCGGCGTAAATCCCCTGGCCAACACGAGAATAGATGCATTGAAGAGGCGCTTGCACTGATCATCGTCGAGATTTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                               <  2:98106/146‑1 (MQ=255)
                                                                                            gCGTAAATACCCGGGCCACCCCGAGAATAGATGCAGTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGCATAAGGCGTTCACGCCGCATGTGGAAACCGATGCCTTATAAGACGCGCCCGCGTCGCATCAGGcatcggtt                                                    <  2:69090/149‑8 (MQ=255)
                                                                                                       cTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                         <  1:380249/149‑1 (MQ=255)
                                                                                                        tgcccACACCAAGAATAGAGCGATTTAAGAGGCGGTAGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        <  1:131442/146‑1 (MQ=21)
                                                                                                        tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCCTCGGCGAAATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        >  1:388804/1‑149 (MQ=25)
                                                                                                        tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        <  2:9804/149‑1 (MQ=255)
                                                                                                         ggCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAGAAGCCGTTCACGCCGCATACGGCAACCGTTGCCTGATGCACCGCGGTCGCGTCTTGATAGGCCTACCGGTTCATGCCGatat                                       >  2:122041/1‑149 (MQ=255)
                                                                                                          gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTCATCATCGTCGATATATAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTCGCGTCTTATCAGTCTTACAGGTCTATGCCGATATg                                      >  2:391557/1‑149 (MQ=37)
                                                                                                                      aTAGATGCAGTGAAGAGGCGCTTGCAGTGAGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATTCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                          <  1:46503/149‑1 (MQ=255)
                                                                                                                          aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                      >  1:205904/1‑149 (MQ=255)
                                                                                                                          aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                      >  2:401088/1‑149 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACTGATGCCTGATGCGACGCGGTCGCGTGTTGTCAGGCGTACAGGTCGATGCAGATATGTGTATCGGATGCGGCAAGTa                  >  1:121846/1‑149 (MQ=25)
                                                                                                                               ttGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                 <  1:68190/149‑1 (MQ=255)
                                                                                                                                 gAAGAGGCGCTTGCAGTGATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCTACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAGCGTATTCGGCAATTattc               <  2:246030/149‑4 (MQ=25)
                                                                                                                                              cAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAGATCGTTTTCGGCACTTAATACATAGCACgattg  >  1:75889/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_pae_pgi/4234520‑4234810

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: