Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707 G→C 12.8% D38E (GAC→GAG parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,7070GC12.8% 135.0 / 10.9 47D38E (GAC→GAGparAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base C (3/3);  total (25/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 = 9.75e-01

TTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGG  >  NC_002947/1558‑1851
                                                                                                                                                     |                                                                                                                                                
ttGGCTGGCAGTAACTGGGAGCGTCCGGGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGTCCTGCTGATCGAgg                                                                                                                                                  <  2:974706/150‑1 (MQ=255)
    ccGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGa                                                                                                                                              >  1:429086/1‑150 (MQ=255)
          gcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgca                                                                                                                                        <  2:600916/150‑1 (MQ=255)
           caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTCCTGGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagcag                                                                                                                                       <  1:751322/150‑1 (MQ=255)
             aCTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGa                                                                                                                                     <  1:850045/150‑1 (MQ=255)
                 gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAgggcaac                                                                                                                                 <  2:429086/150‑6 (MQ=255)
                       tCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                           >  1:843427/1‑150 (MQ=255)
                                 gAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagcc                                                                                                                 >  1:243689/1‑150 (MQ=255)
                                    tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagcag                                                                                                                                       <  2:573679/125‑1 (MQ=255)
                                    tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagcag                                                                                                                                       >  1:573679/1‑125 (MQ=255)
                                        gCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAgg                                                                                                          >  1:403584/1‑150 (MQ=255)
                                              ccTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCga                                                                                                    >  2:190919/1‑150 (MQ=255)
                                                  ggCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGCTGTCTAGCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAtt                                                                                                <  1:1124101/150‑1 (MQ=255)
                                                      aaGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATa                                                                                            <  1:190919/150‑1 (MQ=255)
                                                         tCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                         >  2:575770/1‑150 (MQ=255)
                                                          cGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAGGgcaa                                                                                                                                  >  2:422140/1‑107 (MQ=38)
                                                          cGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAGGgcaa                                                                                                                                  <  1:422140/108‑2 (MQ=38)
                                                             aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                                     <  2:843427/150‑1 (MQ=255)
                                                                  ccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                                <  2:109208/150‑1 (MQ=255)
                                                                  ccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                                <  1:863364/150‑1 (MQ=255)
                                                                   cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCt                                                                               <  2:243689/150‑1 (MQ=255)
                                                                     aTGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTg                                                                             <  1:1181383/150‑1 (MQ=255)
                                                                                 tAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGGTGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCt                                                                 >  1:259823/1‑150 (MQ=255)
                                                                                    cacGACTGCTAGCGCTCCTCCCTGTCCACACCGCTGCCCATTGTGGCGTTGCCCTGCGGATCGAGGTCGATTAGCAGGAAACGACGCTTGGTCGCAGCCAGCGAGTCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTc                                                              <  2:259823/148‑1 (MQ=255)
                                                                                            cTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgc                                                      <  2:730565/150‑1 (MQ=255)
                                                                                             tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgcg                                                     >  1:231819/1‑150 (MQ=255)
                                                                                                   tCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                         <  2:736673/108‑1 (MQ=255)
                                                                                                   tCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                         >  1:736673/1‑108 (MQ=255)
                                                                                                   tCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCg                                               >  1:1186676/1‑150 (MQ=255)
                                                                                                     gTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGaa                                             >  2:1169981/1‑150 (MQ=255)
                                                                                                       gCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCt                                                             >  2:103648/1‑132 (MQ=255)
                                                                                                       gCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCt                                                             <  1:103648/132‑1 (MQ=255)
                                                                                                         ttGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGa                                                    <  2:992577/139‑1 (MQ=255)
                                                                                                         ttGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGa                                                    >  1:992577/1‑139 (MQ=255)
                                                                                                         ttGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAcc                                         >  2:1190044/1‑150 (MQ=255)
                                                                                                           gTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCtt                                       >  2:928139/1‑150 (MQ=255)
                                                                                                               acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGcc                                   >  2:680041/1‑150 (MQ=255)
                                                                                                                 acCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAt                                 >  1:403625/1‑150 (MQ=255)
                                                                                                                  ccGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAtt                                >  2:387664/1‑150 (MQ=255)
                                                                                                                        cccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCg                          >  1:1100599/1‑150 (MQ=255)
                                                                                                                           ttggtggCTTTGGCATGCGGATCGATGTCGATCAGCAGGACACGACGTTTGTTCGCAGCCAGCGAGCCGGCGAGATTGACACAGGTGCTAGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAgt                       <  2:51375/149‑1 (MQ=255)
                                                                                                                           aTGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCAACACCACATTTATGGTTAGCGAGTGCGAATATCTTAGCCATTGTTGCGAgt                       >  2:435863/1‑150 (MQ=255)
                                                                                                                             ggtggCGTTGCCGTGCGGATCGAGGTCGATCAGCAGCACACGACGCATGGTCGCAGCCAGCGAGGCGGCCAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGTTTCGAGATTGCCAATACCTTAGCCATTGTTGCGAgtgt                     <  2:57375/150‑1 (MQ=255)
                                                                                                                                  tcgtgtccTGCTGATCGACCTCGATCCGCAGGGCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATGGATAAAGGTGGCCAGCTAGTTGTCATCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCa                >  1:875625/7‑150 (MQ=255)
                                                                                                                                         cTGCGGATCGAGGCCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGc         >  2:849706/1‑150 (MQ=255)
                                                                                                                                              gATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgc    <  2:1197859/150‑1 (MQ=255)
                                                                                                                                              gATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgc    <  1:1169981/150‑1 (MQ=255)
                                                                                                                                                tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgcgg  >  1:966329/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGG  >  NC_002947/1558‑1851

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: