Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,757,731 C→A 7.5% L263I (CTA→ATA)  PP_1563 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,757,7310CA 7.5% 164.8 / 3.7 53L263I (CTA→ATA) PP_1563terminase large subunit
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base A (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TACGACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGC  >  NC_002947/1757584‑1757878
                                                                                                                                                   |                                                                                                                                                   
tACGACACGATGGAGACCGGGATGGGGGCCCACGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGAcc                                                                                                                                                     >  2:158039/1‑148 (MQ=255)
    acacGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGTTACGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACATATgg                                                                                                                                                >  1:266730/1‑149 (MQ=255)
        gATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccgg                                                                                                                                            <  5:169941/149‑1 (MQ=255)
          tGGCGACCGGGTTGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                          <  6:226943/149‑1 (MQ=255)
          tGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                          <  4:95870/149‑1 (MQ=255)
            gagaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTAGAGGGTGTGATCGACCGGCCGGCCATAGgtccggccg                                                                                                                                         <  2:12555/148‑2 (MQ=255)
             agaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCctc                                                                                                                                       >  5:195947/1‑149 (MQ=255)
             agaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCctc                                                                                                                                       <  5:316232/149‑1 (MQ=255)
              gaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCctct                                                                                                                                      <  6:195947/149‑1 (MQ=255)
                 cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTaca                                                                                                                                   <  7:14107/149‑1 (MQ=255)
                 cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTac                                                                                                                                    >  8:14107/1‑148 (MQ=255)
                     aTGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTAt                                                                                                                               <  1:320738/149‑1 (MQ=255)
                             ccGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAgg                                                                                                                       >  6:213119/1‑149 (MQ=255)
                               gcgAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGc                                                                                                                     >  5:291155/1‑149 (MQ=255)
                                          atcagGCTGGTGATCACCACGTCGTGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCTGCCCTCTACACTATCGACCAGGGCGACGATTGGAc                                                                                                          <  6:312531/149‑1 (MQ=255)
                                             atGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCa                                                                                                        <  2:54882/148‑1 (MQ=255)
                                              tGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCAGCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGTACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGc                                                                                                      <  5:274985/149‑1 (MQ=255)
                                                   gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCTACCGGCCGGACCTATGGCCGGCCCTCTACTCTATCGACCAGGGCGACGATTGGACCAGCGaaa                                                                                                  <  1:39876/148‑1 (MQ=255)
                                                      aTCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCg                                                                                               <  6:214734/148‑1 (MQ=255)
                                                      aTCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCg                                                                                               >  3:293872/1‑148 (MQ=255)
                                                       tCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCa                                                                                             >  5:103638/1‑149 (MQ=255)
                                                         accacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccg                                                                                                                                             <  5:220134/99‑1 (MQ=255)
                                                         accacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccg                                                                                                                                             >  6:220134/1‑99 (MQ=255)
                                                          ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTg                                                                                          <  5:213119/149‑1 (MQ=255)
                                                                     tCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGaa                                                                               >  1:139529/1‑149 (MQ=255)
                                                                     tCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGaa                                                                               >  5:103401/1‑149 (MQ=255)
                                                                       aaGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCaa                                                                                     >  6:89066/1‑141 (MQ=255)
                                                                       aaGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCaa                                                                                     <  5:89066/141‑1 (MQ=255)
                                                                       aaGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACATATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAAt                                                                              <  2:266730/148‑1 (MQ=255)
                                                                       aaGTATTGGCGGCCCCTGCCACCAGGTTCTCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTTGACCAGCGAAATCGCACTGCGCAAGGCGAAt                                                                              <  2:70972/148‑1 (MQ=255)
                                                                        aGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATcc                                                                            <  6:291155/149‑1 (MQ=255)
                                                                                cggcCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                                    <  6:103401/149‑1 (MQ=255)
                                                                                   ccccTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGCCCGGCCGGACCTATGGCCGGCCCTCTCCACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGa                                                                 >  5:314053/1‑149 (MQ=255)
                                                                                     ccTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCATCGAAATCGCACTCCGCAAGTCGAATCCAAACTTTgttat                                                                >  8:228102/1‑144 (MQ=255)
                                                                                     ccTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAATGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTAGGCata                                                               >  6:317628/1‑149 (MQ=255)
                                                                                           accaGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTc                                                         >  8:287461/1‑149 (MQ=255)
                                                                                                   aTCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                                    <  1:106501/130‑1 (MQ=255)
                                                                                                   aTCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                                    >  2:106501/1‑130 (MQ=255)
                                                                                                       gcgATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                                   <  8:231860/127‑1 (MQ=255)
                                                                                                       gcgATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                                   >  7:231860/1‑127 (MQ=255)
                                                                                                       gcgATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                                   >  2:164034/1‑127 (MQ=255)
                                                                                                       gcgATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                                   <  1:164034/127‑1 (MQ=255)
                                                                                                              cGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATc                                                                                                <  4:52146/91‑1 (MQ=255)
                                                                                                              cGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATc                                                                                                >  3:52146/1‑91 (MQ=255)
                                                                                                              cGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTg                                      >  1:310873/1‑149 (MQ=255)
                                                                                                                  cGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCg                                  >  8:15920/1‑149 (MQ=255)
                                                                                                                        cTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACGCTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCg                            <  1:158039/149‑1 (MQ=255)
                                                                                                                            tgggTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGACATCGCACTGCGCTAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgat                        <  5:317628/148‑1 (MQ=255)
                                                                                                                              ggTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgc                      >  8:18348/1‑149 (MQ=255)
                                                                                                                                     tCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCa                <  2:310873/148‑1 (MQ=255)
                                                                                                                                        gagggccggccATAGGTCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCg            >  2:151711/4‑149 (MQ=255)
                                                                                                                                        aCCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTga                         <  3:138981/136‑1 (MQ=255)
                                                                                                                                        aCCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTga                         >  4:138981/1‑136 (MQ=255)
                                                                                                                                                 aCCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAg   <  6:314053/149‑1 (MQ=255)
                                                                                                                                                  ccTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGc  <  6:103638/149‑1 (MQ=255)
                                                                                                                                                   cTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAg   >  3:107027/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACGACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGC  >  NC_002947/1757584‑1757878

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