Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,601,180 C→T 13.1% L219L (CTC→CTT)  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,1800CT13.1% 90.9 / 5.8 38L219L (CTC→CTT) PP_2279head‑to‑tail joining protein
Reads supporting (aligned to +/- strand):  ref base C (21/12);  new base T (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 = 3.81e-01

AACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGG  >  NC_002947/2601034‑2601325
                                                                                                                                                  |                                                                                                                                                 
aaCGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTc                                                                                                                                                   >  2:376055/1‑147 (MQ=255)
      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)
                                                                                                                                                  |                                                                                                                                                 
AACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGG  >  NC_002947/2601034‑2601325

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