Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,508 A→T 7.6% intergenic (‑179/+49) 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,5080AT 7.6% 145.7 / 3.1 53intergenic (‑179/+49)PP_5038/PP_5039lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/22);  new base T (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGG  >  NC_002947/5741373‑5741654
                                                                                                                                       |                                                                                                                                                  
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                                                                                                                                                   <  8:78941/116‑1 (MQ=255)
                     ctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                                                                                   >  7:78941/1‑116 (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                                                                                                                                       >  6:192077/1‑110 (MQ=255)
                                       cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCtt                                                                                                                                       <  5:192077/110‑1 (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)
                                                                                  tgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCg                                                                                           >  1:197882/1‑111 (MQ=255)
                                                                                  tgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCg                                                                                           <  2:197882/111‑1 (MQ=255)
                                                                                   gCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                        >  7:206626/1‑129 (MQ=255)
                                                                                   gCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                        <  8:206626/129‑1 (MQ=255)
                                                                                              caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGc                                                        >  2:88791/1‑134 (MQ=255)
                                                                                              caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGc                                                        <  1:88791/134‑1 (MQ=255)
                                                                                              caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCAATCTCTTCGCTGGTGCTCGGCCGCt                                         >  5:263273/1‑149 (MQ=255)
                                                                                                        aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTgg                                                                                                       <  5:41807/77‑1 (MQ=255)
                                                                                                        aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTgg                                                                                                       >  6:41807/1‑77 (MQ=255)
                                                                                                        aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAg                               >  2:308861/1‑149 (MQ=255)
                                                                                                                    aTAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGCTGCCTGCACAGCAGGCCCCTTGGGCCTCAGCACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcgg                   >  1:154846/1‑149 (MQ=255)
                                                                                                                        aaaaGGCCGCCCGCTATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCTGGCCCGTTGGGCCTCTGAACCGGGACCCAGGCTCCATGAAAAAGTCCATCTGTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCTGGaa                <  5:188262/148‑1 (MQ=255)
                                                                                                                        aaaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAACTCTTCGCTGGTGCTCTGCCGCTCCACCACCACGTTGCGGTGCGGGaa                >  6:188262/1‑148 (MQ=255)
                                                                                                                          aaGGCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCAGGACCCAGGCTCCAGCAAAAAGt                                                                    <  1:86501/94‑1 (MQ=38)
                                                                                                                          aaGGCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCAGGACCCAGGCTCCAGCAAAAAGt                                                                    >  2:86501/1‑94 (MQ=38)
                                                                                                                              ccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                         <  4:303552/133‑1 (MQ=255)
                                                                                                                              ccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                         >  3:303552/1‑133 (MQ=255)
                                                                                                                                 cccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAg                                                                     >  3:87807/1‑86 (MQ=255)
                                                                                                                                 cccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAg                                                                     <  4:87807/86‑1 (MQ=255)
                                                                                                                                 cccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAg      >  1:268477/1‑149 (MQ=255)
                                                                                                                                    ccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCATCAAAAAGTCCATCTCTTCGCTGGTGCTCGGACGCTACAGAACCAGGTTGCGGTGCGTGAAGCGCCCGAAGCgg   >  6:181084/1‑149 (MQ=255)
                                                                                                                                     cAATGGGAGGGGCCTTGCGTCGCGATGCGTGCTCAGCAGGCCCCTTGGGCCTCAGAACCTGGACCCAGGCTCCAGTAAAAAGTCCATCTCTTCTCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCggg  <  6:263273/149‑1 (MQ=255)
                                                                                                                                     cAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTcc                                       <  7:188018/112‑1 (MQ=255)
                                                                                                                                     cAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTcc                                       >  8:188018/1‑112 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
TGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGG  >  NC_002947/5741373‑5741654

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