Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,764,227 A→G 8.5% E45E (GAA→GAG)  PP_1573 → major tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,764,2270AG 8.5% 130.6 / 2.5 47E45E (GAA→GAG) PP_1573major tail protein
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base G (2/2);  total (23/24)
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.96e-01

AATTGGAGAAGCCCATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCGAGAGCAT  >  NC_002947/1764079‑1764373
                                                                                                                                                    |                                                                                                                                                  
aaTTGGAGAAGCCCATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGaa                                                                                                                                                    <  3:47532/149‑1 (MQ=255)
aaTTGGAGAAGCACATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGaa                                                                                                                                                    <  8:265743/149‑1 (MQ=255)
     gagaAGCCCATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGa                                                                                                                                                <  8:230246/148‑1 (MQ=255)
     gagaAGCCCATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGAc                                                                                                                                               <  3:205003/149‑1 (MQ=255)
           cccATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGtct                                                                                                                                         <  5:67709/149‑1 (MQ=255)
              aTGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAg                                                                                                                                      <  5:26297/149‑1 (MQ=255)
                          gcgcAGCGTAGCCAGGGCTCTGCGCTGGTGCCGACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCTCCGCTTTCAGTCGAGGCGGGATCGCCGCTGACCAGGTTGAAGTGACCTGTGTCAGGGACACCGttg                                                                                                                           <  8:103275/145‑3 (MQ=255)
                           cTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCg                                                                                                                              >  5:52663/1‑144 (MQ=255)
                           cTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCg                                                                                                                              <  6:52663/144‑1 (MQ=255)
                           cTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTgcg                                                                                                                          >  3:177999/1‑148 (MQ=255)
                              cgggTACCCAGCTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATACGAAATCATGGAGGTTGAATGCGACACCGCCTTCAGTCCAGGCGGGAACCCCGCTGCCCAGGTAGAAGAGACCTGACTCAGCGACACCGTGCGCAAg                                                                                                                      >  8:163129/2‑149 (MQ=255)
                                     ccAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGc                                                                                                               >  6:66993/1‑149 (MQ=255)
                                      cAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCt                                                                                                                                             >  1:122278/1‑118 (MQ=255)
                                      cAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCt                                                                                                                                             <  2:122278/118‑1 (MQ=255)
                                       aGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGt                                                                                                             >  1:255289/1‑149 (MQ=255)
                                       aGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGt                                                                                                             <  4:177999/149‑1 (MQ=255)
                                           cTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCaa                                                                                                                       >  2:308439/1‑135 (MQ=255)
                                           cTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCaa                                                                                                                       <  1:308439/135‑1 (MQ=255)
                                           cTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGc                                                                                                          >  1:126217/1‑148 (MQ=255)
                                              tGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGaa                                                                                                                                                    >  7:281916/1‑103 (MQ=255)
                                              tGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGaa                                                                                                                                                    <  8:281916/103‑1 (MQ=255)
                                              tGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGcc                                                                                                         <  6:28442/146‑1 (MQ=255)
                                                     gTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCAcc                                                                                                 <  8:34138/147‑1 (MQ=255)
                                                        cccACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCggg                                                                                            >  2:294594/1‑149 (MQ=255)
                                                                  cAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCtt                                                                                    <  2:255289/147‑1 (MQ=255)
                                                                  cAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTc                                                                                   <  5:66993/148‑1 (MQ=255)
                                                                         taccgtaTTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCATCGACACCGTGCGCAAGTACCTGCGTGTCCTTCGCACCCCGGGCCAGGCTTCGCTGAccc                                                                           <  7:163129/149‑1 (MQ=255)
                                                                            cgtaTTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTc                                                                         >  4:135754/1‑148 (MQ=255)
                                                                            cgtaTTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGAGCAAGTACCTGCGTGGCCGTCGCACCCCGGGCCAGGCTACGCTGACCCTc                                                                         >  1:125352/1‑148 (MQ=255)
                                                                               aTTCGAAGTCATGGAGGTTGATTGCGCCACCGCCTTCAGTCCAGGGGGGAACCCCGCTGACCTGGTTGAAGTGACCTGTCTCAGGGACACCGTGCGCATGTACCTGCGTGGGCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCa                                                                        <  1:70742/146‑1 (MQ=255)
                                                                               aTTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACg                                                                     >  1:156653/1‑149 (MQ=255)
                                                                               aTTCGAACTCATGGACGATGAATCCGCCACCCCCTTCAGTCCACGCCGGAACCCCCCTGACCAGGTTGAAGTCACCTGACTCAGCGACCCCCTGCGCAAGTACCAGCGTCGACATCGCACCCCGCGCCAGGCTTCGCTGACCCTCa                                                                        >  2:70742/1‑146 (MQ=255)
                                                                                             aGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAAc                                                       >  6:211985/1‑149 (MQ=255)
                                                                                                   aaTGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCg                                                 <  1:294594/149‑1 (MQ=255)
                                                                                                           aCCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGAc                                                                             <  4:258123/113‑1 (MQ=255)
                                                                                                                        cAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTgc                                          >  8:27007/1‑135 (MQ=255)
                                                                                                                        cAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTgc                                          <  7:27007/135‑1 (MQ=255)
                                                                                                                           gCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTtcatca                                <  6:170775/142‑1 (MQ=255)
                                                                                                                           gCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTtcatca                                >  5:170775/1‑142 (MQ=255)
                                                                                                                               gAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCgg                        >  8:255870/1‑146 (MQ=255)
                                                                                                                                aaCCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAga                    >  6:11671/1‑149 (MQ=255)
                                                                                                                                 aCCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAgac                   <  5:11671/149‑1 (MQ=255)
                                                                                                                                   cccGCTGACCAGGTTGAGGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTc                                                                   >  4:156989/1‑99 (MQ=255)
                                                                                                                                   cccGCTGACCAGGTTGAGGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTc                                                                   <  1:160387/99‑1 (MQ=255)
                                                                                                                                   cccGCTGACCAGGTTGAGGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTc                                                                   <  3:156989/99‑1 (MQ=255)
                                                                                                                                   cccGCTGACCAGGTTGAGGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTc                                                                   >  2:160387/1‑99 (MQ=255)
                                                                                                                                     cGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGaca               >  3:215133/1‑149 (MQ=255)
                                                                                                                                              ggTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCgaga      <  7:255870/149‑1 (MQ=255)
                                                                                                                                              ggTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCgaga      <  5:211985/149‑1 (MQ=255)
                                                                                                                                                 tGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCGAGAGCa   >  7:102399/1‑149 (MQ=255)
                                                                                                                                                  gAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCGAGAGCAt  >  5:239868/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATTGGAGAAGCCCATGGCGATTCTCGCTCAGGGTACCCAGGTCTATGCCCTGGTGCCCACAGCTGCAAACCGTACCGTATTCGAAGTCATGGAGGTTGAATGCGCCACCGCCTTCAGTCCAGGCGGGAACCCCGCTGACCAGGTTGAAGTGACCTGTCTCAGCGACACCGTGCGCAAGTACCTGCGTGGCCTTCGCACCCCGGGCCAGGCTTCGCTGACCCTCAACGTCGACCCACGCAACGCCTCGCACGTGCGCCTTCATCAGATCTCGGAAGACGACACCATCGAGAGCAT  >  NC_002947/1764079‑1764373

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