Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,402,344:1 +A 9.8% coding (66/330 nt) PP_2103 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,402,3441.A 9.8% 146.5 / 11.6 41coding (66/330 nt)PP_2103hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CCACCACGGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCC  >  NC_002947/2402199‑2402487
                                                                                                                                                  |                                                                                                                                               
ccaccaCGGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATgcgc‑aggc                                                                                                                                             <  1:811072/150‑1 (MQ=255)
 caccacGGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATgcgc‑aggcg                                                                                                                                            >  2:262101/1‑150 (MQ=255)
   ccacGGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATgcgc‑aggcgca                                                                                                                                          <  1:1161246/150‑1 (MQ=255)
             gATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                <  1:668035/150‑1 (MQ=255)
                  cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCCTGCGC‑AGGCGCAGCAGCGCCTTGCGCa                                                                                                                           <  2:979764/150‑1 (MQ=255)
                           gcgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgcg                                                                                                                  >  1:1098240/1‑150 (MQ=255)
                                     gATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                   >  2:67341/1‑139 (MQ=255)
                                     gATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                   <  1:67341/139‑1 (MQ=255)
                                         cccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                    >  2:138648/1‑150 (MQ=255)
                                         cccAGGCCTTAATGTAGCGCGCTGCCGATACTGCGCATTCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                    <  2:316341/150‑1 (MQ=255)
                                          ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                   >  2:242759/1‑150 (MQ=255)
                                           cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                             >  2:811065/1‑123 (MQ=255)
                                           cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                             <  1:811065/123‑1 (MQ=255)
                                           cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                  >  2:1045901/1‑150 (MQ=255)
                                           cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCATGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                  <  1:990039/150‑1 (MQ=255)
                                            aGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGc                                                                                                 >  1:652830/1‑150 (MQ=255)
                                                       tAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAg                                                                                      <  2:59649/150‑1 (MQ=255)
                                                       tAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAg                                                                                      <  1:605123/150‑1 (MQ=255)
                                                            aGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatg                                                                                                                         <  1:843425/110‑3 (MQ=38)
                                                            aGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatg                                                                                                                         >  2:843425/1‑108 (MQ=38)
                                                               tGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                              <  1:138648/150‑1 (MQ=255)
                                                                   cATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAAc                                                                          >  1:445557/1‑150 (MQ=255)
                                                                             aCGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTccc                                                                >  1:472287/1‑150 (MQ=255)
                                                                               gcggcgCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTg                                                              >  2:653977/1‑150 (MQ=255)
                                                                                       gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                      <  1:242759/150‑1 (MQ=255)
                                                                                         ggTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCaa                                                    >  2:811324/1‑150 (MQ=255)
                                                                                              aaGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcg                                               >  2:416429/1‑150 (MQ=255)
                                                                                               aGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgt                                              >  2:962316/1‑150 (MQ=255)
                                                                                                gCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgtc                                             <  1:1045901/150‑1 (MQ=255)
                                                                                                 cAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcg                                                                                                                    >  1:885986/1‑78 (MQ=255)
                                                                                                 cAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcg                                                                                                                    <  2:885986/78‑1 (MQ=255)
                                                                                                   aCTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                          <  2:759383/150‑1 (MQ=255)
                                                                                                    cTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                         <  1:653977/150‑1 (MQ=255)
                                                                                                       aCCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTg                                      >  1:417961/1‑150 (MQ=255)
                                                                                                                      gATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                       <  2:1098240/150‑1 (MQ=255)
                                                                                                                                    gCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAacac         <  1:416429/150‑1 (MQ=255)
                                                                                                                                      aTTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccg       >  1:411444/1‑150 (MQ=255)
                                                                                                                                        ttCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgcc     <  1:319701/150‑1 (MQ=255)
                                                                                                                                         gaactGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgccg    >  1:56145/5‑150 (MQ=255)
                                                                                                                                           cATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCACATCTGGCAGGCTCATGGTTACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCgccgcc  >  1:147749/1‑150 (MQ=255)
                                                                                                                                            aTGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgatgat                    >  2:461723/1‑131 (MQ=255)
                                                                                                                                            aTGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgatgat                    <  1:461723/131‑1 (MQ=255)
                                                                                                                                             tGCGCAAGGCGCTGCTGCGCCTT‑CGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                       <  1:384963/127‑1 (MQ=255)
                                                                                                                                             tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                       >  2:384963/1‑127 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCACCACGGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCC  >  NC_002947/2402199‑2402487

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