Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,601,181 A→G 14.8% K220E (AAG→GAG)  PP_2279 → head‑to‑tail joining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,601,1810AG14.8% 59.5 / 3.6 27K220E (AAG→GAG) PP_2279head‑to‑tail joining protein
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

CTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAA  >  NC_002947/2601040‑2601327
                                                                                                                                             |                                                                                                                                                   
cTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTga                                                                                                                                              <  7:76654/149‑1 (MQ=255)
     cATCGTGGTCATGGAGGAGGTCTCCCTGGCGGTACTGCCAGATGAGGCCCTGTCATTAGTCGATCATGGCGATGACTCCGTCGATTTTTTTTTATATGATAATTAGAACTTCTAGATATACACCCACATCACCCTTGAGGGTgatgtgg                                                                                                                                         <  2:1713/149‑6 (MQ=255)
     cATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTTGAGGGTgatgtgg                                                                                                                                         <  8:283417/149‑6 (MQ=255)
             tCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACgtg                                                                                                                                 >  6:197282/1‑149 (MQ=255)
                      aggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCaggag                                                                                                                        <  5:197282/149‑1 (MQ=255)
                        gTCTCCCTGGCTGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGACTTGAGATGACTCCGTAGATTTTTTTTTTTATGATAAGTAGAACTTCTAGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACgt                                                                                                                                  <  4:250091/137‑1 (MQ=255)
                        gTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACAAACATAAACCTCAAGGGTGAGACCTGGAACgt                                                                                                                                  >  3:250091/1‑137 (MQ=255)
                                 gCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGTAACGTGTATCAGGAGCTCAAGGAcc                                                                                                              >  2:259230/1‑148 (MQ=255)
                                            agagGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTcc                                                                                                    <  3:70028/147‑1 (MQ=255)
                                            agagGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGCGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTcc                                                                                                    >  4:70028/1‑147 (MQ=255)
                                                      ctggcATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATCACAAGAACGTCTCGCTCTACACCCACATCACCCTCAAGGGAGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAgg                                                                                        <  7:263502/145‑1 (MQ=255)
                                                               gTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATAAGGAGGTCAAGGGACAGATTGTTCCCGGCTCCCGAGGCACGTAccc                                                                               >  5:209515/1‑149 (MQ=255)
                                                                   aGCCTGGCGATGACTCCGGCTTCACCCCCCGAGAACATAACAAGAACGTCTCGATCTACACCCACCTCAGCCTCAAGGGTGAGACCTGGAACGTGTATCAGTAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGaaa                                                                           <  1:259230/149‑1 (MQ=255)
                                                                               aCTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGcc                                                              >  6:186111/1‑149 (MQ=255)
                                                                                                     ccATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTc                                         >  2:79199/1‑148 (MQ=255)
                                                                                                        tAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCgattcgattc                                            <  3:245455/142‑1 (MQ=255)
                                                                                                        tAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCgattcgattc                                            >  4:245455/1‑142 (MQ=255)
                                                                                                           cAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACaa                                                                     <  5:60696/114‑1 (MQ=255)
                                                                                                            aaGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGCCCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCACGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCAGGCTCCCGCTTCGATTCGTCAAGATCGa                                 >  7:321248/1‑149 (MQ=255)
                                                                                                                   tctcGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTga                           <  1:79199/149‑1 (MQ=255)
                                                                                                                      cGATCTACACCCACATCACCCTCGAGGGTGAGACCTGGAACGTGTATCAGAAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAAc                       >  8:181026/1‑149 (MQ=255)
                                                                                                                      cGATCTACACCCACATCACCCTCGAGGGTGAGACCTGGAACGTGTATCAGAAGGTCAAGGGCCAGATTGTTACCGGCTCCCGACGCACGTACCCGAAA‑GACAACTGCGCCTGGCTCCAGATTAGATTCGTCAAGATCGACGGTGAGAAc                       >  8:175853/1‑149 (MQ=255)
                                                                                                                       gATCTACACACATATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACt                      <  6:209515/149‑1 (MQ=255)
                                                                                                                               cccACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGc              <  5:186111/149‑1 (MQ=255)
                                                                                                                                ccACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCgattcgattc                                            <  6:220905/118‑1 (MQ=255)
                                                                                                                                ccACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCt             >  6:325270/1‑149 (MQ=255)
                                                                                                                                          cTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTgg    >  3:172676/1‑148 (MQ=255)
                                                                                                                                             aaGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGaa  <  4:172676/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAA‑GACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAA  >  NC_002947/2601040‑2601327

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