Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,601,181 A→G 13.1% 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,1810AG13.1% 87.9 / 5.6 38K220E (AAG→GAG) PP_2279head‑to‑tail joining protein
Reads supporting (aligned to +/- strand):  ref base A (21/12);  new base G (3/2);  total (24/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 = 7.50e-01

CTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAAGA  >  NC_002947/2601040‑2601329
                                                                                                                                             |                                                                                                                                                    
cTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCTAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTga                                                                                                                                               <  1:100906/149‑1 (MQ=255)
cTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTga                                                                                                                                               <  1:376055/149‑1 (MQ=255)
               aaGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGaaggccag                                                                                                                                >  2:89720/1‑143 (MQ=255)
                     gaggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCAAAACAAGAACGTCTCGATCTACACCAACATCACCCTCAAGGGTGACACCTGGAACGAGTATCagga                                                                                                                          >  1:98224/1‑149 (MQ=255)
                                          ccAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTca                                                                                                     >  8:373888/1‑148 (MQ=255)
                                           cAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTccc                                                                                                    >  2:53876/1‑149 (MQ=255)
                                                         tCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCAc                                                                                      <  1:53876/149‑1 (MQ=255)
                                                           aTTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCAc                                                                                      <  1:171356/147‑1 (MQ=255)
                                                           aTTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACCTGTATCAGGAGGTCAAGGGCCAGAATGTTCCCGGCTCCCGAGCCAc                                                                                      >  2:171356/1‑147 (MQ=255)
                                                                           gATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATAAGGAGTTCAAGGTCCAGATTGTTCCCGGCTCCCGAGGAACGTACCCGAAAGACAAGt                                                                     >  8:185733/1‑148 (MQ=255)
                                                                                cTCCGGCGACAACCCCCGAGACCATAACAACAACGTCTCGATCTACACCCACATCTCCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCt                                                               <  7:194197/149‑1 (MQ=255)
                                                                                   cGGCGACACCCCCCGAGACCATAACAAGAACGTCACGATCTACACCCACATCTCCGTCAAGCGTGAGACCTGGAACGTGTATCAGGACGTCAAGGGCCAGAATGTTCCCGGCTCCCGAGGGACGTACCCGAAAGCCAACTGCCCCTGGc                                                            >  1:161852/1‑149 (MQ=255)
                                                                                          acCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCt                                                               <  5:79291/139‑1 (MQ=255)
                                                                                          acCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCATGGCTCCagat                                                     >  7:320202/1‑149 (MQ=255)
                                                                                          acCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTACCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCt                                                               >  6:79291/1‑139 (MQ=255)
                                                                                              cccGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattcga                                                 >  2:351683/1‑149 (MQ=255)
                                                                                                         aaCATGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgatt                                                    <  5:415773/135‑1 (MQ=255)
                                                                                                         aaCAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCCGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgatt                                                    >  6:415773/1‑135 (MQ=255)
                                                                                                            aataaCGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAgg                                                                                                                  <  5:348922/67‑1 (MQ=255)
                                                                                                              gAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACg                                 >  7:308795/1‑149 (MQ=255)
                                                                                                                acttctaGATCTACACCCACATCACCCTCAAGGGTGAGACGTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGt                               <  8:320202/146‑1 (MQ=255)
                                                                                                                   tctcGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGt                               >  2:38223/1‑146 (MQ=255)
                                                                                                                    ctcGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgatt                                                    >  8:214365/1‑124 (MQ=255)
                                                                                                                    ctcGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgatt                                                    <  7:214365/124‑1 (MQ=255)
                                                                                                                              acCCACATCACACTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTAAAGGGACAGATTGTTCCCGTCTCCCGAGGCACGTACCCGAAAGACAAGTGAGCCTGGCTCCCGAATCGATTCGTCAAGATAGACGCTGAGAACTATGGcct                 >  4:334467/1‑148 (MQ=255)
                                                                                                                               cccACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAACACAAGTGCACCAGGCTCCCGATTCGAATCGTCAAGATAGACGGTGAGAACTATGTCCg                 >  8:72737/1‑148 (MQ=255)
                                                                                                                                ccACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCt               >  5:150423/1‑149 (MQ=255)
                                                                                                                                 ggtctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGcc                  <  8:93045/141‑1 (MQ=255)
                                                                                                                                 ggtctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGcc                  >  7:93045/5‑145 (MQ=255)
                                                                                                                                  gtctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTcc             >  5:181629/4‑149 (MQ=255)
                                                                                                                                   tctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCaa                                         >  6:197121/3‑120 (MQ=255)
                                                                                                                                   tctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCaa                                         <  5:197121/118‑1 (MQ=255)
                                                                                                                                    aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGc                   >  5:126266/1‑141 (MQ=255)
                                                                                                                                    aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGc                   <  6:126266/141‑1 (MQ=255)
                                                                                                                                    aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCt            >  6:194641/1‑148 (MQ=255)
                                                                                                                                    aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCt            <  5:194641/148‑1 (MQ=255)
                                                                                                                                    aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTa           >  4:110544/1‑149 (MQ=255)
                                                                                                                                           tCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTgg      >  8:61348/1‑147 (MQ=255)
                                                                                                                                             aaGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAaga  <  1:351683/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAAGA  >  NC_002947/2601040‑2601329

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