Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,393,681 C→A 8.7% E431* (GAG→TAG)  PP_3869 ← sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,393,6810CA 8.7% 133.4 / 4.0 46E431* (GAG→TAG) PP_3869sheath protein
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TAGCTGGTTCACCAGATCCGGCGGATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGACTGATGCAGCGTCTCG  >  NC_002947/4393535‑4393823
                                                                                                                                                  |                                                                                                                                               
tAGCTGGTTCACCAGATCCGGCGGATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCgg                                                                                                                                               >  7:134962/1‑149 (MQ=255)
      gTTCACCAGATCCGGCGGATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGa                                                                                                                                         >  1:229544/1‑149 (MQ=255)
             aGATCCGGCGGATACAGCACGTTGACCCGGGCGAGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGtt                                                                                                                                  >  5:201459/1‑149 (MQ=255)
              gATCCGGCGGATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGttt                                                                                                                                 <  2:229551/149‑1 (MQ=255)
                  cggcggATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTgggg                                                                                                                             <  2:229544/149‑1 (MQ=255)
                     cggATACAGCACGTTGACCCGGGCGGGGTCTGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGa                                                                                                                          <  5:347076/149‑1 (MQ=255)
                           cAGCACGTTGACCCGGGCGGTGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATc                                                                                                                     <  5:348647/148‑1 (MQ=255)
                            aGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGCTCAGGTACTTGGCGAAGAGATCTGCGTAAACAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGc                                                                                                                   >  1:306490/1‑149 (MQ=255)
                            aGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGc                                                                                                                   >  1:192908/1‑149 (MQ=255)
                            aGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCCGCCCGGAGGACGTTTGGGGTGCCGATCGc                                                                                                                   >  8:341266/1‑149 (MQ=255)
                            aGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGAGATCAGCTCGGCGCGGATGACGTTT‑GGGTGACGATCGcc                                                                                                                  >  6:155902/1‑149 (MQ=255)
                              cACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCt                                                                                                                 >  2:236549/1‑149 (MQ=255)
                                 gTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGAc                                                                                                              >  6:144329/1‑149 (MQ=255)
                                 gTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGAc                                                                                                              <  8:134962/149‑1 (MQ=255)
                                    gACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCa                                                                                                            <  2:131664/148‑1 (MQ=255)
                                       ccGGGCGGGGTCGGTAGTCGAACGCTCTACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAgc‑gc                                                                                                       >  1:209012/1‑149 (MQ=255)
                                       ccGGGCGGGGTCGGTAGTCGAACGCTCAATGATCAGGTACTAGGCGAAGAGATATGCGTTCTCAGCAATGCCAAGCGTCTCAAGCACAAAGTAACCGGCGATCAGCTCGGAGCGGATGACGTTTGTGGTGACGATCGACTGACCAtc‑g                                                                                                        >  4:257154/1‑148 (MQ=255)
                                       ccGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAgc‑gc                                                                                                       <  2:157916/149‑1 (MQ=255)
                                           gcgggGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGaa                                                                                                   <  7:269506/149‑1 (MQ=255)
                                           gcgggGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGaa                                                                                                   <  1:96903/149‑1 (MQ=255)
                                                           aaCGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTAGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGtcgtcg                                                                                   >  5:104024/1‑149 (MQ=255)
                                                           aaCGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTAGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGtcgtcg                                                                                   >  5:98701/1‑149 (MQ=255)
                                                                tCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAg                                                                              >  8:195015/1‑149 (MQ=255)
                                                                  aaCGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCt                                                                            >  7:62058/1‑149 (MQ=255)
                                                                     gATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTtgt                                                                         <  5:53242/149‑1 (MQ=255)
                                                                      agcAGGGACTTGGCGAAGAGATCTGCCGTCTCAGCAATGACCAGCGTCTCTAGCACGTTGTAACCGGCGATCAGCGCGGCGCGGTTGACGTTTGGTGTGACGATCGCCTGACCAGC‑GCCGAATCGCGTGACGTCGTCGGCCAGCTtgtg                                                                        <  8:325005/147‑1 (MQ=255)
                                                                       tCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGCTGCCGAAGCGCGTGCCGTCGTCGGCCAGCTtgtg                                                                        <  1:236549/149‑1 (MQ=255)
                                                                         aGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTgccg                                                                     >  5:252395/1‑149 (MQ=255)
                                                                          ggTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTgccgc                                                                    <  2:192908/149‑1 (MQ=255)
                                                                                        agaTCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGt                                                      <  6:252395/149‑1 (MQ=255)
                                                                                        agaTCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGt                                                      <  8:62058/149‑1 (MQ=255)
                                                                                                   ctcAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACtt                                                            <  2:155046/132‑1 (MQ=255)
                                                                                                   ctcAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACtt                                                            >  1:155046/1‑132 (MQ=255)
                                                                                                             cccAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAAtt                                 <  2:209012/149‑1 (MQ=255)
                                                                                                                  cGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATtgatga                            <  7:195015/149‑1 (MQ=255)
                                                                                                                     ctcAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGatgatg                          >  4:80711/1‑148 (MQ=255)
                                                                                                                      tcAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGatgatg                          >  5:115003/1‑147 (MQ=255)
                                                                                                                        aaGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAgg                      >  7:279841/1‑149 (MQ=255)
                                                                                                                            aCGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCTTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGa                  <  6:201459/149‑1 (MQ=255)
                                                                                                                                    acgggCGATCTGGTGGGCGCGGATGAGGTTTGGGGTGACGATCTCCTGACCAGC‑GCCGAAGCGCGTGCCGTCTTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATTTAGGCCGACTGATGc           <  2:306490/145‑1 (MQ=255)
                                                                                                                                    aCCGGCGATCAGCTAGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATTTAGGCCGACTGATGCa          <  6:98701/149‑1 (MQ=255)
                                                                                                                                    aCCGGCGATCAGCTAGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGACTGATGCa          <  6:104024/149‑1 (MQ=255)
                                                                                                                                      cGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGaaa                                   <  6:295849/122‑1 (MQ=255)
                                                                                                                                      cGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGaaa                                   >  5:295849/1‑122 (MQ=255)
                                                                                                                                           aTCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGACTGATGCAGCGtctc   >  4:363982/1‑149 (MQ=255)
                                                                                                                                            tCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGACTGATGCAGCGTCTCg  >  8:279604/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TAGCTGGTTCACCAGATCCGGCGGATACAGCACGTTGACCCGGGCGGGGTCGGTAGTCGAACGCTCAACGATCAGGTACTTGGCGAAGAGATCTGCGTTCTCAGCAATGCCCAGCGTCTCAAGCACGTAGTAACCGGCGATCAGCTCGGCGCGGATGACGTTTGGGGTGACGATCGCCTGACCAGC‑GCCGAAGCGCGTGCCGTCGTCGGCCAGCTTGTGCCGCCCGTACTTGCTGGTGATTCGCCCCTTAAGGAAATTGATGATGTAGGCCGACTGATGCAGCGTCTCG  >  NC_002947/4393535‑4393823

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