Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,800,755 C→T V606I (GTT→ATT)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,800,7550CT98.0% 170.0 / ‑4.1 51V606I (GTT→ATT) thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/26);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

AGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTT  >  NC_000913/1800661‑1800849
                                                                                              |                                                                                              
aGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGc                                                                                           <  1:149826/100‑1 (MQ=255)
   cTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTg                                                                                        <  1:1273579/100‑1 (MQ=255)
    tGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGa                                                                                       >  2:691869/1‑100 (MQ=255)
     gCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGAt                                                                                      >  2:525278/1‑100 (MQ=255)
      cGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGAtt                                                                                     <  1:693260/100‑1 (MQ=255)
           gCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCAcc                                                                                <  1:1136261/100‑1 (MQ=255)
           gCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCAcc                                                                                <  2:159627/100‑1 (MQ=255)
               aTCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCtt                                                                            <  2:1255551/100‑1 (MQ=255)
                tctcTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCttt                                                                           >  2:217778/1‑100 (MQ=255)
                tctcTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCttt                                                                           <  1:449382/100‑1 (MQ=255)
                 ctctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTa                                                                          <  1:374203/100‑1 (MQ=255)
                 ctctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTa                                                                          <  1:374302/100‑1 (MQ=255)
                  tctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTAt                                                                         >  1:467320/1‑100 (MQ=255)
                   ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATc                                                                        >  2:445282/1‑100 (MQ=255)
                   ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATc                                                                        <  2:59916/100‑1 (MQ=255)
                     tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcac                                                                      >  1:414/1‑100 (MQ=255)
                     tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcac                                                                      <  2:1212190/100‑1 (MQ=255)
                          ctgCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAg                                                                 <  2:789345/98‑1 (MQ=255)
                          cAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAAGGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAg                                                                 <  2:698928/100‑1 (MQ=255)
                               tctcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAg                                                            <  2:358272/100‑1 (MQ=255)
                                ctcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGc                                                           >  2:1198287/1‑100 (MQ=255)
                                 tcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCa                                                          <  2:398480/100‑1 (MQ=255)
                                 tcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCa                                                          <  1:970359/100‑1 (MQ=255)
                                  cGATCACTTCGTTTACGTCCATGCGTCCCTGGGCTTTACCACGGCGGGTGCGAACGGCAATTTTGCATGATTCCACCTCTTTATCACCACGGTCCAGCat                                                         <  2:292183/100‑1 (MQ=255)
                                  cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                         >  2:698870/1‑100 (MQ=255)
                                  cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                         >  1:81754/1‑100 (MQ=255)
                                  cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                         >  2:427346/1‑100 (MQ=255)
                                   gATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCata                                                        >  1:1143664/1‑100 (MQ=255)
                                   gATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCata                                                        <  1:812790/100‑1 (MQ=255)
                                       aCTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATg                                                    >  1:818911/1‑100 (MQ=255)
                                        cTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATgg                                                   <  1:768403/100‑1 (MQ=255)
                                         ttCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATggg                                                  <  1:1031036/100‑1 (MQ=255)
                                            aTTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacg                                               <  2:414/100‑1 (MQ=255)
                                               tACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcga                                            >  2:400975/1‑100 (MQ=255)
                                                 cGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacg                                          <  1:691873/100‑1 (MQ=255)
                                                      aTGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAg                                     >  1:595309/1‑100 (MQ=255)
                                                      aTGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAg                                     <  1:1198298/100‑1 (MQ=255)
                                                              cAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcg                             >  1:96265/1‑100 (MQ=255)
                                                              cAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCAcg                             >  2:520162/1‑100 (MQ=255)
                                                                ggTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                           >  2:1062695/1‑100 (MQ=255)
                                                                ggTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                           <  1:1062705/100‑1 (MQ=255)
                                                                 gTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCgg                          <  2:598376/100‑1 (MQ=255)
                                                                    tttACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAtt                       >  2:59665/1‑100 (MQ=255)
                                                                      tACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAtttt                     >  1:255403/1‑100 (MQ=255)
                                                                        ccACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaa                   >  1:415113/1‑100 (MQ=255)
                                                                             gcggGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCa              >  1:461463/1‑100 (MQ=255)
                                                                              cggGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCaa             >  2:1129476/1‑100 (MQ=255)
                                                                                  tGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtctt         >  1:574575/1‑100 (MQ=255)
                                                                                      aaCGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCa     <  2:312722/100‑1 (MQ=255)
                                                                                        cGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCAtt   >  1:1191847/1‑100 (MQ=255)
                                                                                         ggCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCAttt  >  2:40899/1‑100 (MQ=255)
                                                                                              |                                                                                              
AGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTT  >  NC_000913/1800661‑1800849

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