Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,887,120 C→A 8.0% intergenic (‑49/‑81) PP_5161 ← / → PP_5162 lysophospholipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,887,1200CA 8.0% 149.8 / 2.2 50intergenic (‑49/‑81)PP_5161/PP_5162lysophospholipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base A (2/2);  total (27/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 = 8.82e-01

GAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCT  >  NC_002947/5886981‑5887267
                                                                                                                                            |                                                                                                                                                   
gAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTc                                                                                                                                           >  1:108757/1‑150 (MQ=255)
     gCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCcaca                                                                                                                                      >  2:219183/1‑150 (MQ=255)
      cAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAg                                                                                                                                     <  2:768231/150‑1 (MQ=255)
       aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGTGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTcccgcaat                                                                                                                                    <  1:790645/150‑6 (MQ=38)
          gggCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGTGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTCCcgcaattag                                                                                                                                 >  2:519014/1‑145 (MQ=37)
           ggCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATc                                                                                                                                >  1:545066/1‑150 (MQ=255)
                   gTACATCCACATTGAGGCAAGTATCGATATGTAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAg                                                                                                                        <  2:108757/150‑1 (MQ=255)
                         ccACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGCACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTa                                                                                                                  <  2:613202/150‑1 (MQ=255)
                          cacaTTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGACGCCAGCGTCAGGGATCGGCGCCAGCTTTTaa                                                                                                                 >  2:108588/1‑149 (MQ=255)
                           acaTTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGa                                                                                                                <  2:545066/150‑1 (MQ=255)
                             aTTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCAACTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAgg                                                                                                                                    >  1:50244/1‑128 (MQ=255)
                             aTTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCAACTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAgg                                                                                                                                    <  2:50244/128‑1 (MQ=255)
                              ttGATGCAATTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCAGCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAt                                                                                                             <  1:338112/150‑1 (MQ=255)
                                                 ggAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                          >  2:482666/1‑150 (MQ=255)
                                                     tCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                             >  1:669388/1‑111 (MQ=255)
                                                     tCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                             <  2:669388/111‑1 (MQ=255)
                                                         cAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                          >  1:436031/1‑142 (MQ=255)
                                                         cAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                          <  2:436031/142‑1 (MQ=255)
                                                             tctGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAc                                                                              <  1:482666/150‑1 (MQ=255)
                                                             tctGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAc                                                                              >  2:592808/1‑150 (MQ=255)
                                                                 gCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGAc                                                                          >  1:654760/1‑150 (MQ=255)
                                                                   cGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAGTTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGActc                                                                        >  1:817442/1‑150 (MQ=255)
                                                                       gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTccc                                                                                                                                         <  1:822883/81‑1 (MQ=21)
                                                                       gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTccc                                                                                                                                         >  2:822883/1‑81 (MQ=21)
                                                                             gCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATa                                                              <  1:219183/150‑1 (MQ=255)
                                                                               aaGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTc                                                                                                                                           >  1:579134/1‑71 (MQ=255)
                                                                               aaGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTc                                                                                                                                           <  2:579134/71‑1 (MQ=255)
                                                                                       tCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGcc                                                           >  2:383262/1‑143 (MQ=255)
                                                                                       tCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGcc                                                           <  1:383262/143‑1 (MQ=255)
                                                                                       tCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTACTCCAGTAAGGGCCACCGACTCTTCATGAATAGCCGAatcat                                                    >  1:261965/1‑150 (MQ=37)
                                                                                       tCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACCCACGGCCGGTCTACCGCGCCACGGCCGGAGCCTACGGAAATTTAACAACGCTGTTTCGTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAatcat                                                    >  2:736481/1‑150 (MQ=255)
                                                                                         aCCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAg                                                                                                         >  2:390910/1‑95 (MQ=255)
                                                                                         aCCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAg                                                                                                         <  1:390910/95‑1 (MQ=255)
                                                                                         aCCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGAACACACCAAAACGTTGGAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCGCCGACTCTTCATGAATAGCCGAAtcatca                                                  >  2:839859/1‑150 (MQ=255)
                                                                                           cTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCCCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCaaa                                                >  2:327415/1‑150 (MQ=255)
                                                                                              gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCTAAGCATCAAACCg                                             >  1:390425/1‑150 (MQ=255)
                                                                                                  agccAGCGTCGAACTACTCGCGGGCGCGGCCTTTCGCGCCCCCAGGCCGTCCCAACAGGGATCGGCGCCAGCTTTTAGAAATTGGGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgt                                         <  2:817442/148‑1 (MQ=255)
                                                                                                     cAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCa                                      >  1:82664/1‑150 (MQ=255)
                                                                                                       gCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACg                                    <  2:492035/150‑1 (MQ=255)
                                                                                                         aTCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAg                                                             >  1:479787/1‑123 (MQ=255)
                                                                                                         aTCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAg                                                             <  2:479787/123‑1 (MQ=255)
                                                                                                         aTCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGt                                  >  1:335406/1‑150 (MQ=255)
                                                                                                            gAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCAt                               >  1:547120/1‑150 (MQ=255)
                                                                                                                  aTTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAatcat                                                    <  2:767219/123‑1 (MQ=255)
                                                                                                                  aTTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAatcat                                                    >  1:767219/1‑123 (MQ=255)
                                                                                                                   ttGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                                   <  1:327246/91‑1 (MQ=255)
                                                                                                                   ttGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                                   >  2:327246/1‑91 (MQ=255)
                                                                                                                            cGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgc               <  1:592808/150‑1 (MQ=255)
                                                                                                                                 tttCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCAgcgaccc          <  2:82664/150‑5 (MQ=255)
                                                                                                                                          cGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCctct  <  2:390425/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCT  >  NC_002947/5886981‑5887267

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: