Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,844,171 C→T 11.1% A246V (GCC→GTC)  PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,844,1710CT11.1% 113.6 / 5.6 45A246V (GCC→GTC) PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCT  >  NC_002947/5844026‑5844311
                                                                                                                                                 |                                                                                                                                            
ccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGGCGTCCAg                                                                                                                                          <  2:196181/150‑1 (MQ=255)
ccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGGCGTCCAg                                                                                                                                          <  2:196285/150‑1 (MQ=255)
ccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                          <  2:398927/150‑1 (MQ=255)
 cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCTCAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAgg                                                                                                                                         <  1:506513/150‑1 (MQ=255)
 cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAgg                                                                                                                                         <  1:55658/150‑1 (MQ=255)
  ggTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAggg                                                                                                                                        >  2:459166/1‑150 (MQ=255)
  ggTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAggg                                                                                                                                        <  1:359013/150‑1 (MQ=255)
  ggTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAggg                                                                                                                                        >  2:367842/1‑150 (MQ=255)
            aaaCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCt                                                                                                                              >  2:167454/1‑150 (MQ=255)
               cAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTAcc                                                                                                                           >  2:208820/1‑150 (MQ=255)
                  tCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggc                                                                                                                        >  2:339263/1‑150 (MQ=255)
                        cTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcagg                                                                                                                     >  1:370345/1‑147 (MQ=255)
                        cTGGCGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcagg                                                                                                                     <  2:370345/147‑1 (MQ=255)
                              aCCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGc                                                                                                            <  1:600205/150‑1 (MQ=255)
                                ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCATGCCCTGGGCGTCCAGGGcctgg                                                                                                                                   >  2:240722/1‑121 (MQ=255)
                                ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCATGCCCTGGGCGTCCAGGGcctgg                                                                                                                                   <  1:240722/125‑5 (MQ=255)
                                ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                                          >  2:250892/1‑150 (MQ=255)
                                               ccTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcag                                                                                                                      >  2:615117/1‑123 (MQ=255)
                                               ccTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcag                                                                                                                      <  1:615117/123‑1 (MQ=255)
                                                cTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGGCGGTTTTGTCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAc                                                                                          >  1:195253/1‑150 (MQ=255)
                                                 tCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTAGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCTAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcc                                                                                         <  2:565243/150‑1 (MQ=255)
                                                    aGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctcc                                                                                      >  2:389158/1‑150 (MQ=255)
                                                     gCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctcct                                                                                     <  1:389158/150‑1 (MQ=255)
                                                               gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTCGc                                                                           <  1:616406/150‑1 (MQ=255)
                                                               gAACTGCTGTATCTGGCGCGGGCAGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCGCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTAGCCGGTGGc                                                                           >  2:114107/1‑150 (MQ=255)
                                                                     ctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATc                                                                     >  2:78026/1‑150 (MQ=255)
                                                                      tgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCa                                                                    >  1:31559/1‑150 (MQ=255)
                                                                       gTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcc                                                                                         <  1:321625/128‑1 (MQ=255)
                                                                       gTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcc                                                                                         >  2:321625/1‑128 (MQ=255)
                                                                        tATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAAc                                                                  <  1:339263/150‑1 (MQ=255)
                                                                        tATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAAc                                                                  <  1:167454/150‑1 (MQ=255)
                                                                         aTCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACa                                                                 <  1:459166/150‑1 (MQ=255)
                                                                          tCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAg                                                                >  1:326424/1‑150 (MQ=255)
                                                                                gcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGc                                                          >  1:316198/1‑150 (MQ=255)
                                                                                       gAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctgctg                                                   >  2:27873/1‑150 (MQ=255)
                                                                                             ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTgcgc                                             <  1:250892/150‑1 (MQ=255)
                                                                                             ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTgcgc                                             >  2:575090/1‑150 (MQ=255)
                                                                                             ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATGACCTCCGGGCCGGTGTCGTTATCAGCAGCCCTGCCCTGCTCCTgcgc                                             >  2:436647/1‑150 (MQ=255)
                                                                                                                   ttCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATg                       <  1:208820/150‑1 (MQ=255)
                                                                                                                             gACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTggg             <  1:367842/150‑1 (MQ=255)
                                                                                                                                   gtatgCCCTGGGCGTCCAGGGCCTACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCg       >  1:310828/5‑150 (MQ=255)
                                                                                                                                   ccAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                       <  1:195737/70‑1 (MQ=255)
                                                                                                                                   ccAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                       >  2:195737/1‑70 (MQ=255)
                                                                                                                                     aGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGaa     >  1:53215/1‑150 (MQ=255)
                                                                                                                                        cccTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCt  <  2:53215/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCT  >  NC_002947/5844026‑5844311

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: ‑