Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,429,067 C→T 12.9% S244S (AGC→AGT PP_3044 → lambda family portal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,429,0670CT12.9% 72.6 / 4.0 31S244S (AGC→AGTPP_3044lambda family portal protein
Reads supporting (aligned to +/- strand):  ref base C (12/15);  new base T (2/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CCCCAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCTT  >  NC_002947/3428920‑3429216
                                                                                                                                                   |                                                                                                                                                     
ccccAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCt                                                                                                                                                     <  2:82658/150‑1 (MQ=255)
               ttCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGAGAACTacgacga                                                                                                                                      >  1:57558/1‑150 (MQ=255)
               ttCGCTGAACGCCGGTTACAACCCGCTNGTGCGTGTCCCGNCCNATCNGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTacgacgg                                                                                                                                      >  1:457261/1‑149 (MQ=255)
                         gCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGTCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGAcgcgagacat                                                                                                                            <  2:279872/150‑7 (MQ=255)
                            ggTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCTCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTc                                                                                                                         <  1:54052/150‑1 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc                                                                                                                        <  1:279225/145‑1 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc                                                                                                                        <  1:279301/145‑1 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc                                                                                                                        >  2:279301/1‑145 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc                                                                                                                        >  2:279225/1‑145 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAg                                                                                                                   >  1:343537/1‑150 (MQ=255)
                                  aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAg                                                                                                                   >  1:343580/1‑150 (MQ=255)
                                                 gtCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTg                                                                                                    <  2:448255/150‑1 (MQ=255)
                                                     cGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCg                                                                                                <  1:249343/150‑1 (MQ=255)
                                                           aTCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCt                                                                                          <  2:343537/150‑1 (MQ=255)
                                                           aTCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCt                                                                                          <  2:343580/150‑1 (MQ=255)
                                                                aTGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCAAACTTGTTCGCCGGCTTTGtt                                                                                     >  2:124025/1‑150 (MQ=255)
                                                                 tgctgcACATCTTCGAGCCGTTGGAGCCTCGCGAATCCCGTGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTc                                                                                    <  2:40888/150‑1 (MQ=255)
                                                                    tgcACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCa                                                                                 >  2:156416/1‑150 (MQ=255)
                                                                                   cGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATg                                                                  <  1:124025/150‑1 (MQ=255)
                                                                                     ttGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCgg                                                                                            >  1:330846/1‑122 (MQ=255)
                                                                                     ttGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTCGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCgg                                                                                            <  2:330846/122‑1 (MQ=255)
                                                                                             tGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAACTGGCCACCTTGTTAGACGGCTTTCTATACAGACCTGCGCCGGATGGGCTGGGTgg                                                        >  2:310303/1‑150 (MQ=255)
                                                                                                      gcccggTGGCCCCCGGCTTGCCCCGATCCTCTATTGCCTGCCCATCCTGGAAACCTACGAAGCCCCCTATAACTACCGGCAGTAAGTGGCCAACTGGTGCGCCGGCTTTGTTCGCAAACCGGCGCCGGATCGGCTGGGTGGCCCGcccct                                               <  2:57558/146‑4 (MQ=255)
                                                                                                                  ccGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCAc                                   <  1:224159/150‑1 (MQ=255)
                                                                                                                        ggCGCCGATCCTTAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACcggcgg                             <  2:75489/150‑1 (MQ=255)
                                                                                                                           gcCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGc                          >  1:320782/1‑150 (MQ=255)
                                                                                                                              gATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGgt                       >  2:34027/1‑150 (MQ=255)
                                                                                                                               aTCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGgtt                      <  2:172040/150‑1 (MQ=255)
                                                                                                                                   tCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTTTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTc                  >  2:222983/1‑150 (MQ=255)
                                                                                                                                   tCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTc                  <  1:310303/150‑1 (MQ=255)
                                                                                                                                          ccTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCg           >  2:382988/1‑150 (MQ=255)
                                                                                                                                                   ccTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCtt  <  1:382988/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCCAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCTT  >  NC_002947/3428920‑3429216

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: