Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,062,550 G→A 16.7% intergenic (‑3/+18) yccJ ← / ← wrbA uncharacterized protein/NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,062,5500GA16.7% 48.0 / 4.2 24intergenic (‑3/+18)yccJ/wrbAuncharacterized protein/NAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (2/2);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

TCCCAAATCTTTTCCGCCATTTTTTCGTCGTACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCGTCACCGCCTGCGATGGTGGTTGCGCCGTA  >  CP009273/1062417‑1062687
                                                                                                                                     |                                                                                                                                         
tcCCAAATCTTTTCCGCCATTTTTTCGTCGTACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGaa                                                                                                                            >  1:129455/1‑149 (MQ=255)
       ccctctCCGCCAGTTGTTCGTGTGTCCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTCCGTTTTTGTTGGGTGCATGCGTATCCTCCTGTTGAAGa                                                                                                                          <  2:345198/139‑1 (MQ=255)
               gCCATTTTTTCGTCGTACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTGGCTGCTTGAGTTGGCATGCGTATCCTCCTGTTGAGGATTTGGCGTTAAGt                                                                                                             >  1:250330/1‑149 (MQ=255)
                  aTTTTCTCGTCCTACCCGGCGCCTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCCTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTa                                                                                                          <  1:162222/149‑1 (MQ=255)
                   ttttttCGTCGTACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTaa                                                                                                         <  1:143426/149‑1 (MQ=255)
                     ttttCGTCGAACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTACTGCGCAGAGATGCCCATTCAACGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACt                                                                                                       <  1:396585/149‑1 (MQ=255)
                           tcgtACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTg                                                                                                      >  2:245243/1‑144 (MQ=255)
                                 ccGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGTAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCg                                                                                           >  1:256261/1‑149 (MQ=255)
                                         ttCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGTGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCATATCCTCCTGTTGAAGATTGGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTc                                                                                   >  1:1202/1‑149 (MQ=255)
                                         ttCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGTGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCATATCCTCCTGTTGAAGATTAGCCGTTACGTTTGACTTCCAGACCTGCGACATATTc                                                                                   >  1:107827/1‑149 (MQ=255)
                                         ttCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCGTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACa                                                                                        >  2:119930/1‑144 (MQ=255)
                                              aaaTGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTc                                                                                   <  2:35812/144‑1 (MQ=255)
                                                   gCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGCGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCCGACCTGCGACATATTCCCCTTGATac                                                                         >  1:43070/1‑148 (MQ=255)
                                                    ccTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAAc                                                                        >  1:110902/1‑149 (MQ=255)
                                                     cTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGATTCGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGa                                                                            <  2:20568/144‑1 (MQ=255)
                                                      tCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGa                                                                      <  1:164243/149‑1 (MQ=255)
                                                                    cGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAGCGAGCAATAgcc                                                             >  2:39223/1‑142 (MQ=255)
                                                                                      tGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCg                                      <  1:282340/149‑1 (MQ=255)
                                                                                      tGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCg                                      <  1:220194/149‑1 (MQ=255)
                                                                                           aTTCACCGACGTGATGAGCTTTCGCATCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAAc                                 >  1:238174/1‑149 (MQ=255)
                                                                                                  gACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCg                               <  2:69374/144‑1 (MQ=255)
                                                                                                                    ttcttGGGTTGGCATGCATATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCGTCACCGCCTGCGAtggtg             <  2:982/144‑1 (MQ=255)
                                                                                                                    gtcttGAGTTGGCATGCATATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCGTCACCGCCTGCGAtggtg             <  2:96175/143‑1 (MQ=255)
                                                                                                                               gCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCGTCACCGCCTGCGATGGTGGTTGCGCCGTa  <  2:114831/144‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
TCCCAAATCTTTTCCGCCATTTTTTCGTCGTACCCGGCGACTTCGAAAATGGCCTCGGCTATTTCCGGCGACGTATTGCGCAGAGATGCCCATTCACCGACGTGATGAGCTTTCGCTTCTTGAGTTGGCATGCGTATCCTCCTGTTGAAGATTAGCCGTTAAGTTTAACTGCCAGACCTGCGACATATTCCCCTTGATAACGAGCAATAGACAGTTCTTCCTGGCTGGGCTGGCGTGAACCGTCACCGCCTGCGATGGTGGTTGCGCCGTA  >  CP009273/1062417‑1062687

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