Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,192,822 T→A 14.3% intergenic (+155/‑47) PP_4572 → / → PP_4573 hypothetical protein/AAA family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,192,8220TA14.3% 86.2 / 9.2 35intergenic (+155/‑47)PP_4572/PP_4573hypothetical protein/AAA family ATPase
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base A (3/2);  total (18/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 = 8.73e-02

GACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTGGTCAAGGG  >  NC_002947/5192676‑5192966
                                                                                                                                                  |                                                                                                                                                
gACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAGCCTGACAGCTTGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTAc                                                                                                                                               <  2:401573/150‑1 (MQ=255)
            aGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAAACCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                   >  1:230716/1‑150 (MQ=255)
                       aCCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGccc                                                                                                                        <  1:192108/150‑1 (MQ=255)
                         ccGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTAccc                                                                                                                                             <  2:60828/127‑1 (MQ=255)
                         ccGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTAccc                                                                                                                                             >  1:60828/1‑127 (MQ=255)
                               tGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACCAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCTCGGTTAACCCCGCCCCTACAACCTTATTAGGCGCCCGTTTACCCGCGCAGAGGCCCGCAAAGCCCTAAACaga                                                                                                                <  2:35307/150‑1 (MQ=255)
                                 cGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAgg                                                                                                              >  2:664216/1‑150 (MQ=255)
                                   tttGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGAc                                                                                                            >  2:576901/1‑150 (MQ=255)
                                                     cgcgTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCg                                                                                          >  2:305743/1‑150 (MQ=255)
                                                                  gAAACAACTCCAGATTTTTTTCTTCCAGCTGAGGCCGGTTCGCGGGTAAACCCGCCCCTACAACAATCTCCGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGa                                                                             <  1:352563/150‑1 (MQ=255)
                                                                         aTCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGa                                                                                                             >  2:690462/1‑111 (MQ=255)
                                                                         aTCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGa                                                                                                             <  1:690462/111‑1 (MQ=255)
                                                                          tCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCg                                                                     >  1:86293/1‑150 (MQ=255)
                                                                                     tgctgcCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGAc                                                          <  1:305743/150‑1 (MQ=255)
                                                                                            tGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAAc                                                   >  1:128643/1‑150 (MQ=255)
                                                                                                 gggCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCg                                              <  2:128643/150‑1 (MQ=255)
                                                                                                     ctctTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTaa                                          >  2:285550/1‑150 (MQ=255)
                                                                                                              ggTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCg                                              <  1:6079/137‑1 (MQ=255)
                                                                                                              ggTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCg                                              >  2:6079/1‑137 (MQ=255)
                                                                                                                   aCCCGCCCCTACAACCTTATTAGGCCCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTAcc                            <  1:285550/150‑1 (MQ=255)
                                                                                                                     ccGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCt                          >  2:465147/1‑150 (MQ=255)
                                                                                                                       gCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACg                                                                        <  2:27990/102‑1 (MQ=255)
                                                                                                                       gCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACg                                                                        >  1:27990/1‑102 (MQ=255)
                                                                                                                         cccTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGaa                      >  2:43519/1‑150 (MQ=255)
                                                                                                                               aaCCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTCCACTCGAACGCCCg                >  2:388548/1‑150 (MQ=255)
                                                                                                                                  cTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTg             >  1:165979/1‑150 (MQ=255)
                                                                                                                                      ttaGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTgg         <  2:464773/150‑1 (MQ=255)
                                                                                                                                      ctcttcgcgGGTAAACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGcc                  >  2:415523/6‑141 (MQ=255)
                                                                                                                                      ctcttcgcgGGTAAACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGcc                  <  1:415523/136‑1 (MQ=255)
                                                                                                                                        cttcgcgGGTAAACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTGGTc       >  2:201916/4‑150 (MQ=255)
                                                                                                                                         ttcgcgGGTAAACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGaaa                                                    <  2:564013/102‑1 (MQ=255)
                                                                                                                                         ttcgcgGGTAAACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGaaa                                                    >  1:564013/3‑104 (MQ=255)
                                                                                                                                          gcgcGTGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAANCTGGCGNTAAACGCGGCCANTACCCTCGAACGCCCGCTGTTGGTCaa     <  1:388548/150‑1 (MQ=255)
                                                                                                                                          gcgcGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTGGTCaa     <  2:557055/150‑1 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTGGTCAAggg  <  1:465147/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGCCCGCTGTTGGTCAAGGG  >  NC_002947/5192676‑5192966

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