Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,081,045 T→A 16.8% intergenic (+36/+26) PP_5332 → / ← PP_5333 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,081,0450TA16.8% 83.0 / 11.2 37intergenic (+36/+26)PP_5332/PP_5333hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base A (3/3);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAG  >  NC_002947/6080902‑6081191
                                                                                                                                               |                                                                                                                                                   
cAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                               <  1:469540/150‑1 (MQ=255)
     aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          <  1:444662/150‑1 (MQ=255)
                          gAATGTGTGGGATGACCGGGTGGAGTTGAAGCGTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAgg                                                                                                                     <  2:379795/150‑1 (MQ=255)
                           aaTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGt                                                                                                                    >  2:462844/1‑150 (MQ=255)
                             tgtgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCATAGGGTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCt                                                                                                                  <  1:7697/150‑1 (MQ=255)
                               tgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                >  2:573801/1‑150 (MQ=255)
                               tgtgGGATGACCTGGTGGAGTCGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                >  2:325139/1‑150 (MQ=255)
                                   ggATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTTTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCt                                                                                                            <  1:256199/150‑1 (MQ=255)
                                   ggATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCt                                                                                                            <  1:45224/150‑1 (MQ=255)
                                   ggATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCt                                                                                                            >  1:721162/1‑150 (MQ=255)
                                          cTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGCGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTccacca                                                                                                     >  2:389708/1‑150 (MQ=255)
                                              tggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCt                                                                                                 >  2:445976/1‑150 (MQ=255)
                                                    tGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTTCTCTTTTTTGACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCtt                                                                                           <  1:13522/150‑1 (MQ=255)
                                                     gAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTg                                                                                          >  1:747453/1‑150 (MQ=255)
                                                                    tCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCa                                                                           <  1:325139/150‑1 (MQ=255)
                                                                                      gATCGGTTGTGAAAATATCTGAGCCCTGCTGCATTATTGGTCGCTTCTCCGTGAGATCGGCGGCGCACCCGTCTGCGCAATTATCAAAGCTCTGCGTCTTCCACGACTCTCACCTTGTCTGCACCCAGCGCATAGGCCGCATCGGCCAgg                                                         >  1:14477/1‑150 (MQ=255)
                                                                                       aTCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                             <  2:260965/65‑1 (MQ=37)
                                                                                       aTCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                             >  1:260965/1‑65 (MQ=37)
                                                                                           gTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCAtt                                                   >  2:413652/1‑150 (MQ=255)
                                                                                             tgtgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGc                                                 <  2:144706/150‑1 (MQ=255)
                                                                                              gtgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCt                                                <  1:462844/150‑1 (MQ=255)
                                                                                                    atatCTGAGCCCCTCTCTTTTTGTACCGGCCTCTTCGTGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCtt                                          <  2:747453/150‑1 (MQ=255)
                                                                                                      atCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtct                                        <  2:721162/150‑1 (MQ=255)
                                                                                                       tCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGATGGGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCATCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtctt                                       >  2:754250/1‑149 (MQ=255)
                                                                                                       tCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtctc                                       >  1:9898/1‑150 (MQ=255)
                                                                                                         tGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCa                                     >  2:731897/1‑150 (MQ=255)
                                                                                                             cccTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtctc                                       >  1:683383/1‑144 (MQ=255)
                                                                                                             cccTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtctc                                       <  2:683383/144‑1 (MQ=255)
                                                                                                                           tACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGa                   <  2:14477/150‑1 (MQ=255)
                                                                                                                             ccGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTGTCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAcc                 >  2:578457/1‑150 (MQ=255)
                                                                                                                                  ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAg            >  1:188263/1‑150 (MQ=255)
                                                                                                                                   tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc           >  2:96645/1‑150 (MQ=255)
                                                                                                                                    ctTCGCGGGCAAGCCCTCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg          >  2:797059/1‑150 (MQ=255)
                                                                                                                                     ttCGCGGGCAAGCCCGCGAAGAGGCCGGTGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGcc                                                                     >  1:262383/1‑91 (MQ=38)
                                                                                                                                     ttCGCGGGCAAGCCCGCGAAGAGGCCGGTGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGcc                                                                     <  2:262383/91‑1 (MQ=255)
                                                                                                                                         cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgtg     <  1:573801/150‑1 (MQ=255)
                                                                                                                                            gCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCTCCTTGTCTGCATCCAGCGAATCGGCAGCATCGGCGTATTCACGTCCTGCGTTGTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAg  <  1:413652/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAG  >  NC_002947/6080902‑6081191

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: