Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,699 G→T 14.6% P41Q (CCG→CAG)  parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,6990GT14.6% 126.3 / 11.0 48P41Q (CCG→CAG) parAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base T (3/4);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTCACGGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGC  >  NC_002947/1551‑1844
                                                                                                                                                    |                                                                                                                                                 
gtcaCGGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGCGTAGTGCATAGCCTGTGCCAAGTCGCACTCCCCGATGAGCAAGTCGTAGACCGAATGCTGGAACTCTTGCTTGTCCACACCGCTCCCCGTGGTGGCGTTGCCCTGCgg                                                                                                                                                  >  2:1181383/1‑150 (MQ=255)
gtcaCGGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCgg                                                                                                                                                  >  2:1124101/1‑150 (MQ=255)
  caCGGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGAt                                                                                                                                                <  2:243375/150‑1 (MQ=255)
       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)
                               ccGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATc                                                                                                                                               >  2:240852/1‑122 (MQ=255)
                               ccGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATc                                                                                                                                               <  1:240852/122‑1 (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                                                                                                                           <  1:422140/108‑2 (MQ=38)
                                                                 cGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAGGgcaa                                                                                                                           >  2:422140/1‑107 (MQ=38)
                                                                    aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                              <  2:843427/150‑1 (MQ=255)
                                                                         ccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                         <  1:863364/150‑1 (MQ=255)
                                                                         ccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                         <  2:109208/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)
                                                                                                                                                    |                                                                                                                                                 
GTCACGGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGC  >  NC_002947/1551‑1844

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: ‑