Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,597,312 A→T 8.7% I267N (ATT→AAT)  PP_3174 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,597,3120AT 8.7% 133.3 / 3.4 46I267N (ATT→AAT) PP_3174oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base T (2/2);  total (24/22)
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.84e-01

ATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCTTCCGGCACCGCGCCCGG  >  NC_002947/3597164‑3597460
                                                                                                                                                    |                                                                                                                                                     
aTGCCCTGGCGTTACCAACGAATAGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCaa                                                                                                                                                       >  4:219366/1‑149 (MQ=255)
  gCCCTGGCGTTACCAACGAATCGATGAACGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCATTg                                                                                                                                                     >  3:74373/1‑149 (MQ=255)
   cccTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATgg                                                                                                                                                    >  8:206014/1‑149 (MQ=255)
   cccTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATg                                                                                                                                                     <  2:112310/148‑1 (MQ=255)
             aCCAACGAATAGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGa                                                                                                                                           <  3:219366/148‑1 (MQ=255)
              ccAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGaa                                                                                                                                                >  1:145455/1‑142 (MQ=255)
              ccAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGaa                                                                                                                                                <  2:145455/142‑1 (MQ=255)
              ccAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGAc                                                                                                                                          >  3:197173/1‑148 (MQ=255)
                 aCGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATg                                                                                                                                      <  6:91976/149‑1 (MQ=255)
                  cGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGc                                                                                                                                     <  7:206014/149‑1 (MQ=255)
                     aTCTATGAACGCATGATCCTCAGCGGTATTCTTCGCAGCCAATGATGCGCCATTGTTATCCCACTGCTCCTCGGTCCGCGGCCGGACTATCGCGGAAGTGACCCATTGGTTATTCAACACCCAGGCTTTTTCGAATTAGCCATTGcctg                                                                                                                                  <  4:74373/149‑3 (MQ=255)
                         aTGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCAcccc                                                                                                                              >  6:212911/1‑149 (MQ=255)
                                 cTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCct                                                                                                                      >  7:53186/1‑149 (MQ=255)
                                  tCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCctc                                                                                                                     >  7:53818/1‑149 (MQ=255)
                                       ctccGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCAt                                                                                                                <  3:349185/149‑1 (MQ=255)
                                           gcgGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGt                                                                                                            >  5:288316/1‑149 (MQ=255)
                                               tAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGa                                                                                                        <  2:33560/149‑1 (MQ=255)
                                                    ttCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTg                                                                                                  >  4:307453/1‑149 (MQ=255)
                                                       gccgccAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGcc                                                                                                <  5:212911/148‑1 (MQ=255)
                                                                 cgccgcCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGt                                                                                      <  3:307453/148‑1 (MQ=255)
                                                                   ccgccATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAg                                                                                   <  5:343787/149‑1 (MQ=255)
                                                                       cATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCATTGTCGAATTCGCCATTGcctgg                                                                                                                                 >  2:139195/1‑97 (MQ=38)
                                                                       cATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCATTGTCGAATTCGCCATTGcctgg                                                                                                                                 <  1:139195/100‑4 (MQ=38)
                                                                                   aCTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGccac                                                                   <  4:52422/149‑1 (MQ=255)
                                                                                        tcctcGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGccactccac                                                              >  8:235819/1‑149 (MQ=255)
                                                                                          ctcGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAAGGATTCCGGGCAATTGTAAGCGACGACATGCGccactacaat                                                             >  1:126364/1‑147 (MQ=255)
                                                                                          ctcGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCt                                                                           <  7:62796/134‑1 (MQ=255)
                                                                                          ctcGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCt                                                                           >  8:62796/1‑134 (MQ=255)
                                                                                          ctcGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACtt                                                            <  6:288316/149‑1 (MQ=255)
                                                                                            cGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCATGCAATGGCGAATTCGCCAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTcc                                                          <  8:53818/149‑1 (MQ=255)
                                                                                                      cccGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGtt                                                >  6:318055/1‑149 (MQ=255)
                                                                                                       ccGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGttt                                               >  7:156411/1‑149 (MQ=255)
                                                                                                                           ccATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGgc                                                                        <  6:228190/104‑1 (MQ=255)
                                                                                                                           ccATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGgc                                                                        >  5:228190/1‑104 (MQ=255)
                                                                                                                           ccATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGccc                                   <  1:116686/141‑1 (MQ=255)
                                                                                                                           ccATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGccc                                   >  2:116686/1‑141 (MQ=255)
                                                                                                                           ccATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGc                           >  2:302538/1‑149 (MQ=255)
                                                                                                                              ttGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTAc                        >  5:342933/1‑149 (MQ=255)
                                                                                                                                 gTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCgg                     >  6:326475/1‑149 (MQ=255)
                                                                                                                                 gTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCg                      <  3:338398/148‑1 (MQ=255)
                                                                                                                                    aTTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCtt                  >  4:252331/1‑149 (MQ=255)
                                                                                                                                       cAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCTTcc                <  5:318055/148‑1 (MQ=255)
                                                                                                                                         acacCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTg                                                                                                                   <  5:274651/48‑1 (MQ=255)
                                                                                                                                         acacCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTg                                                                                                                   >  6:274651/1‑48 (MQ=255)
                                                                                                                                             ccAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGttt                                               >  8:164900/1‑111 (MQ=255)
                                                                                                                                             ccAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGttt                                               <  7:164900/111‑1 (MQ=255)
                                                                                                                                                  cAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCTTCCGGCACCGCGccc    >  6:367592/1‑149 (MQ=255)
                                                                                                                                                    aTGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCTTCCGGCACCGCGCCCgg  <  4:197173/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATGCGCCGCCATAGGTATCCCACTGCTCCTCGGTCCGCGGCCCGACAATCGCGGAACTGACCCATTGGTTATTCAACACCCAGGCAATGGCGAATTCGACAATGCCCACCCCCTTGGCCTCTGCATACGTCTGAA‑TCTGCCGGGCAATGGTCAGCGACTCCTGGCGCCACTCCACTTCCATGATGCGTTTGTCCTGCCGCCCTGCTCGGCTACCGGCTTCCGGCACCGCGCCCGG  >  NC_002947/3597164‑3597460

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