Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,222,617 C→T 7.3% E94K (GAA→AAA)  U712_11155 ← Uncharacterized protein ypbG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,222,6170CT 7.3% 176.4 / 3.3 55E94K (GAA→AAA) U712_11155Uncharacterized protein ypbG
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base T (2/2);  total (32/23)
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.55e-01

GGACTGACTCATTTCTCAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGAAATAAAAAAGATT  >  CP006881/2222479‑2222753
                                                                                                                                          |                                                                                                                                        
ggACTGACTCATTTCTCAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTc                                                                                                                                          >  5:220034/1‑139 (MQ=255)
  actgactCATTTCTCAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCt                                                                                                                                        <  2:27935/139‑1 (MQ=255)
   cagactCATTTCTCAGAGTGATGACGCCGGGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTc                                                                                                                                       <  7:107104/137‑1 (MQ=255)
           tttCTCAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCg                                                                                                                               >  1:113938/1‑139 (MQ=255)
              ctcAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCa                                                                                                                             <  2:154114/138‑1 (MQ=255)
                           cGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGcc                                                                                                               <  1:253834/139‑1 (MQ=255)
                           cGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGcc                                                                                                               <  3:73202/139‑1 (MQ=255)
                             ccGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATATTCATTATTCCCCCACACAAACACGATAGGGACCCCGCAATGTACGCGCCGTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCt                                                                                                             >  6:67276/1‑139 (MQ=255)
                              cGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCt                                                                                                             >  1:264432/1‑138 (MQ=255)
                                tGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGCCCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAg                                                                                                          >  2:217261/1‑139 (MQ=255)
                                       tAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATc                                                                                                   >  4:279428/1‑139 (MQ=255)
                                         aaTATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCa                                                                                                  <  3:279428/138‑1 (MQ=255)
                                             tGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCAccgc                                                                                              <  6:220034/138‑1 (MQ=255)
                                               gAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGa                                                                                           >  4:128711/1‑139 (MQ=255)
                                                    cAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCAccc                                                                                                                   <  4:42241/110‑1 (MQ=255)
                                                    cAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCAccc                                                                                                                   >  3:42241/1‑110 (MQ=255)
                                                     aGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTTCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATaa                                                                                     >  1:62402/1‑139 (MQ=255)
                                                     aGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTTCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATaa                                                                                     >  1:66761/1‑139 (MQ=255)
                                                         tATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAAcc                                                                                   >  8:282419/1‑137 (MQ=255)
                                                          aTGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCCCTGGGCACCCCGCCATCAGCTAAATCACCGCCGATAATAAACAAg                                                                                >  8:209526/1‑139 (MQ=255)
                                                             ctgcCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTgg                                                                             >  7:182594/1‑139 (MQ=255)
                                                             ctgcCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTg                                                                              <  8:82551/138‑1 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      <  6:34093/127‑1 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      >  3:4447/1‑127 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      >  5:34093/1‑127 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      <  4:4447/127‑1 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTggggg                                                                          >  7:124918/1‑139 (MQ=255)
                                                                ccTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTggggg                                                                          >  4:30690/1‑139 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCAcc                                                                                                                    <  1:97414/96‑1 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCAcc                                                                                                                    >  2:97414/1‑96 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTgggggg                                                                         >  2:15170/1‑139 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTgggggg                                                                         >  8:17282/1‑139 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTgggggg                                                                         >  7:5399/1‑139 (MQ=255)
                                                                 cTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTgggggg                                                                         >  4:23313/1‑139 (MQ=255)
                                                                                        ccACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTaaa                                                  <  1:15170/139‑1 (MQ=255)
                                                                                        ccACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTaaa                                                  <  7:17282/139‑1 (MQ=255)
                                                                                        ccACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTaaa                                                  <  3:23313/139‑1 (MQ=255)
                                                                                         cacacaAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAAc                                                                                    <  4:74995/104‑1 (MQ=255)
                                                                                         cacacaAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAAc                                                                                    >  3:74995/1‑104 (MQ=255)
                                                                                         cacacaAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTaa                                                   >  5:273662/1‑137 (MQ=255)
                                                                                         cacacaAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTaa                                                   <  6:273662/137‑1 (MQ=255)
                                                                                          acacacacACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTTCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAg                                                <  2:66761/139‑1 (MQ=255)
                                                                                          acacaAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTTCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAg                                                <  2:62402/139‑1 (MQ=255)
                                                                                                        gggACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATg                                   <  3:30690/138‑1 (MQ=255)
                                                                                                              ccGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGtt                            >  5:29787/1‑139 (MQ=255)
                                                                                                              ccGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGtt                            >  5:147591/1‑139 (MQ=255)
                                                                                                                       aCGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGGGGGGGGCATGGCTTCGCGCTGTTTCTTAAAGGTGATGATCGATGAGCCGTTTGTGAATGt                   >  2:242164/1‑139 (MQ=255)
                                                                                                                              ttttGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGaaa             >  2:67926/1‑138 (MQ=255)
                                                                                                                               tttGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAATGTCTTGATCTATGAGAAGTTTGTGAATGTcataaaa            >  1:91472/1‑132 (MQ=255)
                                                                                                                                ttGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      <  4:31523/63‑1 (MQ=255)
                                                                                                                                ttGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATtata                                                                                      >  3:31523/1‑63 (MQ=255)
                                                                                                                                      ttttCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGAAATAAAAAAGa    <  1:67926/139‑1 (MQ=255)
                                                                                                                                       tttCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGAAATAAAAAAGAt   >  4:165091/1‑139 (MQ=255)
                                                                                                                                        ttCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAAc                                                                                    <  2:193786/57‑1 (MQ=255)
                                                                                                                                        ttCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAAc                                                                                    >  1:193786/1‑57 (MQ=255)
                                                                                                                                        ttCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGAAATAAAAAAGAtt  >  1:176407/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GGACTGACTCATTTCTCAGAGTGATGACGCCGTGCGCTTTAAATATGGAATACAGCTTATGCTGCCTGACTTCATAATCATTATTCCCCCACACAAACACGATAGGGACCCCGAAATGTACGAGCCTTTTGATGTTTTCCTCGATTCTCGCACTGGGCACCCCGCCCTCAGCTAAATCACCGCCGATTATAACCAAGTGGGGGGCATGGCTTCGCGCTTTTTCTAAAAGGTCTTGATCGATGAGACGTTTGTGAATGTCTGAAATAAAAAAGATT  >  CP006881/2222479‑2222753

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