Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,963,907 C→A 10.7% I217I (ATC→ATA PP_5229 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,963,9070CA10.7% 131.5 / 6.2 47I217I (ATC→ATAPP_5229hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (3/2);  total (24/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 = 6.03e-01

GTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGCGACCT  >  NC_002947/5963768‑5964055
                                                                                                                                           |                                                                                                                                                    
gTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGc                                                                                                                                            >  1:356859/1‑150 (MQ=255)
  gaagaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCtt                                                                                                                                          >  1:665109/1‑150 (MQ=255)
        cAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCAGTACCCTTTTCCTCGGCTATACCGCCGAGGTGCTTGGCcgc                                                                                                                                    >  2:573254/1‑150 (MQ=255)
          gcgcCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGt                                                                                                                                  >  2:130472/1‑150 (MQ=255)
           cgcCGGGAAGTAGGCACCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTCTGTCTCGGCTATATCGCCGAGGTGCTTGGCCGCGtt                                                                                                                                 <  1:698876/150‑1 (MQ=255)
             ccGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACGCTACAAGAGCAGCGTCGGTGCCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                               <  1:446942/150‑1 (MQ=255)
                    gTAGGCTCCGGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAcgc                                                                                                                        <  2:1101311/150‑1 (MQ=255)
                     aaGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                               >  2:557174/2‑142 (MQ=255)
                     aaGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                               <  1:557174/141‑1 (MQ=255)
                         cTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAc                                                                                                                   <  2:705213/150‑1 (MQ=255)
                                 ccGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGAccc                                                                                                           <  2:665109/150‑1 (MQ=255)
                                      gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTAGAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCtt                                                                                                                                          >  2:867277/1‑114 (MQ=255)
                                      gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCtt                                                                                                                                          <  1:867277/114‑1 (MQ=255)
                                               gAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGc                                                                                             <  2:787055/150‑1 (MQ=255)
                                                    ccAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtgatg                                                                                        >  2:765105/1‑150 (MQ=255)
                                                       ggggCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGaaa                                                                                     <  2:356859/150‑1 (MQ=255)
                                                                             aTCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTgcg                                                               <  2:319770/150‑1 (MQ=255)
                                                                                     ccGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTg                                                       >  2:368724/1‑150 (MQ=255)
                                                                                       gcgACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCGCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTgcg                                                     >  2:807737/1‑150 (MQ=255)
                                                                                              gCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAg                                              >  2:891238/1‑150 (MQ=255)
                                                                                               cAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                             <  1:368724/150‑1 (MQ=255)
                                                                                                   acTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCa                                         <  1:765105/150‑1 (MQ=255)
                                                                                                    cTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtga                                                                                          >  1:305932/1‑100 (MQ=255)
                                                                                                    cTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtga                                                                                          <  2:305932/100‑1 (MQ=255)
                                                                                                     tACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                               <  1:292400/62‑1 (MQ=255)
                                                                                                     tACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                               >  2:292400/1‑62 (MQ=255)
                                                                                                     tACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAtt                                                                  >  1:124979/1‑123 (MQ=255)
                                                                                                     tACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAtt                                                                  <  2:124979/123‑1 (MQ=255)
                                                                                                      aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                                      <  1:130472/150‑1 (MQ=255)
                                                                                                      aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                                      <  1:891238/150‑1 (MQ=255)
                                                                                                        aaGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGTGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTg                                                                 <  2:213234/121‑1 (MQ=255)
                                                                                                        aaGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGAGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTg                                                                 >  1:213234/1‑121 (MQ=255)
                                                                                                                     gTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTc                       >  2:1122710/1‑150 (MQ=255)
                                                                                                                     gTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTc                       <  1:1122710/150‑1 (MQ=255)
                                                                                                                               ccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTa                                                                    <  2:984321/95‑1 (MQ=255)
                                                                                                                               ccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTa                                                                    >  1:984321/1‑95 (MQ=255)
                                                                                                                               ccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGa                                >  2:979183/1‑131 (MQ=255)
                                                                                                                               ccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTAc             >  2:633156/1‑150 (MQ=255)
                                                                                                                               ccTCGGCGATATAGCCGAGGTGCTTGGCCGAGTTTTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGa                                <  1:979183/131‑1 (MQ=255)
                                                                                                                                cTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTAcc            >  1:56466/1‑150 (MQ=255)
                                                                                                                                 tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCg           >  1:492376/1‑150 (MQ=255)
                                                                                                                                    gCTATATCGCCGAGGTGCTGGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTcg        >  2:170425/1‑150 (MQ=255)
                                                                                                                                      tataTCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGa                                                                          >  1:857845/1‑82 (MQ=255)
                                                                                                                                      tataTCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGa                                                                          <  2:857845/82‑1 (MQ=255)
                                                                                                                                      tataTCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTcgcg      >  1:679157/1‑150 (MQ=255)
                                                                                                                                         atCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGCGAcc   >  2:1085631/1‑150 (MQ=255)
                                                                                                                                          tCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGCGACCt  <  1:170425/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGCGACCT  >  NC_002947/5963768‑5964055

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