Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,044 Δ1 bp 14.7% intergenic (+130/+59) PP_5440 → / ← PP_t72 hypothetical protein/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,0440A.14.7% 88.8 / 12.0 34intergenic (+130/+59)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base . (3/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAA  >  NC_002947/343908‑344189
                                                                                                                                        |                                                                                                                                                  
cGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTc                                                                                                                                      <  2:92485/150‑1 (MQ=255)
                                 tGCGGTTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAAccacca                                                                                                     >  2:454721/1‑150 (MQ=255)
                                      gTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGc                                                                                                >  1:329119/1‑150 (MQ=255)
                                           tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTGGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTGGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                           >  1:583857/1‑150 (MQ=255)
                                              aTTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggt                                                                                        <  2:267601/150‑1 (MQ=255)
                                               ttGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACCAGCAAAACGCACATCGCGCGTTGCTTGGGGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtg                                                                                       >  1:73064/1‑150 (MQ=255)
                                                   aTGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                                   <  2:329119/150‑1 (MQ=255)
                                                   aTGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                                   <  1:65716/150‑1 (MQ=255)
                                                      gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                                    >  2:80337/1‑98 (MQ=255)
                                                      gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                                    <  1:80337/98‑1 (MQ=255)
                                                            cTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCCTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                          <  1:423052/150‑1 (MQ=255)
                                                              tCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                        <  1:454721/150‑1 (MQ=255)
                                                                  gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCg                                                                    >  2:238492/1‑150 (MQ=255)
                                                                             aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTAGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGaca                                                         >  1:333628/1‑150 (MQ=255)
                                                                                tctcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTAGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGa                                                      <  2:333628/150‑1 (MQ=255)
                                                                                    gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAt                                                  >  1:48009/1‑150 (MQ=255)
                                                                                      gtTCTTCGACGCCGCCTGAAGTTAAAACTATCAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAttt                                                >  1:475963/1‑150 (MQ=255)
                                                                                       ttcttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                             <  2:296289/88‑1 (MQ=255)
                                                                                       ttcttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                             >  1:296289/1‑88 (MQ=255)
                                                                                       ttcttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAgaga                                                                            <  1:593207/121‑1 (MQ=255)
                                                                                       ttcttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAgaga                                                                            >  2:593207/1‑121 (MQ=255)
                                                                                        tcttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAAAAAGCGCTGCTCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                              <  2:475963/150‑1 (MQ=255)
                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCg                                                                                                                                            >  2:323614/1‑54 (MQ=255)
                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCg                                                                                                                                            <  1:323614/54‑1 (MQ=255)
                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAg                                                                                                 <  2:142561/97‑1 (MQ=255)
                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAg                                                                                                 >  1:142561/1‑97 (MQ=255)
                                                                                            cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAGTTTCAATc                                          <  1:238492/150‑1 (MQ=255)
                                                                                                          gTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGa                            <  2:48009/150‑1 (MQ=255)
                                                                                                                        acaacaAGAAACCCGCCTAGTGCGGGGTTTTTGTTGTCGCTTCGGGTAACGTTGCAACCACCCGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg              >  2:604878/1‑150 (MQ=255)
                                                                                                                            caaGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac          >  2:165884/1‑150 (MQ=255)
                                                                                                                              agaaaCCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                              <  2:407616/126‑1 (MQ=39)
                                                                                                                              agaaaCCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                              >  1:407616/3‑128 (MQ=39)
                                                                                                                                 aaCCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGt     <  1:165884/150‑1 (MQ=255)
                                                                                                                                    ccGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACGCGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGaa  <  2:583857/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
CGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAA  >  NC_002947/343908‑344189

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