Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,402,030 A→C 14.0% L375R (CTG→CGG) ‡ PP_3881 ← terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,402,0300AC14.0% 116.8 / 12.1 49L375R (CTG→CGG) ‡PP_3881terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base C (4/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

CACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCT  >  NC_002947/4401883‑4402178
                                                                                                                                                   |                                                                                                                                                    
cACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                     <  2:269374/149‑1 (MQ=255)
   gccgcGTAAGGCACCTTGTCGCGGGCTTCGCCCTTCTCGAGCTCGTGGTCGGGTATCCACAAATACGGCAACAGCCGCCAGTGCGGATCCACGAATCTCGGGTAGAACAGCAGGACAAATAAAGTAAGCTCCCTCGTACAGGAAAGGTc                                                                                                                                                  >  4:68920/1‑149 (MQ=255)
           aGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACCTGTCCAGTACga                                                                                                                                          <  5:160712/149‑2 (MQ=255)
           aGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGGACTGGACCTGTCCAGTACga                                                                                                                                          <  5:37690/149‑2 (MQ=255)
                    gTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTg                                                                                                                                  <  8:356675/148‑1 (MQ=255)
                    gTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTg                                                                                                                                  <  4:177974/148‑1 (MQ=255)
                             ttCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTgc                                                                                                                        >  8:315783/1‑149 (MQ=255)
                                      gTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGc                                                                                                                <  3:118756/148‑1 (MQ=255)
                                      gTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGc                                                                                                                <  4:191365/148‑1 (MQ=255)
                                      gTCGAGCTCGTGGGCGGGAATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGCCAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGc                                                                                                                <  5:365936/148‑1 (MQ=255)
                                         gAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACa                                                                                                            >  8:107807/1‑149 (MQ=255)
                                             tCGTGGTCGGGTATCCAGAAATACTGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCAAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCgg                                                                                                        >  2:253118/1‑149 (MQ=255)
                                             tCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCgg                                                                                                        >  2:219867/1‑149 (MQ=255)
                                               gTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTa                                                                                                      <  7:351325/149‑1 (MQ=255)
                                               gTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTa                                                                                                      >  8:351325/1‑149 (MQ=255)
                                               gTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGAGCACCAGCGACCTCGGTa                                                                                                      >  8:236176/1‑149 (MQ=255)
                                                   tCGGGTATCCAGAATTACGGCAACAGCCGCCAGGGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACCGGTGCAGTCCGCCAACGCTGGGACGATTGCGCAAGAGCGACGTAGGTAggcg                                                                                                  >  4:283200/1‑149 (MQ=255)
                                                      ggTATCCAGAAATACGGCAACAGCCGCCAGTGGGGATCCTCGTATGTCGGGTAGTACAGCAGGACAATTGAAGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGttc                                                                                               <  1:119576/149‑1 (MQ=255)
                                                      ggTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATCATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGttc                                                                                               <  3:68920/149‑1 (MQ=255)
                                                            cAGAAATACGTCAACAGCCGCCAGTGCGGATCCTGGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGc                                                                                         <  6:383654/149‑1 (MQ=255)
                                                              gAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCt                                                                                       <  3:169350/149‑1 (MQ=255)
                                                              gAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCt                                                                                       >  1:322911/1‑149 (MQ=255)
                                                                aaTACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCATCGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTgg                                                                                     <  1:93594/149‑1 (MQ=255)
                                                                     ggCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTc                                                                                  >  3:297252/1‑147 (MQ=255)
                                                                      gCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAg                                                                               <  2:59392/149‑1 (MQ=255)
                                                                                ccAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAgg                                                                     <  3:283200/149‑1 (MQ=255)
                                                                                   gTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAg                                                                      >  3:183368/1‑145 (MQ=255)
                                                                                         aTCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACAt                                                                           >  1:143321/1‑134 (MQ=255)
                                                                                         aTCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACAt                                                                           <  2:143321/134‑1 (MQ=255)
                                                                                          tCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACgg                                                           <  7:107807/149‑1 (MQ=255)
                                                                                                      gggTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCAc                                                <  7:315783/148‑1 (MQ=255)
                                                                                                      gggTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGACTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCAc                                                >  4:356542/1‑148 (MQ=255)
                                                                                                         tAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGt                                            >  4:364101/1‑149 (MQ=255)
                                                                                                           caaCAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCCCCAACACTTGGACGATTGCGCAACAGCGACATCTCTACTCGTTCTAAAACCAGCCACCATACATCCCACGAGATCCACGGCGACTCAGCACCAGTc                                           >  6:117540/2‑148 (MQ=255)
                                                                                                           caaCAGCAGGACAAATGATCTCAGGTCCGTCGTACTGGACAGGTCCACTCCGCCAACCCTGGGACGATTGCGCAACACCGACTACGCTACGCGTTCTACACCCTGGCTCCAGACATCCCAGGAGATCCCCGGCCACTCCCCACCAGAcc                                          >  6:112374/2‑149 (MQ=255)
                                                                                                                 cagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCg                                     >  6:286133/1‑148 (MQ=255)
                                                                                                                       aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGt                              >  3:74631/1‑149 (MQ=255)
                                                                                                                                    tCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCAc                 >  4:80689/1‑149 (MQ=255)
                                                                                                                                       gtcgtAGTGGAGAGGTCCAGTCCGCCTACGCTGGGACGATTGCGCAACAGCGAGATCGGTACGCGTTGTTCAGCGTGGCTCCAGACATCCCAGGAGATCCTCGGCGACTCAGCAGCAGTGCATTCGCAGAAGTTGAGGCGGCGCACCAt              <  5:112374/149‑1 (MQ=255)
                                                                                                                                       gtcgtACTGGACAGGTCCTGTCCGCCAACGCTGGGACGATTGCGAAACAGCGACATGGGTTCGCGTTCTTCAGCCTGGCTCCAGACATCCGAGGTGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCAt              <  5:117540/149‑1 (MQ=255)
                                                                                                                                       gtcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGCCGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACGAt              >  5:197748/1‑149 (MQ=255)
                                                                                                                                       gtcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCAt              >  5:139774/1‑149 (MQ=255)
                                                                                                                                        tcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTc                                      >  6:112512/1‑124 (MQ=255)
                                                                                                                                        tcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTc                                      <  5:112512/124‑1 (MQ=255)
                                                                                                                                         cggaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATgg            >  2:82374/4‑149 (MQ=255)
                                                                                                                                          ggaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGc           >  6:168007/3‑149 (MQ=255)
                                                                                                                                           gaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCa                                              <  2:295836/112‑1 (MQ=255)
                                                                                                                                           gaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCa                                              >  1:295836/2‑113 (MQ=255)
                                                                                                                                                gACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCTGCACCATGGCTTCCtt     >  2:102978/1‑149 (MQ=255)
                                                                                                                                                   aGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCt  >  1:60312/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCT  >  NC_002947/4401883‑4402178

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