Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,004,649 C→T 11.8% A318V (GCA→GTA)  PP_4410 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,004,6490CT11.8% 83.0 / 3.2 34A318V (GCA→GTA) PP_4410hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (2/2);  total (17/17)
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.32e-01

GGAATCGTTAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCTGTCGAAGTTTAAGTACCCCTCAACACCT  >  NC_002947/5004519‑5004797
                                                                                                                                  |                                                                                                                                                    
ggAATCGTTAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTAttt                                                                                                                                       >  1:336368/1‑146 (MQ=255)
ggAATCGTTAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTAttt                                                                                                                                       <  2:336368/146‑1 (MQ=255)
       ttAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAAc                                                                                                                            >  2:666400/1‑150 (MQ=255)
       ttAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAAc                                                                                                                            >  1:461086/1‑150 (MQ=255)
         aGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCa                                                                                                                          <  1:538358/150‑1 (MQ=255)
          gATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGCTGACTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAg                                                                                                                         <  2:590479/150‑1 (MQ=255)
               cGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGaa                                                                                                                    <  1:666400/150‑1 (MQ=255)
                               gcagcCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCa                                                                                                    >  2:713446/1‑150 (MQ=255)
                                              aaaGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGa                                                                                     >  2:594552/1‑150 (MQ=255)
                                                 gTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGTAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGaa                                                                                                 >  2:221489/1‑135 (MQ=255)
                                                 gTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGTAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGaa                                                                                                 >  2:221721/1‑135 (MQ=255)
                                                 gTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGTAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGaa                                                                                                 <  1:221721/135‑1 (MQ=255)
                                                 gTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGTAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGaa                                                                                                 <  1:221489/135‑1 (MQ=255)
                                                     aGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCaa                                                                              <  2:461086/150‑1 (MQ=255)
                                                                        attcacattcaTCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACttattt                                                                 <  2:243078/144‑1 (MQ=255)
                                                                        attcacattcaTCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACttattt                                                                 >  1:243078/1‑144 (MQ=255)
                                                                                   tcaAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCaa                                                                              >  2:57632/1‑120 (MQ=255)
                                                                                   tcaAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCaa                                                                              <  2:145757/120‑1 (MQ=255)
                                                                                   tcaAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCaa                                                                              <  1:57632/120‑1 (MQ=255)
                                                                                   tcaAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCaa                                                                              >  1:145757/1‑120 (MQ=255)
                                                                                         gCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCtt                                          >  1:592093/1‑150 (MQ=255)
                                                                                         gCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCtt                                          <  2:592093/150‑1 (MQ=255)
                                                                                            ggATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATgg                                                                                              >  1:494078/1‑95 (MQ=255)
                                                                                            ggATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATgg                                                                                              <  2:494078/95‑1 (MQ=255)
                                                                                                     cAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCt                              <  2:751691/150‑1 (MQ=255)
                                                                                                                tGGATGATTGGCTGAATGCAGAGCGTGTCTGGTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCTGTCGAAGTTTa                   >  1:19673/1‑150 (MQ=255)
                                                                                                                     gatTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGa                                              <  1:88761/118‑1 (MQ=255)
                                                                                                                     gatTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGa                                              >  2:88761/1‑118 (MQ=255)
                                                                                                                      atTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCTGTCGAAGTTTAAGTAcc             <  1:713446/150‑1 (MQ=255)
                                                                                                                           cTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCttt                                         <  2:423953/117‑1 (MQ=255)
                                                                                                                           cTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCttt                                         >  1:423953/1‑117 (MQ=255)
                                                                                                                                tGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGagaa                                                                                  <  2:385974/71‑1 (MQ=255)
                                                                                                                                tGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGagaa                                                                                  >  1:385974/1‑68 (MQ=255)
                                                                                                                                 gCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCTGTCGAAGTTTAAGTACCCCTCAACACCt  >  2:605046/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                                    
GGAATCGTTAGATATCGCAAATATGGAATCAGCAGCCTTTCATGCGAAAGTATAGCCGACTATAGGCTTGACATTCACATTCATCAAGGGCTGGATATATTCAAGGCAGCAATGGATGATTGGCTGAATGCAGAGCGTGCTTATTTGAGGCATCAACCAGTTGAAAAATTAAAAGCACGCAGAATGGCAATCGGGAAAATCAATATATACTTATTTATATATCTACCCCTTTGGATCTTTCAAAACCCCCTGTCGAAGTTTAAGTACCCCTCAACACCT  >  NC_002947/5004519‑5004797

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