Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,989 A→C 10.8% noncoding (91/109 nt) spf → small regulatory RNA Spot 42

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,9890AC10.8% 77.8 / 2.6 35noncoding (91/109 nt)spfsmall regulatory RNA Spot 42
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base C (2/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

AATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAA  >  NC_000913/4049848‑4050127
                                                                                                                                             |                                                                                                                                          
aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCTGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCGCCAGTCAGTAAt                                                                                                                                           >  2:305492/1‑143 (MQ=255)
   tACAAAAAGTGCGTTCAAAAAAAAAAAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                                                        <  2:265050/143‑1 (MQ=255)
     cAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTCATGACTg                                                                                                                                      >  2:335847/1‑143 (MQ=255)
      aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                     >  2:33526/1‑143 (MQ=255)
              cTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                                                             <  1:620762/143‑1 (MQ=255)
                ttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                           >  2:393183/1‑143 (MQ=255)
                     aactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTCATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                      >  1:313388/1‑143 (MQ=255)
                               aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTATCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGTGCGTTTTTTCTTAGGCgacagaaa                                                                                                            >  2:65827/1‑143 (MQ=255)
                                    agTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                                                                       <  2:76755/143‑1 (MQ=255)
                                         agtgagtgGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                                                                  <  2:145616/143‑1 (MQ=255)
                                         aGTTAGTCGCGTATGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                                                                  <  2:621659/143‑1 (MQ=255)
                                             aGTCGCGGAGGGTCCAGAGGTAAGATGTCCTATCTGTCAGCCCGTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                                                              <  1:15829/143‑1 (MQ=255)
                                                 gcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCACAAGAAAAGATCCGTACTGCCTGag                                                                                          >  2:512619/1‑142 (MQ=255)
                                                      gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCGTGATGCGCt                                                                                     >  2:235212/1‑143 (MQ=255)
                                                      gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                                                                     >  2:494615/1‑143 (MQ=255)
                                                         ttcAGAGGTAAGTTGTTCTATCTTTCAGACTTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTTGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                                                                  <  2:561787/141‑1 (MQ=255)
                                                            agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCAGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAACCGTCATGCCTGATGCGCGATGttt                                                                               >  2:89807/1‑143 (MQ=255)
                                                                 tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                          <  2:25902/143‑1 (MQ=255)
                                                                 tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                          <  1:95141/143‑1 (MQ=255)
                                                                             tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTGTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                                                              <  2:512066/143‑1 (MQ=255)
                                                                                          ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt                                                 <  2:365305/143‑1 (MQ=255)
                                                                                             cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa                                              <  1:33526/143‑1 (MQ=255)
                                                                                                 aCGTAAGAGGTTTAGTCTGCATTTTCTAGCCTCCCCAGTCAGTAATGACTGGGGCATTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTa                                          <  2:435914/143‑1 (MQ=255)
                                                                                                   gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAACAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACAGTACAATTGTAAACTATTGAATTTAcc                                        >  1:328003/1‑142 (MQ=255)
                                                                                                          gATTTGGCTGAAGAGTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGa                                 <  1:372067/143‑1 (MQ=255)
                                                                                                           aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGAt                                >  2:289193/1‑143 (MQ=255)
                                                                                                           aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGAt                                >  2:434069/1‑143 (MQ=255)
                                                                                                                 cTGGATATTTTAGCCGCCCCGGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCca                          <  2:238762/143‑1 (MQ=255)
                                                                                                                       atTTTATCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                    <  2:230303/143‑1 (MQ=255)
                                                                                                                                 ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGTATTTACTGGGCGGATACGCGACGTCCCGCAGAAggct          >  2:244280/1‑140 (MQ=255)
                                                                                                                                  cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCCGGCCAACGGTAGAATTGTAATCTATAGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCAtt         >  1:232801/1‑143 (MQ=255)
                                                                                                                                   ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTa        >  1:268576/1‑143 (MQ=255)
                                                                                                                                   ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAGTTTACGGGCCGGATACGCCACATCCGGCACAGGCGTTa        >  2:390355/1‑143 (MQ=255)
                                                                                                                                     agtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg      <  2:268576/143‑1 (MQ=255)
                                                                                                                                         agtCATGACTGGGGCGTTTTTTATTGGCAGCGCGAAGAGATCCGTACTGCCAGATCCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTGACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa  <  1:140957/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
AATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAA  >  NC_000913/4049848‑4050127

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: