Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,448 A→C 8.6% intergenic (‑119/+109) PP_5038 ← / ← PP_5039 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,4480AC 8.6% 124.2 / 2.0 44intergenic (‑119/+109)PP_5038/PP_5039lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base C (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGCGAGGGAACGCAACGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCC  >  NC_002947/5741307‑5741590
                                                                                                                                             |                                                                                                                                              
aGGGCGAGGGAACGCAACGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACATGAAAATCCCTCCCAGGTTGCCTGCCCCACCCCCTATTTCCCAACCCCATTGGTCGGATGGCCCCTACtgt                                                                                                                                         >  4:119706/1‑149 (MQ=255)
  ggCGAGGGAACGCAACGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCactacaga                                                                                                                                         >  5:2840/1‑144 (MQ=255)
  ggCGAGGGAACGCAACGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCacta                                                                                                                                             <  6:2840/143‑1 (MQ=255)
               accGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAAt                                                                                                                                <  8:92079/141‑1 (MQ=255)
               aaCGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACAGTGCGATTAATcaga                                                                                                                            >  7:92079/1‑145 (MQ=255)
               aaCGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTAAACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACaaa                                                                                                                          >  3:88455/1‑149 (MQ=255)
                       gAAGGGACACTCCTGAAAGATGTTACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTAAGTCCCTACACCAATGGTCGGAAGGCCACGACAGTGCGATTAAACACAAAGCaataat                                                                                                                  >  7:199440/1‑149 (MQ=255)
                        cagcgACACTCCTGTAAGTTGTTACTGGCGTGGGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATc                                                                                                                 <  8:199440/145‑1 (MQ=255)
                            gACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATgc                                                                                                             <  2:227224/149‑1 (MQ=255)
                            gACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCAATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATgc                                                                                                             <  4:76209/149‑1 (MQ=255)
                                 tCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCAtt                                                                                                                     >  3:121211/1‑136 (MQ=255)
                                 tCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCAtt                                                                                                                     <  4:121211/136‑1 (MQ=255)
                                             atACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGgcc                                                                                            <  2:100752/149‑1 (MQ=255)
                                                              gCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGAt                                                                                                                                    <  8:57040/92‑1 (MQ=255)
                                                              gCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGAt                                                                                                                                    >  7:57040/1‑92 (MQ=255)
                                                                  tGACCTTCGACATGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCtt                                                                       <  3:119706/149‑1 (MQ=255)
                                                                         cGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCCCCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAATGCCGCCCCCAATGGGAGCGGCCTTGCGTcgc                                                                >  2:287288/1‑149 (MQ=255)
                                                                             gTGAAAATCCCTCCCAGGTTCCCTGCCCCCCCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACGGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCCCCCCCAATGGGAGCGGCCTTGCGGCGCGCTg                                                            >  4:209405/1‑149 (MQ=255)
                                                                                      cctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcag                                                   >  8:242739/1‑149 (MQ=255)
                                                                                       ctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                   >  7:78941/1‑116 (MQ=255)
                                                                                       ctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                   <  8:78941/116‑1 (MQ=255)
                                                                                          ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCAt                                                                                                               <  5:126678/85‑1 (MQ=255)
                                                                                          ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCAt                                                                                                               >  6:126678/1‑85 (MQ=255)
                                                                                           cAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCATTAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                              >  3:38940/1‑149 (MQ=255)
                                                                                                 cccTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCCCTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTg                                        >  6:190589/1‑149 (MQ=255)
                                                                                                 cccTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTg                                        >  6:186207/1‑149 (MQ=255)
                                                                                                  ccTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTCCGTCGCGATGCCTGCCCCGCAGGCCCCTTgg                                       >  2:225751/1‑149 (MQ=255)
                                                                                                       cccACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTc                                  <  7:242739/149‑1 (MQ=255)
                                                                                                         cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCtt                                                                       <  5:192077/110‑1 (MQ=255)
                                                                                                         cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCtt                                                                       >  6:192077/1‑110 (MQ=255)
                                                                                                         cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCt                                   >  2:28621/1‑146 (MQ=255)
                                                                                                              ccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTgg                                       >  8:105985/1‑137 (MQ=255)
                                                                                                              ccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTgg                                       <  7:105985/137‑1 (MQ=255)
                                                                                                               cTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGTAGGCCCCTTGGGCCTCAGAACCgg                          <  2:241471/149‑1 (MQ=255)
                                                                                                               cTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCgg                          <  8:29035/149‑1 (MQ=255)
                                                                                                                  tGTCCCTACACCATTGGTCGGATGGCCCCTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                        <  5:190589/148‑1 (MQ=255)
                                                                                                                  tGTCCCTACACCATTGGTCGGATGGCCCCTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                        <  5:186207/148‑1 (MQ=255)
                                                                                                                       cTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAgg                  >  3:276625/1‑149 (MQ=255)
                                                                                                                          cacCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCt                <  3:144677/148‑1 (MQ=255)
                                                                                                                          cacCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTc               <  1:287288/149‑1 (MQ=255)
                                                                                                                             cATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAg            >  7:133562/1‑149 (MQ=255)
                                                                                                                                   tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTc               >  6:291144/1‑140 (MQ=255)
                                                                                                                                   tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTc               <  5:291144/140‑1 (MQ=255)
                                                                                                                                      gATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAGAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc   <  4:88455/149‑1 (MQ=255)
                                                                                                                                       aTGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTcc  <  4:276625/149‑1 (MQ=255)
                                                                                                                                       aTGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc   >  2:301256/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                              
AGGGCGAGGGAACGCAACGACATGAAGGGACACTCCTGAAAGATGATACAGGCGAGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCC  >  NC_002947/5741307‑5741590

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