Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 81,515 C→T Q148Q (CAG→CAA leuB ← 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091381,5150CT98.0% 154.7 / ‑4.2 49Q148Q (CAG→CAAleuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/24);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

TCGATCGACGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGA  >  NC_000913/81378‑81645
                                                                                                                                         |                                                                                                                                  
tcgatcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTg                                                                                                                                  >  2:752031/1‑140 (MQ=255)
 cgatcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGa                                                                                                                                 <  2:599566/140‑1 (MQ=255)
    tcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGGTTTTTCATGTTGACCGCTACCTTCGCGGCCTTTTGGTTGACCg                                                                                                                              >  1:628906/1‑140 (MQ=255)
    tcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTTTCCGCTACCTTCGCGGCCTTTTGGTTGACCg                                                                                                                              >  1:1007252/1‑140 (MQ=255)
    tcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCg                                                                                                                              >  1:893453/1‑140 (MQ=255)
    tcgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTAATATTGTTCGCGACGTTCGCGGCCTTTTGGTTGACCg                                                                                                                              >  1:302259/1‑140 (MQ=255)
     cgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGa                                                                                                                             <  1:7786/140‑1 (MQ=255)
            aCTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTGTTCATATTGTCCCCTAGCTTCGCGGACTTTAGGTTGACCGAAATAgct                                                                                                                      >  2:349600/1‑138 (MQ=255)
                       aCGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCa                                                                                                           >  2:216154/1‑140 (MQ=255)
                        cGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAg                                                                                                          <  1:21774/140‑1 (MQ=255)
                            tGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTcg                                                                                                      <  1:216154/140‑1 (MQ=255)
                              cGAGCAGATTCAAACTCGATGCTGGTGATCCGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTcgcg                                                                                                    <  1:316857/140‑1 (MQ=255)
                              cGAGCAGATTCAAACGCGATGCGGGCGATACTTTCTATCACAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTcgcg                                                                                                    <  2:412828/140‑1 (MQ=255)
                                  cAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGcaca                                                                                                <  1:672429/140‑1 (MQ=255)
                                  cAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGcaca                                                                                                >  1:885757/1‑140 (MQ=255)
                                    gATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGcacaca                                                                                              <  2:851438/140‑1 (MQ=255)
                                       tCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAg                                                                                           >  1:246536/1‑140 (MQ=255)
                                               gATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGa                                                                                   >  2:782146/1‑140 (MQ=255)
                                                 tGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAg                                                                                 <  1:596233/140‑1 (MQ=255)
                                                          tACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCg                                                                        >  1:117916/1‑140 (MQ=255)
                                                             gTTCGATCTCAAAACGGTGATATACCTCGGTTTCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCa                                                                     <  2:441622/140‑1 (MQ=255)
                                                              ttCGATCTCAAAACGGTTATACTCCTCGGTATCAACGGCCTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCaa                                                                    <  2:430128/140‑1 (MQ=255)
                                                                     tcAAAACGGTGATACACTTCGGTATCAAAGGTTTTGTCAAATTGTCCGTTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGc                                                             <  2:346906/140‑1 (MQ=255)
                                                                          aCGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCa                                                        <  2:793806/140‑1 (MQ=255)
                                                                          aCGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCa                                                        <  2:893453/140‑1 (MQ=255)
                                                                           cGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAg                                                       <  1:835210/140‑1 (MQ=255)
                                                                              tGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCgg                                                    >  1:849931/1‑140 (MQ=255)
                                                                              tGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAACCGTTTTCGGCAAAGTCTGCACGCAGCgg                                                    >  1:444337/1‑140 (MQ=255)
                                                                                atgccccTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACCCAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGAc                                                  <  1:744693/135‑1 (MQ=255)
                                                                                    acCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGaa                                              <  2:516914/140‑1 (MQ=255)
                                                                                        cGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGTCAACGTCTGCACGCAGCGGACAGAATGCt                                          >  1:383652/1‑140 (MQ=255)
                                                                                            aTCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCa                                      <  2:885757/140‑1 (MQ=255)
                                                                                              cAAAGGCTTTTTCATATTGTCCTCTACCTTCGCGGCCCTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGc                                    <  1:275355/140‑1 (MQ=255)
                                                                                              cAAAGGCTTTTTCATATAGTCCTCTACCTTCACGGACTTCTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGc                                    <  2:417631/140‑1 (MQ=255)
                                                                                                aaGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGTCGTTTGCGGCAATGTCTGCACGCAGCGGa                                                   <  1:241854/123‑1 (MQ=255)
                                                                                                aaGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGccc                                  <  2:117916/140‑1 (MQ=255)
                                                                                                aaGGCTTTTTCATATTGTCCGCTACCTTCGCGCCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGTCGTTTGCGGCAATGTCTGCACGCAGCGGa                                                   >  2:241854/1‑123 (MQ=255)
                                                                                                 aGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGcccc                                 <  1:958616/140‑1 (MQ=255)
                                                                                                     tttttCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAATCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGAt                             >  1:221870/1‑140 (MQ=255)
                                                                                                     tttttCATATTGTCCGCTACCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGAt                             >  1:11529/1‑140 (MQ=255)
                                                                                                       tttCATATTGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTGTGCACGCAGCGGACAGAATGCGTCCAGCCCCTGATAc                           >  1:417631/1‑140 (MQ=255)
                                                                                                               tGTCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGc                   >  2:229412/1‑140 (MQ=255)
                                                                                                                 tCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCg                 >  1:574343/1‑140 (MQ=255)
                                                                                                                 tCCGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCg                 >  2:916672/1‑140 (MQ=255)
                                                                                                                   cGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCggg               >  2:535499/1‑140 (MQ=255)
                                                                                                                   cGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCggg               >  2:827431/1‑140 (MQ=255)
                                                                                                                   cGCTACCTTCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCggg               >  1:191524/1‑140 (MQ=255)
                                                                                                                           tCGCGGCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGtt       >  1:501798/1‑140 (MQ=255)
                                                                                                                                gCCTTTTGGTTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGa  <  1:916672/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
TCGATCGACGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGA  >  NC_000913/81378‑81645

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