Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,821 T→G 12.9% intergenic (+150/+56) loiP → / ← speB Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,8210TG12.9% 88.3 / 7.2 31intergenic (+150/+56)loiP/speBPhe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base G (2/2);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

GACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGT  >  NC_000913/3082721‑3082962
                                                                                                    |                                                                                                                                             
gACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgccttat                                                                                                                                     >  1:286442/1‑107 (MQ=37)
gACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgccttat                                                                                                                                     <  2:286442/111‑5 (MQ=37)
 aCATGACTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc                                                                                                                       <  2:97604/124‑1 (MQ=255)
 aCATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc                                                                                                                       >  1:97604/1‑124 (MQ=255)
  cATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCa                                                                                                                                             <  1:310664/101‑1 (MQ=255)
  cATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCa                                                                                                                                             >  2:310664/1‑101 (MQ=255)
             caTTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt                                                                                                                  <  2:471964/117‑1 (MQ=255)
             caTTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt                                                                                                                  >  1:471964/1‑117 (MQ=255)
                        ttGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgcctta                                                                                                                                      >  2:775055/1‑83 (MQ=21)
                        ttGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCATACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgcctta                                                                                                                                      <  1:775055/86‑4 (MQ=21)
                                             atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                >  3:518272/1‑87 (MQ=255)
                                             atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                <  4:518272/87‑1 (MQ=255)
                                             atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCAtt                                                   >  2:945979/1‑148 (MQ=255)
                                                  ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat                                                             <  3:671290/133‑1 (MQ=255)
                                                  ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat                                                             >  4:671290/1‑133 (MQ=255)
                                                    cgcATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAata                                                            >  3:502437/1‑132 (MQ=255)
                                                    cgcATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAata                                                            <  4:502437/132‑1 (MQ=255)
                                                      cATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgcc                                          <  1:945979/148‑1 (MQ=255)
                                                                          tGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGc                      >  3:300992/1‑148 (MQ=255)
                                                                                 cGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGAtt               <  4:300992/148‑1 (MQ=255)
                                                                                       tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat                                                             <  4:876377/96‑1 (MQ=255)
                                                                                       tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat                                                             >  3:876377/1‑96 (MQ=255)
                                                                                       tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTc                                                 <  2:819405/108‑1 (MQ=255)
                                                                                       tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTc                                                 >  1:819405/1‑108 (MQ=255)
                                                                                           aTAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAAtata                                                          >  2:562831/1‑95 (MQ=255)
                                                                                           aTAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAAtata                                                          <  1:562831/95‑1 (MQ=255)
                                                                                            tAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATc    >  3:288406/1‑148 (MQ=255)
                                                                                                 cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGc                                                                                <  4:72603/67‑1 (MQ=255)
                                                                                                 cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGc                                                                                >  3:72603/1‑67 (MQ=255)
                                                                                                 cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                  <  2:706588/81‑1 (MQ=255)
                                                                                                 cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                  >  1:706588/1‑81 (MQ=255)
                                                                                                   ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTgg                                                                                           <  1:301093/54‑1 (MQ=255)
                                                                                                   ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTgg                                                                                           >  2:301093/1‑54 (MQ=255)
                                                                                                   ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGt  <  2:424144/143‑1 (MQ=255)
                                                                                                   ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGt  >  1:424144/1‑143 (MQ=255)
                                                                                                    tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc             >  2:57416/1‑131 (MQ=255)
                                                                                                    tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc             <  1:57416/131‑1 (MQ=255)
                                                                                                    |                                                                                                                                             
GACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGT  >  NC_000913/3082721‑3082962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: