Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,969 A→C 9.7% T4T (ACT→ACG) ‡ PP_4272 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,857,9690AC 9.7% 122.3 / 2.9 45T4T (ACT→ACG) ‡PP_4272membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base C (3/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA  >  NC_002947/4857820‑4858116
                                                                                                                                                     |                                                                                                                                                    
aaCCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa                                                                                                                                                      <  2:129107/150‑1 (MQ=255)
      cTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCaaa                                                                                                                                                >  1:917378/1‑149 (MQ=255)
           cagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCATTCTGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaa                                                                                                                                          <  1:660051/150‑1 (MQ=255)
            agcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaag                                                                                                                                         <  1:19358/150‑1 (MQ=255)
                aCCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa                                                                                                                                                      >  2:855627/1‑134 (MQ=255)
                aCCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa                                                                                                                                                      <  1:855627/134‑1 (MQ=255)
                    ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGc                                                                                                                                 >  2:946408/1‑150 (MQ=255)
                    ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGc                                                                                                                                 >  1:377580/1‑150 (MQ=255)
                     cacaGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCAAAAAAGTTGAAGTCAACTACAAGTCAAGGGcc                                                                                                                                >  2:308441/1‑150 (MQ=255)
                           tacATCAGTAAAGTGACACCAAAGCAAATCGGCGATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAATCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAAt                                                                                                                          <  2:232218/148‑1 (MQ=255)
                             cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                        <  1:1292023/150‑1 (MQ=255)
                             cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                        <  1:789316/150‑1 (MQ=255)
                              aTCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCAAGTCAAGGGCCCGAAATGGa                                                                                                                       >  2:1157412/1‑150 (MQ=255)
                                cACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACa                                                                                                                     <  1:308441/150‑1 (MQ=255)
                                             acctggCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGc                                                                                                        >  2:895134/6‑150 (MQ=255)
                                              cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                        >  2:1125554/1‑133 (MQ=255)
                                              cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                        <  1:1125554/133‑1 (MQ=255)
                                              aaaTGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCCCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGcc                                                                                                       <  2:568306/149‑1 (MQ=255)
                                                     aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa                                                                                                    <  2:1205114/146‑1 (MQ=255)
                                                     aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa                                                                                                    >  1:1205114/1‑146 (MQ=255)
                                                        cGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTg                                                                                             <  1:946408/150‑1 (MQ=255)
                                                              tGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTAc                                                                                       >  2:403666/1‑150 (MQ=255)
                                                                   cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCg                                                                                                               >  1:1211173/1‑121 (MQ=255)
                                                                   cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCg                                                                                                               <  2:1211173/121‑1 (MQ=255)
                                                                    tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaagtcaa                                                                                                                                     >  2:564219/1‑98 (MQ=255)
                                                                    tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaagtcaa                                                                                                                                     <  1:564219/98‑1 (MQ=255)
                                                                          cATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGacac                                                                           <  2:775838/150‑1 (MQ=255)
                                                                            tCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                         >  2:774272/1‑150 (MQ=255)
                                                                                aGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGtt                                                                     >  2:1244146/1‑150 (MQ=255)
                                                                                              gAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCACAGTAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGaaa                                                       <  1:1157412/150‑1 (MQ=255)
                                                                                              gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat                                                                                                                                          <  1:1053355/68‑4 (MQ=38)
                                                                                              gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat                                                                                                                                          >  2:1053355/1‑65 (MQ=38)
                                                                                                 cAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc                                                    <  1:1244146/150‑1 (MQ=255)
                                                                                                 cAAGCCCGCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc                                                    <  1:567834/150‑1 (MQ=255)
                                                                                                    gCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTg                                                 >  2:361757/1‑150 (MQ=255)
                                                                                                       cTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTc                                              <  2:377580/150‑1 (MQ=255)
                                                                                                        tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCa                                             >  2:236083/1‑150 (MQ=255)
                                                                                                          gAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTAAAGGGCCCGAAATGGACATTGCCGACGTCTCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGACGAAAGGCTTGCTCAAg                                           >  2:1366650/1‑150 (MQ=255)
                                                                                                                      tCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCggg                               >  1:817944/1‑150 (MQ=255)
                                                                                                                                 ccGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaagg                    <  1:774272/150‑4 (MQ=255)
                                                                                                                                   gcgcAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTc                  <  2:817944/150‑1 (MQ=255)
                                                                                                                                           ttgAAGTTGCCTTCAACTTCAAGTCAAGGGCCCGAGATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGttt          >  2:59822/3‑150 (MQ=255)
                                                                                                                                            tgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                         <  2:324015/85‑1 (MQ=255)
                                                                                                                                            tgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                         >  1:324015/2‑86 (MQ=255)
                                                                                                                                            aGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTg         >  1:45114/1‑150 (MQ=255)
                                                                                                                                                gTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc                                  <  2:450076/121‑1 (MQ=255)
                                                                                                                                                gTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc                                  >  1:450076/1‑121 (MQ=255)
                                                                                                                                                   gAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGa  <  1:361757/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA  >  NC_002947/4857820‑4858116

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: