Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,009,341 (TGCCT)7→5 intergenic (+189/‑63) unknown → / → unknown Threonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,009,3140C.100.0% 80.1 / NA 20intergenic (+162/‑99)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/8);  total (12/8)
*11,009,3150T.100.0% 80.1 / NA 20intergenic (+163/‑98)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/8);  total (12/8)
*11,009,3160T.100.0% 80.1 / NA 20intergenic (+164/‑97)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/8);  total (12/8)
*11,009,3170G.100.0% 80.1 / NA 20intergenic (+165/‑96)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/8);  total (12/8)
*11,009,3180C.100.0% 80.1 / NA 20intergenic (+166/‑95)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/8);  total (12/8)
*11,009,3190C.100.0% 80.1 / NA 20intergenic (+167/‑94)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/8);  total (12/8)
*11,009,3200T.100.0% 80.1 / NA 20intergenic (+168/‑93)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/8);  total (12/8)
*11,009,3210T.100.0% 80.1 / NA 20intergenic (+169/‑92)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/8);  total (12/8)
*11,009,3220G.100.0% 80.1 / NA 20intergenic (+170/‑91)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/8);  total (12/8)
*11,009,3230C.100.0% 80.1 / NA 20intergenic (+171/‑90)unknown/unknownThreonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/8);  total (12/8)

TTCACCGTGGAAGAATTTTTAGACATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATT  >  1/1009075‑1009534
                                                                                                                                                                                                                                               ||||||||||                                                                                                                                                                                                                   
ttCACCGTGGAAGAATTTTTAGACATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCATTAATAGGAGTAAAAAATATGcttgccttgcct                                                                                                                                                                                                                   <  1:483973/251‑1 (MQ=255)
                      aCATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgccttgcctaaacta                                                                                                                                                                                              <  1:109985/250‑7 (MQ=255)
                                                            ttattGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTaaaa                                                                                                                                             >  1:762537/1‑251 (MQ=255)
                                                                        gCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgccttgcctaaact                                                                                                                                                                                               <  1:563741/199‑6 (MQ=255)
                                                                        gCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCggag                                                                                                                                 >  1:651693/1‑251 (MQ=255)
                                                                        gCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCggag                                                                                                                                 >  1:515156/1‑251 (MQ=255)
                                                                        gCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCAAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgccttgcctaaac                                                                                                                                                                                                >  1:446700/1‑194 (MQ=255)
                                                                            aGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCAc                                                                                                                                                        >  1:230110/1‑224 (MQ=255)
                                                                                     aaTTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgcc                                                                                                                                                                                                                    >  1:416748/1‑165 (MQ=255)
                                                                                                  gAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACtttt                                                                                                        >  1:99769/1‑250 (MQ=255)
                                                                                                            tAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgccttgcctaaac                                                                                                                                                                                                >  1:82494/1‑158 (MQ=255)
                                                                                                                  gAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgccttgcct                                                                                                                                                                                                    >  1:125110/1‑152 (MQ=255)
                                                                                                                  gAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCAc                                                                                                                                                        <  1:509339/186‑1 (MQ=255)
                                                                                                                      catcatTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGc                                                                                   <  1:320664/251‑1 (MQ=255)
                                                                                                                             cAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAAct                                                                              >  1:330672/1‑249 (MQ=255)
                                                                                                                                    tttATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAAc                                                                     <  1:697985/251‑1 (MQ=255)
                                                                                                                                                               ttCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTacccaacttaaagctc                                          <  1:137326/251‑16 (MQ=255)
                                                                                                                                                                              cTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTacccaacttaaagc                                            <  1:275169/234‑14 (MQ=255)
                                                                                                                                                                                  atatTGTATAAAACAAAGTACACTTTGTACTTGATAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTCATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTacccaacttaaagc                                            >  1:34500/1‑217 (MQ=255)
                                                                                                                                                                                  atatTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATaaa                                                                                                                           >  1:42710/1‑151 (MQ=255)
                                                                                                                                                                                  atatTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacccaacttacccaacttaaagctcttttcccgaattgtttcga                       >  1:22655‑M1/1‑201 (MQ=255)
                                                                                                                                                                                       gTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGcttgccttgccttgccttgcct                                                                                                                                                                                                         <  1:79286/78‑1 (MQ=255)
                                                                                                                                                                                       gTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacccaacttacccaacttaaagctcttttcccgaattgtttcgatgtaa                  >  1:170108‑M1/1‑196 (MQ=255)
                                                                                                                                                                                       gTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTacccaacttaaagctcttttcc                                    <  1:734897/233‑22 (MQ=255)
                                                                                                                                                                                               cAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacccaacttacccaacttaaagctcttttcccgaattgtttcgatgtaagcggtcaa          >  1:511262‑M1/1‑188 (MQ=255)
                                                                                                                                                                                                   gTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTAcc                                                                >  1:71705/1‑193 (MQ=255)
                                                                                                                                                                                                       acTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacccaacttacccaacttaaagctcttttcccgaattgtttcgatgtaagcggtcaatttttatt  <  1:843693‑M1/251‑72 (MQ=255)
                                                                                                                                                                                                            gTACTTGACAACCAGTAATAGGAGTAAAAAATATG‑‑‑‑‑‑‑‑‑‑CTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTT‑AAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacccaacttacccaacttaaagctcttttcccgaattgtttcgatgtaagcggtcaattttta    <  1:636407‑M1/243‑70 (MQ=255)
                                                                                                                                                                                                                                               ||||||||||                                                                                                                                                                                                                   
TTCACCGTGGAAGAATTTTTAGACATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATT  >  1/1009075‑1009534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑