Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,628 A→G 11.1% I237I (ATT→ATC PP_0636 ← cold shock DNA‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947743,6280AG11.1% 95.2 / 3.4 36I237I (ATT→ATCPP_0636cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (12/20);  new base G (2/2);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTGGGCGTTTTCATTTTCGATGCGAGCTAA  >  NC_002947/743483‑743774
                                                                                                                                                 |                                                                                                                                                  
ccGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATgg                                                                                                                                                <  1:215063/150‑1 (MQ=255)
   gggAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGAATCAAGCTAGCTTGACGGTCACGCTAAGCCGAGTGAGGCATTGGATCCTGCTGGCGACGCATTTGCTTAGCTAAGTCCACAATGCCGa                                                                                                                                             >  2:340494/1‑150 (MQ=255)
         ctactTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCa                                                                                                                                       >  2:705337/1‑150 (MQ=255)
                          gtgGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAAtgtg                                                                                                                      <  2:589257/150‑1 (MQ=255)
                           tgGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTgg                                                                                                                     <  1:50618/150‑1 (MQ=255)
                                       tGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGTTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATTGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCa                                                                                                         <  2:65765/150‑1 (MQ=255)
                                       tGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGGTGGGCGGTCACCCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCa                                                                                                         <  2:188840/150‑1 (MQ=255)
                                          aGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGa                                                                                                      <  1:436837/150‑1 (MQ=255)
                                             ggAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCa                                                                                                   <  2:780791/150‑1 (MQ=255)
                                                gCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGaa                                                                                                <  1:753044/150‑1 (MQ=255)
                                                           ggCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCgtgt                                                                                     <  2:397198/150‑1 (MQ=255)
                                                            gCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACGATGGCGAAGCCCAGGATTTTCATGa                                                                                                                           >  1:408829/1‑111 (MQ=255)
                                                            gCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACGATGGCGAAGCCCAGGATTTTCATGa                                                                                                                           >  1:408993/1‑111 (MQ=255)
                                                            gCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACGATGGCGAAGCCCAGGATTTTCATGa                                                                                                                           <  2:408993/111‑1 (MQ=255)
                                                            gCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACGATGGCGAAGCCCAGGATTTTCATGa                                                                                                                           <  2:408829/111‑1 (MQ=255)
                                                                 ttCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAAt                                                                               >  1:780354/1‑150 (MQ=255)
                                                                     agctagctTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCgtgt                                                                                     <  1:118987/140‑1 (MQ=255)
                                                                     agctagctTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCgtgt                                                                                     >  2:118987/1‑140 (MQ=255)
                                                                           ctTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGggtg                                                                     <  1:69137/150‑1 (MQ=255)
                                                                            ttGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGggt                                                                      >  2:504155/1‑148 (MQ=255)
                                                                            ttGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGggt                                                                      <  1:504155/148‑1 (MQ=255)
                                                                                               gAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCAt                                                                                                                             >  2:505754/1‑74 (MQ=255)
                                                                                               gAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCAt                                                                                                                             <  1:505754/74‑1 (MQ=255)
                                                                                                aGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCgg                                                >  1:705542/1‑150 (MQ=255)
                                                                                                  agAGGAATTCGCGCCTCCTAGCGACGCAGCTGCTTAGCTAACGCAAGAATGCGGGGGCCGAGGATTTTCATGGATGTGGCCGTCGAGCTCAGGAGCAGAATAGCTCCGTGTGTGAATGCGCCGGGTGGGACCAGGCTGGTTACGCCGGTa                                              >  2:311679/2‑150 (MQ=255)
                                                                                                     ggCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATcc                                           <  2:705542/150‑1 (MQ=255)
                                                                                                         ttCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAAc                                       <  1:720382/150‑1 (MQ=255)
                                                                                                             aTCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGc                                   >  2:138462/1‑150 (MQ=255)
                                                                                                                   cTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCACAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGTGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTg                             <  1:576036/150‑1 (MQ=255)
                                                                                                                         aCGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGcc                                                  >  2:270798/1‑123 (MQ=255)
                                                                                                                         aCGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGcc                                                  <  1:270798/123‑1 (MQ=255)
                                                                                                                           gCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTGGGCGtttt                     <  2:780354/150‑1 (MQ=255)
                                                                                                                                 gCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTc                                                                                                                >  1:527391/1‑53 (MQ=255)
                                                                                                                                 gCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTc                                                                                                                <  2:527391/53‑1 (MQ=255)
                                                                                                                                    tAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTGGGCGTTTTCATTTTCGa            <  1:705337/150‑1 (MQ=255)
                                                                                                                                             ccACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTGGGCGTTTTCATTTTCGATGCGAGCTa   >  1:742260/1‑150 (MQ=255)
                                                                                                                                              cacaATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGACTGATGCGGGTGGTACCAAGACGCTTAAGGATCTCGCCATACCCTCAGTTTGGGCGTTTTCATTTTCGATGCGAGCTaa  >  1:575296/1‑150 (MQ=255)
                                                                                                                                              cacaATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGtt                                                       >  2:587451/1‑97 (MQ=255)
                                                                                                                                              cacaATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGtt                                                       <  1:587451/97‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACGCATTTGCTTAGCTAAGGCCACAATGGCGAAGCCCAGGATTTTCATGAATGTGGCCGTCCAGCTCAGGAGCAGAATAGCTCCGTGTTTGAATGCTCCGGGTGGTACCAGGCTGGTTACGCCGGTATCCCAACCTGCGGGCTGGGCGTTTTCATTTTCGATGCGAGCTAA  >  NC_002947/743483‑743774

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