Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,392,567 G→T 19.6% intergenic (+20/‑16) glyX → / → glyY tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,5670GT19.6% 76.0 / 2.5 42intergenic (+20/‑16)glyX/glyYtRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (25/7);  new base T (7/2);  total (32/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTT  >  NC_000913/4392418‑4392716
                                                                                                                                                     |                                                                                                                                                     
gcccccgttCCGGCTCACAAATTCGAAAATTTCGGCACCAAACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAg                                                                                                                                                       <  2:343685/143‑1 (MQ=14)
       ttCCCGCTCCTAAATTTGAAAAGTGCTGCACAGAACATACCACCCAAGCGGTAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGa                                                                                                                                                <  1:418942/150‑1 (MQ=18)
        tCCCGCTCCAAAATTTGAAAAGTGCTTCACAGCACAGACCTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGac                                                                                                                                               <  1:1474261/150‑1 (MQ=255)
          ccGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGATAGTGCTGTAATGCACAgcacc                                                                                                                                             >  1:1468248/1‑146 (MQ=255)
            gCTCCAAAATTTGAAAAGTTCTGCAAAGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGaccacc                                                                                                                                           <  1:733641/150‑1 (MQ=255)
             cTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTTCTGTAAGGGACAGAGCAccc                                                                                                                                          >  1:1800492/1‑150 (MQ=255)
             cTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGAAATGGCTGAAAGGCACACaccgacc                                                                                                                                          >  1:489594/1‑146 (MQ=17)
             cTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAATAGCTATAATCAACAGACCGccc                                                                                                                                          >  1:1333735/1‑150 (MQ=14)
                 aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                                      >  1:1523154/1‑150 (MQ=255)
                       tGAAAAGTTCTGCAAACAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt                                                                                                                                <  2:1698109/150‑1 (MQ=255)
                        gAAAAGTGCTCCAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGAGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATa                                                                                                                               <  1:816508/150‑1 (MQ=255)
                                              ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGTGATGTAATGCACAGACCCCCCCAGCGGGAATAGCGCAGTTGGTATAGAACGAcc                                                                                                         >  2:1111526/1‑150 (MQ=14)
                                                                cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCACGCTCCAAAATATGAAAGTGATGTAAGACACACAACACCCAAGCGGGAATAGCTAAGTTTGTAGAGCAAGACAATGACAATGTCgtagacc                                                                                       >  1:1650894/1‑144 (MQ=11)
                                                                          tgagCCCTACCTTGCCAAGGTCTGGGTCGCGAGTTCGAGTCTCGTTTCCCGCCCCAAAATTTTAAAGTCCTTAAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                             <  1:1061870/149‑1 (MQ=11)
                                                                          agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTTTAATGAACAGACCACCCAAGAGGGAATATCTCAGTTGGTAGAGCACGACCATAACAAGGTAGTGGTCGCGAGTACGAg                                                                             >  2:569974/1‑150 (MQ=11)
                                                                             gcccgcccTTGCCAAGGTCGGGGTGGCGATTGCGAGTCTCGTTTCCCGCTCCAAATTTTTAAATTGCTGTAAGTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                                          <  2:1067112/144‑1 (MQ=11)
                                                                             gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGAGATTTAAGAAACAGACCACCCAGACGGGATTAGCTAAGTTTTTATAGAAAGACCTTGCCAAGGTCATGCTCTCTAGTTCGAGtca                                                                          >  2:1551421/1‑149 (MQ=2)
                                                                             gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGTTATGTAATGCACACACCCCCCAAGCGGGATTATCTCAGTTGTTAGAGCAACACCATGCCAATGTAGGGGTCGCGATTTCTAggca                                                                          >  2:1830180/1‑147 (MQ=11)
                                                                              cccGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGATGAAATCAACAGACCACCCACGCGGGAAAAGGACAGTTTGTCGAGCCCCACCATGCCAATGGACCGGTCGCGAGTACGACtccc                                                                         >  2:1223650/3‑148 (MQ=2)
                                                                                     ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGTAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCATGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                  >  2:1721719/1‑150 (MQ=255)
                                                                                        ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGTGCTGTAAGTAACAGACCACCCAAGCGGGAATAGCTAAGTTGGTATAGCAAGACCATGCCAAGGTAAGGGTCGCGAGTACGAGTCACGTTTcactct                                                               >  2:1759880/1‑145 (MQ=11)
                                                                                        ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTTCTTTAAGGCACAGACCACCCAAGCGGGAATAGCTAAGTTGGTATAGCACGACCATGCCAAGTTCGGGGTCGCGAGTTCGAGTCTACTTTCACGCt                                                               >  2:1665009/1‑150 (MQ=11)
                                                                                           aGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAATGAACACACCACCCAAGCGGGAATATCTAAGTTTTTATAGCAACACCTAGCCAAGTTCGGGGTCGCGACTACGACTCGCGTTTCCCGCTGCa                                                            >  2:1117742/1‑150 (MQ=11)
                                                                                                ggggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGGAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCAAGACCTTGACAAGGTCGGGGTCGCGAGTTCGAGTCGCGTTTCCCGATACAAAttt                                                       >  2:1701008/1‑149 (MQ=14)
                                                                                                             cGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCACAGTTCGTAGAGCACGACCTAGCCAAGGTCGGCGACGCGAGTTCGAGTCTCGTTTCCCGATACAAATTCTTCTATCAATaaa                                          >  1:1074371/1‑150 (MQ=17)
                                                                                                                    ggTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCca                                   <  1:635458/149‑1 (MQ=255)
                                                                                                                     gTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATCTCcac                                  >  1:481391/1‑150 (MQ=255)
                                                                                                                      tttCCCGCCCCAAATTTTAAAATTGCTGAAACTCACATACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                 <  1:595344/150‑1 (MQ=17)
                                                                                                                         cccGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                              >  2:1164208/1‑150 (MQ=255)
                                                                                                                               ccAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGCCGCGa                        >  2:1242101/1‑150 (MQ=255)
                                                                                                                                cAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACCGCGACGCGAt                       >  2:768797/1‑150 (MQ=255)
                                                                                                                                cAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAt                       >  1:1659878/1‑150 (MQ=255)
                                                                                                                                 aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGATTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCACGCTCCAAATTCTTATCTCAATAAAATATCCACAGCGACGCGATg                      >  1:1755175/1‑150 (MQ=255)
                                                                                                                                 aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTGTCAATAAAATATCCACAGCGACGCGATg                      >  1:1799798/1‑150 (MQ=255)
                                                                                                                                       tGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                >  2:1799355/1‑150 (MQ=255)
                                                                                                                                       tGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAGGCGTTAt                >  2:1229748/1‑150 (MQ=255)
                                                                                                                                       tGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGGTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                >  2:1223444/1‑150 (MQ=255)
                                                                                                                                        gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAtt               >  1:851015/1‑150 (MQ=255)
                                                                                                                                        gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAtt               >  1:418067/1‑150 (MQ=255)
                                                                                                                                         aaaGTGCTGTAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg              <  1:364064/150‑1 (MQ=255)
                                                                                                                                            gTGCTGTAAGGCACAGACCACCCAAGCAGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg           >  1:850073/1‑150 (MQ=255)
                                                                                                                                               cTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt        >  2:1421458/1‑150 (MQ=255)
                                                                                                                                                 gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtttt      >  2:72726/1‑150 (MQ=255)
                                                                                                                                                     ggCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttgtt  >  1:1213459/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTT  >  NC_000913/4392418‑4392716

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