Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,306,399 C→T W119* (TGG→TAG)  mqo ← malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,306,3990CT96.3% 175.8 / ‑1.7 54W119* (TGG→TAG) mqomalate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (26/26);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACGTTATCCTCGCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCCAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTAGTTCAGTTC  >  NC_000913/2306306‑2306493
                                                                                             |                                                                                              
aCGTTATCCTCGCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACt                                                                                          >  2:526675/1‑100 (MQ=255)
       ccTCGCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGa                                                                                   >  1:913285/1‑100 (MQ=255)
        cTCGCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGaa                                                                                  >  1:344752/1‑100 (MQ=255)
           gCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATc                                                                               <  1:295273/100‑1 (MQ=255)
            ccccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCt                                                                              <  1:132907/100‑1 (MQ=255)
            ccccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCt                                                                              >  2:401496/1‑100 (MQ=255)
            ccccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCt                                                                              >  1:940662/1‑100 (MQ=255)
             cccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTg                                                                             <  2:923725/100‑1 (MQ=255)
             cccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTg                                                                             <  2:273741/100‑1 (MQ=255)
             cccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTg                                                                             >  1:290107/1‑100 (MQ=255)
              ccAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGa                                                                            >  2:69141/1‑100 (MQ=255)
                aGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGaaa                                                                          <  1:547377/100‑1 (MQ=255)
                   cAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGc                                                                       >  2:890950/1‑100 (MQ=255)
                     aaGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCtt                                                                     >  1:1186954/1‑100 (MQ=255)
                     aaGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCtt                                                                     >  1:526441/1‑100 (MQ=255)
                       gCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCg                                                                   >  2:127502/1‑100 (MQ=255)
                        cTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGt                                                                  >  1:81596/1‑100 (MQ=255)
                        cTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGt                                                                  <  1:642430/100‑1 (MQ=255)
                         tCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGtt                                                                 >  1:170969/1‑100 (MQ=255)
                         tCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGtt                                                                 >  1:246281/1‑100 (MQ=255)
                          cATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTa                                                                >  2:718735/1‑100 (MQ=255)
                             atGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATg                                                             <  2:445894/100‑1 (MQ=255)
                              tGCGGAACGGTATTGATAAATGAACGCGGCGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATgg                                                            <  2:787086/100‑1 (MQ=255)
                                cGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCg                                                          <  2:58183/100‑1 (MQ=255)
                                cGGAACGGTATTGATAAATGAACGAAGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTCGAAGTGGCGGGAAATCTGAAATGCTTCGTTAATgg                                                            >  2:80565/1‑98 (MQ=255)
                                   aaCGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACt                                                       <  1:69142/100‑1 (MQ=255)
                                   aaCGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACt                                                       <  1:436507/100‑1 (MQ=255)
                                         aTTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCtttt                                                 <  1:890957/100‑1 (MQ=255)
                                         aTTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCtttt                                                 >  2:314715/1‑100 (MQ=255)
                                         aTTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCCAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCtttt                                                 >  2:81238/1‑100 (MQ=255)
                                            gATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCa                                              <  1:498106/100‑1 (MQ=255)
                                            gATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCa                                              <  1:178433/100‑1 (MQ=255)
                                              tAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAAt                                            <  1:268410/100‑1 (MQ=255)
                                              tAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCCAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAAt                                            >  1:664189/1‑100 (MQ=255)
                                               aaaTGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATa                                           >  1:1048125/1‑100 (MQ=255)
                                                    aaCGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGCAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGAt                                      >  1:1290218/1‑100 (MQ=255)
                                                              aCGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCg                            <  1:750972/100‑1 (MQ=255)
                                                               cGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGTGAAATCTGAAATGCTTCGTGAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCgg                           >  2:709885/1‑100 (MQ=255)
                                                               cGCAGCACGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCATTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCgg                           >  2:184391/1‑100 (MQ=255)
                                                                gcagcaCGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGc                          <  1:785499/100‑1 (MQ=255)
                                                                  agcaCGCCGCGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGt                        <  2:253043/100‑1 (MQ=255)
                                                                        ccgcGCTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCg                  <  2:1141693/100‑1 (MQ=255)
                                                                          gcgcTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCggg                <  2:1179402/100‑1 (MQ=255)
                                                                          gcgcTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCggg                <  1:640074/100‑1 (MQ=255)
                                                                           cgcTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCgggg               >  1:547382/1‑100 (MQ=255)
                                                                             cTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGgtg             <  1:49319/100‑1 (MQ=255)
                                                                             cTCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGgtg             >  2:766128/1‑100 (MQ=255)
                                                                              tCAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGgtgt            <  1:1196479/100‑1 (MQ=255)
                                                                               cAACCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGgtgt            >  2:102638/1‑99 (MQ=255)
                                                                                 aCCTGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTagt         >  2:578250/1‑100 (MQ=255)
                                                                                    tGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTagttca      <  2:919135/100‑1 (MQ=255)
                                                                                    tGGTGCGCCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTagttca      <  2:940654/100‑1 (MQ=255)
                                                                                        gcgcCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTagttcagttc  <  2:913277/100‑1 (MQ=255)
                                                                                        gcgcCTAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTagttcagtt   >  2:834749/1‑99 (MQ=255)
                                                                                             |                                                                                              
ACGTTATCCTCGCCCCAGACAAAGCTCATATGCGGAACGGTATTGATAAATGAACGCGGAGTACGCAGCACGCCGCGCTCAACCTGGTGCGCCCAGAACTGGCGGGAAATCTGAAATGCTTCGTTAATGGCGACTGCTTTTTCAATACTGATGCTGCCATCGGCGTTTTGCGGGGTGTAGTTCAGTTC  >  NC_000913/2306306‑2306493

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: