Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 64,794 T→C A329A (GCA→GCG polB ← DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091364,7940TC97.1% 112.5 / NA 34A329A (GCA→GCGpolBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/21);  minor base A (1/0);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTAT  >  NC_000913/64549‑65035
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
aGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAg                                                                                                                                                                                                                                               <  1:129805/251‑1 (MQ=255)
                           ccgccAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGATCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAAcggcg                                                                                                                                                                                                                    >  2:120790/1‑251 (MQ=255)
                                    atGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCCGTCACCGTCGCGCGTTCGAGGAAAAAAGGCATGATTTCAGGTTTATGGGAAATCTGCTTCACCAGCTCGCAATCTTTCAGGTTCTAAGTCGCCAGCGCCGGTGTATCTTCGGAGAAACGGCGGTACAAttt                                                                                                                                                                                                           >  2:208897/2‑251 (MQ=255)
                                              gcggcACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGa                                                                                                                                                                                                <  1:231923/251‑1 (MQ=255)
                                                             gaTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCg                                                                                                                                                                                 <  1:203293/251‑1 (MQ=255)
                                                                 aGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAg                                                                                                                                                                             <  1:246336/251‑1 (MQ=255)
                                                                     gcgACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAAATTCGGCCATTCGATCCCACGGGGTTTCGATaacatt                                                                                                                                                                         >  2:220305/1‑246 (MQ=255)
                                                                       gACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTc                                                                                                                                                                       <  2:52046/251‑1 (MQ=255)
                                                                        aCATAACCAGCGCGGTGCATACGAGGGAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCGCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTGTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTGCGTCCATTCGATCCCACGGGTTATCGATAGATTTTcc                                                                                                                                                                      <  2:94827/251‑1 (MQ=255)
                                                                        aCATAACAAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGACAATCTGCGTCACCAGCTCGCACTCTTTCAGGTTATCAGTAGCCCGCGCAGGTTTATCTTCGGCGAAACGGCGGGCCATTTCGTCCATTTTATCCCACGGGGTATCTATAGATTTTcc                                                                                                                                                                      >  2:53938/1‑251 (MQ=255)
                                                                           tAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCCTTCGATCCCACGGGTTATCGATAGATTTTCCTTc                                                                                                                                                                   >  1:151641/1‑251 (MQ=255)
                                                                                  cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAAt                                                                                                                                                            <  1:214615/251‑1 (MQ=255)
                                                                                  cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGGTATCGATAGATTTTCCTTCTCCTAAt                                                                                                                                                            >  1:82465/1‑251 (MQ=255)
                                                                                              gaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTTAATTTCGTCCATTCGATCCCACGGGTTATCGATAGAATTTCCTTCGCCTACTAACTCCTGAGCg                                                                                                                                                >  2:90344/3‑251 (MQ=255)
                                                                                                     aTAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcgggtt                                                                                                                                         >  1:226213/1‑246 (MQ=255)
                                                                                                       agagATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTc                                                                                                                                        <  1:148857/250‑1 (MQ=255)
                                                                                                        gagaTGACCAAACGCCGACACCGATTCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTAGATCGGTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCGTCGGCGAAACGGCGGTCGATTGCGTCCATTCGATCCCACGTGTGATCGATAGATTTTCCTTCGCCTAATAACGCCTGAGCgcagggttcta                                                                                                                                      <  2:178985/251‑11 (MQ=255)
                                                                                                                           aCCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAAtt                                                                                                                   <  1:166707/251‑1 (MQ=255)
                                                                                                                            ccgatcccccgTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTc                                                                                                                  <  1:70926/243‑1 (MQ=255)
                                                                                                                                                                            cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACTGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCATAATAACTCCTTAGCGTCGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGTATTTCAGCGCCTCGATACCTTCGATAATTAGCCGAcc                                                                   <  1:180493/250‑1 (MQ=255)
                                                                                                                                                                            cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGAcc                                                                   <  1:180482/250‑1 (MQ=255)
                                                                                                                                                                                 aaaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGc                                                             >  2:242759/1‑251 (MQ=255)
                                                                                                                                                                                  aaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGcc                                                            <  1:236712/251‑1 (MQ=255)
                                                                                                                                                                                    aaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTg                                                          <  2:209172/251‑1 (MQ=255)
                                                                                                                                                                                         cATGATTTCAGTTGTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTGAGCCTGGGCaa                                                     <  2:25253/251‑1 (MQ=255)
                                                                                                                                                                                            gATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCaaaaa                                                  <  1:242759/251‑1 (MQ=255)
                                                                                                                                                                                               ttCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAg                                                <  1:26468/250‑1 (MQ=255)
                                                                                                                                                                                                          ggAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGttttt                                     <  2:191660/250‑1 (MQ=255)
                                                                                                                                                                                                                     gTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCCATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTc                         >  1:106304/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGCGCTCGCGCCACTcc               >  1:204400/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTc          >  2:234046/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTcactgt     <  2:235439/251‑6 (MQ=255)
                                                                                                                                                                                                                                         cAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTcactgt     <  1:90426/251‑6 (MQ=255)
                                                                                                                                                                                                                                            gTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTATCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCCGCTcactgtttt  >  1:66544/1‑243 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
AGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTAT  >  NC_000913/64549‑65035

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