Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,289,324 G→T Q56K (CAG→AAG)  ybbP ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,289,3240GT96.7% 91.5 / ‑4.7 30Q56K (CAG→AAG) ybbPhypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (19/10);  minor base C (0/1);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAATTTT  >  CP000730/2289191‑2289450
                                                                                                                                     |                                                                                                                              
cTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCa                                                                                                                         >  2:57456/1‑141 (MQ=255)
     tctACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAAc                                                                                                                    >  2:377020/1‑141 (MQ=255)
     tctACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAAc                                                                                                                    >  2:399175/1‑141 (MQ=255)
               tgttCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAat                                                                                                          >  2:202753/1‑141 (MQ=255)
                   cTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTACTATcc                                                                                                      >  2:354482/1‑141 (MQ=255)
                   cTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATcc                                                                                                      >  2:404528/1‑141 (MQ=255)
                             tCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTaa                                                                                            >  1:107318/1‑141 (MQ=255)
                             tCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTaa                                                                                            >  1:206708/1‑141 (MQ=255)
                               tAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAAtt                                                                                          >  1:262311/1‑141 (MQ=255)
                               tAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAAtt                                                                                          >  1:262524/1‑141 (MQ=255)
                                 aTTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGt                                                                                        <  1:207961/141‑1 (MQ=255)
                                     cTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCg                                                                                    <  1:216639/141‑1 (MQ=255)
                                                          aaGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCg                                                               <  2:375507/141‑1 (MQ=255)
                                                               agcaCCCCCCATGGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGtgatg                                                          <  2:358854/138‑1 (MQ=255)
                                                                     ccccATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGATCGTGATGagaaga                                                    >  2:373315/1‑141 (MQ=255)
                                                                     ccccATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGagaaga                                                    >  2:44871/1‑141 (MQ=255)
                                                                            gAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTa                                             <  2:23784/141‑1 (MQ=255)
                                                                                aaCGATATCGAATAATTTAGATGTTGCAGCCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAAACCGCGATGAGAAGATAAAGTACATa                                         >  2:176373/1‑141 (MQ=255)
                                                                                      atCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACAAAAc                                   >  2:57055/1‑141 (MQ=255)
                                                                                       tCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACt                                  <  2:138576/141‑1 (MQ=255)
                                                                                       atgAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACt                                  <  2:212751/139‑1 (MQ=255)
                                                                                           ataatTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTAtta                              <  2:262311/141‑3 (MQ=255)
                                                                                           ataatTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCa                              <  2:262524/141‑1 (MQ=255)
                                                                                                   gATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTAttaataagtta                      >  1:399848/1‑131 (MQ=255)
                                                                                                            gTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTc             >  2:338998/1‑141 (MQ=255)
                                                                                                                  ttCAATATCATACTTACCCCATGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACa       <  1:114560/141‑1 (MQ=255)
                                                                                                                    cAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAAt     >  1:436845/1‑141 (MQ=255)
                                                                                                                      atatCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAAttt   >  2:78794/1‑141 (MQ=255)
                                                                                                                      atatCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAAttt   >  2:237705/1‑141 (MQ=255)
                                                                                                                       tatCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAAtttt  <  1:256361/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                              
CTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAATTTT  >  CP000730/2289191‑2289450

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