Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,402,031 G→T 12.3% L375M (CTG→ATG) ‡ PP_3881 ← terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,402,0310GT12.3% 129.6 / 10.0 49L375M (CTG→ATG) ‡PP_3881terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base T (3/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCTC  >  NC_002947/4401883‑4402179
                                                                                                                                                    |                                                                                                                                                    
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)
                                                                                                                                                    ggTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCTc  <  4:183368/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTGCTC  >  NC_002947/4401883‑4402179

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: ‑