Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 783,000 C→A 8.7% Q306K (CAG→AAG)  nadA → quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913783,0000CA 8.7% 156.5 / 4.9 46Q306K (CAG→AAG) nadAquinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (26/16);  new base A (2/2);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAA  >  NC_000913/782856‑783145
                                                                                                                                                |                                                                                                                                                 
cTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAgg                                                                                                                                                >  3:106974/1‑148 (MQ=255)
  tATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGcc                                                                                                                                              <  3:215723/148‑1 (MQ=255)
  tATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGcc                                                                                                                                              >  3:181438/1‑148 (MQ=255)
    ttGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAtt                                                                                                                                           <  2:80544/149‑2 (MQ=255)
        ggCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAttcat                                                                                                                                        <  4:12226/148‑5 (MQ=255)
          cAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCaga                                                                                                                                      >  3:70672/1‑148 (MQ=255)
            aCCGATCGGGCTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAgg                                                                                                                                    >  3:350642/1‑148 (MQ=255)
                   ggggTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGaa                                                                                                                            <  2:407030/149‑1 (MQ=255)
                     ggTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACa                                                                                                                          <  1:279280/149‑1 (MQ=255)
                          tttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaa                                                                                                                      <  4:253370/148‑1 (MQ=255)
                           ttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaag                                                                                                                     >  3:313005/1‑148 (MQ=255)
                             cTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGACGCATTAGAACaggaagga                                                                                                                   <  3:247358/148‑1 (MQ=255)
                               aCAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGgaaggaag                                                                                                                 >  4:320570/1‑148 (MQ=255)
                               aCAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGc                                                                                                                <  1:292308/149‑1 (MQ=255)
                                 aaaaTGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCa                                                                                                               <  4:181438/148‑1 (MQ=255)
                                     tGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCAc                                                                                                          <  2:64562/149‑1 (MQ=255)
                                     tGCAGCAGGCGGTGCCAAAAAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCAc                                                                                                          <  2:405882/149‑1 (MQ=255)
                                       cagcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACg                                                                                                         <  4:106974/148‑1 (MQ=255)
                                        agcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                                       >  1:146063/1‑149 (MQ=255)
                                         gcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAgg                                                                                                      >  2:307296/1‑149 (MQ=255)
                                                   ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtga                                                                                            <  1:89069/149‑1 (MQ=255)
                                                   ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtga                                                                                            >  1:407206/1‑149 (MQ=255)
                                                       aTAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtgatga                                                                                         >  3:107910/1‑148 (MQ=255)
                                                        tAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCACATTGTCCGTGGATAGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTACATGTTGATGaa                                                                                        <  4:350642/148‑1 (MQ=255)
                                                         aaaGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGCCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGt                                                                                               <  3:383319/140‑1 (MQ=255)
                                                         aaaGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGCCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGt                                                                                               >  4:383319/1‑140 (MQ=255)
                                                                  cTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgaga                                                                             >  1:234106/1‑149 (MQ=255)
                                                                   tGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgaga                                                                             >  4:322048/1‑148 (MQ=255)
                                                                      aaGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAgg                                                                         <  1:167942/149‑1 (MQ=255)
                                                                           ccAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCt                                                                    <  1:307296/149‑1 (MQ=255)
                                                                            cAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTg                                                                   >  2:221657/1‑149 (MQ=255)
                                                                             aaCCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTg                                                                   >  3:189905/1‑148 (MQ=255)
                                                                              aCCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTgg                                                                  >  3:303165/1‑148 (MQ=255)
                                                                                 gCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGc                                                               >  3:313236/1‑148 (MQ=255)
                                                                                            tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGt                                                    >  3:145184/1‑148 (MQ=255)
                                                                                               aaCCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGc                                                >  2:279680/1‑149 (MQ=255)
                                                                                                       gCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGAtttt                                        >  1:40580/1‑149 (MQ=255)
                                                                                                        cAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTg                                       >  2:299740/1‑149 (MQ=255)
                                                                                                               gcgcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTac                                >  1:383653/1‑149 (MQ=255)
                                                                                                                cgcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTac                                >  3:61601/1‑148 (MQ=255)
                                                                                                                 gcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTaca                               <  4:342878/148‑1 (MQ=255)
                                                                                                                                  tctgcgaTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaata             >  2:141244/7‑149 (MQ=255)
                                                                                                                                   ctgcgaTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaata             >  4:210538/6‑148 (MQ=255)
                                                                                                                                    aTGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataa            >  4:172297/1‑148 (MQ=255)
                                                                                                                                            ccTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTa   >  1:180446/1‑149 (MQ=255)
                                                                                                                                              ttCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTaa  >  3:228100/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAA  >  NC_000913/782856‑783145

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