Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 399,531 A→G 12.9% V19A (GTG→GCG) ‡ PP_0332 ← acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947399,5310AG12.9% 72.7 / 2.5 31V19A (GTG→GCG) ‡PP_0332acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base G (2/2);  total (17/14)
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.90e-01

GGGTAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACAAA  >  NC_002947/399382‑399679
                                                                                                                                                     |                                                                                                                                                    
gggTAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTCGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCa                                                                                                                                                      <  2:257727/150‑1 (MQ=255)
   gtgggATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACgg                                                                                                                                                   <  2:318413/148‑1 (MQ=255)
           cacaGGGCGAAGATGGCTTTGTCGATCTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatc                                                                                                                                           <  1:279669/150‑3 (MQ=255)
            acaGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCa                                                                                                                                                      >  2:48150/1‑138 (MQ=255)
            acaGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCa                                                                                                                                                      <  1:48150/138‑1 (MQ=255)
                ggCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatct                                                                                                                                          >  1:147316/1‑143 (MQ=255)
                ggCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatct                                                                                                                                          <  2:147316/146‑4 (MQ=255)
                              tgttcagttCCGACCCGGTTTCGCCCACCCACATCACCGCCCAGGGTGCCTAGTCGCGGCCCCGGCCTGGCAACACCAGCTTGAGCTGTGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTga                                                                                                                        <  1:25860/143‑1 (MQ=255)
                                             aGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAggg                                                                                                                                 >  2:196447/1‑126 (MQ=255)
                                             aGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAggg                                                                                                                                 <  1:196447/126‑1 (MQ=255)
                                                tttCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGc                                                                                                      >  2:350278/1‑150 (MQ=255)
                                                          ccacATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATaaa                                                                                            >  2:370783/1‑150 (MQ=255)
                                                               tCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAag                                                                                           <  2:370473/146‑1 (MQ=255)
                                                               tCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAag                                                                                           >  1:370473/1‑146 (MQ=255)
                                                                       cACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttgttg                                                                               >  1:356214/1‑150 (MQ=255)
                                                                        acgacgCGTAGTCGCGGCCCCGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACGACGGTATTGAGCAGCAGCAGGGTGGTAGTTAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTtgttgt                                                                              <  2:214463/150‑1 (MQ=255)
                                                                               gTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                       <  2:356214/150‑1 (MQ=255)
                                                                                       ggCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGtt                                                               >  1:109781/1‑150 (MQ=255)
                                                                                                  ggCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCg                                                                                                                <  2:261979/90‑1 (MQ=255)
                                                                                                  ggCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCg                                                                                                                >  1:261979/1‑90 (MQ=255)
                                                                                                     aaCACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCa                                                 >  1:155635/1‑150 (MQ=255)
                                                                                                      acacCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAg                                                <  2:155635/150‑1 (MQ=255)
                                                                                                          cAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCg                                            >  1:103429/1‑150 (MQ=255)
                                                                                                               tGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGc                                       <  1:350278/150‑1 (MQ=255)
                                                                                                                    aGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACg                                  >  1:95372/1‑150 (MQ=255)
                                                                                                                           aGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGc                                                  <  2:148386/127‑1 (MQ=255)
                                                                                                                           aGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGc                                                  >  1:148386/1‑127 (MQ=255)
                                                                                                                            gACCAGCAGCGGGTAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGt                          >  1:381545/1‑150 (MQ=255)
                                                                                                                               cagcagCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTtgtg                                                                             <  2:144499/96‑1 (MQ=255)
                                                                                                                               cagcagCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTtgtg                                                                             >  1:144499/1‑96 (MQ=255)
                                                                                                                                            ctCAATACCGTGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAcggcg          >  1:331869/2‑150 (MQ=255)
                                                                                                                                            aTCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTtgttgt                                                                              >  1:97702/1‑82 (MQ=255)
                                                                                                                                            aTCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTtgttgt                                                                              <  2:97702/82‑1 (MQ=255)
                                                                                                                                                    cACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACaaa  >  1:60061/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACAAA  >  NC_002947/399382‑399679

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