Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,166,707 C→A 10.2% intergenic (‑551/‑230) PP_5483 ← / → PP_1921 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,166,7070CA10.2% 146.8 / 6.9 49intergenic (‑551/‑230)PP_5483/PP_1921hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base A (2/3);  total (23/26)
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.28e-01

CCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCG  >  NC_002947/2166560‑2166854
                                                                                                                                                   |                                                                                                                                                   
cccTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCt                                                                                                                                                    >  2:45014/1‑149 (MQ=255)
      aTCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAgg                                                                                                                                              <  3:154514/149‑1 (MQ=255)
          gAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGc                                                                                                                                          <  8:165156/149‑1 (MQ=255)
            aaGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAg                                                                                                                                        <  1:45014/149‑1 (MQ=255)
            aaGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGCATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGGCGCTGGGTGAGTGAAGCTCAATCAACCTGGAAGCCAGGTTgatt                                                                                                                                        <  6:269592/149‑4 (MQ=255)
                  ctctCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGATACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgta                                                                                                                                  <  1:281324/149‑1 (MQ=255)
                   tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGAAGCCAGGTTgattgag                                                                                                                                     <  4:298934/145‑7 (MQ=255)
                   tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGAAGCCAGGTTgattgag                                                                                                                                     >  3:298934/1‑139 (MQ=255)
                     tcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                >  1:313627/1‑148 (MQ=255)
                     tcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtag                                                                                                                                 <  2:4684/147‑1 (MQ=255)
                       aGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                <  5:151193/146‑1 (MQ=255)
                        gCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagt                                                                                                                            <  6:323793/149‑1 (MQ=255)
                          aaCCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtag                                                                                                                          <  2:313627/149‑1 (MQ=255)
                          aaCCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtag                                                                                                                          >  4:150525/1‑149 (MQ=255)
                                  gCAGAGTGTTTAAGGATCGCCTTTGGCCATTATGTGAAACTACCCGCTCCTGGAGCGGTATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGATTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGTTCgtagggagtagtgagtag                                                                                                                   <  8:203290/148‑1 (MQ=255)
                                   cAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagg                                                                                                                  <  6:235285/148‑1 (MQ=255)
                                   cAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagg                                                                                                                  <  3:150525/148‑1 (MQ=255)
                                   cAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagg                                                                                                                  >  5:235285/1‑148 (MQ=255)
                                            tAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagtaggga                                                                                                         <  8:43780/148‑1 (MQ=255)
                                            tAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagtaggga                                                                                                         <  1:115770/148‑1 (MQ=255)
                                                  tCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAgggg                                                                                                  >  8:79851/1‑149 (MQ=255)
                                                         tGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGga                                                                                            <  8:6129/148‑1 (MQ=255)
                                                                 aTGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgcg                                                                                   >  1:232352/1‑149 (MQ=255)
                                                                   gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGcgcg                                                                                 >  7:279627/1‑149 (MQ=255)
                                                                             ccGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgc                                                                                    >  4:310298/1‑136 (MQ=255)
                                                                             ccGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgc                                                                                    <  3:310298/136‑1 (MQ=255)
                                                                                      gAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAg                                                                                               >  7:294047/1‑116 (MQ=255)
                                                                                            gatgatCCGAGGGCTGTGGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCa                                                        <  2:232352/149‑1 (MQ=255)
                                                                                             atgatCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGCAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAg                                                       <  1:45993/149‑1 (MQ=255)
                                                                                                    gAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGagagag                                                <  7:79851/149‑1 (MQ=255)
                                                                                                      gggCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGa                                                                <  4:172120/131‑1 (MQ=255)
                                                                                                      gggCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGa                                                                >  3:172120/1‑131 (MQ=255)
                                                                                                          tGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGt                                          >  2:104830/1‑149 (MQ=255)
                                                                                                                 gtgtCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGtt                                   >  6:304171/1‑149 (MQ=255)
                                                                                                                  tgtCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagta                                                                                                             >  8:225698/1‑74 (MQ=255)
                                                                                                                  tgtCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagta                                                                                                             <  7:225698/74‑1 (MQ=255)
                                                                                                                      gCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTAtttt                              >  4:102772/1‑149 (MQ=255)
                                                                                                                                       caatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATgcg             >  1:127387/1‑149 (MQ=255)
                                                                                                                                        ctgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTAtt                                >  7:91016/4‑129 (MQ=255)
                                                                                                                                        ctgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTAtt                                <  8:91016/126‑1 (MQ=255)
                                                                                                                                        aatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGcta                          <  2:284562/135‑1 (MQ=255)
                                                                                                                                        aatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGcta                          >  1:284562/1‑135 (MQ=255)
                                                                                                                                         atcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTAGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa           >  5:262613/1‑149 (MQ=255)
                                                                                                                                           caaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTAATATGACTATgcgc            >  7:137553/1‑146 (MQ=255)
                                                                                                                                             aCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGctacta                       <  4:239933/133‑1 (MQ=255)
                                                                                                                                             aCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGctacta                       >  3:239933/1‑133 (MQ=255)
                                                                                                                                              ccTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCacgacg      >  5:156539/1‑149 (MQ=255)
                                                                                                                                                  gCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGGAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTTTTTTGCTACTATGACTATGCGCACGACGGGCg  >  8:4005/1‑149 (MQ=255)
                                                                                                                                                  gCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCg  <  6:156539/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCG  >  NC_002947/2166560‑2166854

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