Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,472,568 Δ1 bp 14.8% coding (86/897 nt) ynbB → putative CDP‑diglyceride synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,472,5680A.14.8% 53.0 / 8.0 27coding (86/897 nt)ynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base . (2/2);  total (12/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 = 7.05e-01

GGTCCATTTGGCAAAAGCGATCGCGCGTTAATATTTGGTCTGTGGGGTCTGGCCGTTGCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTTCCGCGTTGGTTATATTGGGGTATTCCTCTTAACTATTTGCTTATCGGTTTTAA  >  CP009273/1472307‑1472784
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                        
ggTACATTTGTTAAAAGGGTTCGCGGGTTTATATTTTGTCTGTGGGGTCTGGCCGTTGCCATTTATCCACAGTGGATGCAGTTGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAggggg                                                                                                                                                                                                              <  2:441615/274‑1 (MQ=255)
 gTCCATTTGGCAAAAGCGATCGCGCGTTAATATTTGGTCTGTGGGGTCTGGCCGTTGCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGATTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTg                                                                                                                                                                                                            >  1:441628/1‑275 (MQ=255)
     aTTTGGCTAAAGCGTTCGCGCGTTAATATTTGGTCTGTGGGGTCTGGCCGTTGCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAgg                                                                                                                                                                                                         <  1:94179/274‑1 (MQ=255)
                   attgtgcgTTTATGTTTGGTCTGTGGGGTTTGTCCGTTTCTATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGTGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTc                                                                                                                                                                                          <  1:242962/270‑1 (MQ=255)
                         cgTTAAGATTTGGTCGGTGGGTTCTGGCCGGTGCCATTTATGCTCCGTGGGTGCAGTGGAATAATCTTTTATTGTGCATTGCTGAAATTCTGCTTCTCTGGACTGCGATTAAGCGTTGTCGGGGTGTGCTTCTTATGAGTGATGAAATATAATGCTGGAAAAATCTCTGGCAACACTTGTCGGGTTGTTAATTGTAGCAAAGCTGATAGACCGTTGTCTGCTAGGGCGGTGAACGGGGAGAAAAGGGGGTGAGGGTACATTACGTATTCGTAcc                                                                                                                                                                                     <  2:135732/274‑1 (MQ=255)
                                                    ccGTTGCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTTGAAAAATCTCTGGCAACACTTTTCGCGTTGTTACTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGTAGAGAAAAGCGGGTGAGGTTACATTACGTATTCTTACCTGTCTGGGCATCGTCATTTTTTTTTCaa                                                                                                                                                         >  1:5601/1‑275 (MQ=255)
                                                        tgtcATTTATCCCCAGTGGTTGAAGGGGAATAATCTTTTAGTGAGCATTGGTTCAATTCGGCTTCTTGTGAATGCGATTAATCGTTGTCGGAGTGTGCTTCTTAAGAGCGCAGAAATATAATGAAGGAAAAATCTCTGGCAACAATTTTAGCGTTGTTAATTTTAGCAACGCGGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGt                                                                                                                                                     <  2:5601/272‑1 (MQ=255)
                                                         gCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGNTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGTCGGTTACCGgagag                                                                                                                                                                                                                       >  1:426988/1‑208 (MQ=255)
                                                                 tCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAg                                                                                                                                             >  1:23717/1‑274 (MQ=255)
                                                                                gTGGAATAATCTTTTATGGAGCATTGCTTCAAGTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGCGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAATTGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTACAGGTCCACGTTGGAtgcc                                                                                                                             >  1:259038/1‑273 (MQ=255)
                                                                                             ttATGGCGCATTGCTTCAATTCTGCTTCTCTGAACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAGATATAATGCTGGAACAATCTCTGGCCACACTTTTCGCGTTGTTAATTTTAGCAACGATGATACACCGTTTTCTGCTATTGCGTTTACCGGAGAGAACAgg                                                                                                                                                                                                                 >  2:221070/1‑178 (MQ=255)
                                                                                                 ggAGTATTGCTTCAATTCTGCTTCTCTGGCCTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGt                                                                                                                                                                                                             <  1:221077/178‑1 (MQ=255)
                                                                                                                           tGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGAGAAAAggg                                                                                                                                                                                                                >  1:212213/1‑148 (MQ=255)
                                                                                                                               cTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGAGAAAAGGGGGTg                                                                                                                                                                                                            <  2:212206/148‑1 (MQ=255)
                                                                                                                                tgggtgttattgTTGGCTGAGTGTGCTTGTGTTGTGTGCTGATATATAATGCTGGAAAAATTTCTTGCACCACTTTTCGATTTGTTAATTTTAGCAACGCTTATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGATGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                                                                             <  2:23716/264‑1 (MQ=255)
                                                                                                                                                                    cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGgagaga                                                                                                                                                                                                                      >  2:266616/1‑102 (MQ=255)
                                                                                                                                                                    cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAtt                                                                                         >  1:480352/1‑227 (MQ=255)
                                                                                                                                                                      cTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGt                                                                                                                                                                                                             >  1:146946/1‑109 (MQ=255)
                                                                                                                                                                        gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAg                                                                                                                                                                                                                  <  1:266624/102‑1 (MQ=255)
                                                                                                                                                                        gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGaaa                                                                                     <  2:480339/227‑1 (MQ=255)
                                                                                                                                                                          aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAgg                                                                                                                                                                                                         <  2:146940/109‑1 (MQ=255)
                                                                                                                                                                              taATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgcgc                                                                                                            >  1:218894/1‑198 (MQ=255)
                                                                                                                                                                                  gCTGGAAATATCTCTGGCAAATCTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGCTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGAt                                                                                                        <  2:218887/198‑1 (MQ=255)
                                                                                                                                                                                      gAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTAtt                                                                                                                                                                                           >  1:439900/1‑111 (MQ=255)
                                                                                                                                                                                          aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTa                                                                                                                                                                                       <  2:439887/111‑1 (MQ=255)
                                                                                                                                                                                                           agtgatatgttgtttttTGCAATGTTGTTATACCGTTTTCTTCTATGGCGGTTACCG‑‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCTTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTTCCGCGTTGGTTATATTGGGGTATTCCTCTTAACTATTTGCTTATCGGTTTTaa  <  2:345616/260‑1 (MQ=255)
                                                                                                                                                                                                                              cAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGt                                                                                                                                                                                                        >  2:426090/1‑58 (MQ=255)
                                                                                                                                                                                                                              cAACGCTGATAAACCATTTTCTGCAATGATGGTTACCG‑‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGNTTATATCTGTACATTTTCCGCGTTGGTTATATTGGGGTATTCCTCTTAACTATTTGCTTATCGGtt      >  1:345628/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACa                                                                                                                                                                                                    <  1:426103/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                        
GGTCCATTTGGCAAAAGCGATCGCGCGTTAATATTTGGTCTGTGGGGTCTGGCCGTTGCCATTTATCCACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTTCCGCGTTGGTTATATTGGGGTATTCCTCTTAACTATTTGCTTATCGGTTTTAA  >  CP009273/1472307‑1472784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: