Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,818:1 (A)6→7 100% intergenic (+166/‑267) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,8121.A88.5% 177.1 / 17.8 61intergenic (+160/‑273)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base . (2/5);  new base A (24/30);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTC  >  NC_000913/781750‑781910
                                                               |                                                                                                  
gCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                                                                <  2:646479/100‑1 (MQ=16)
gCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                                                                >  2:216146/1‑100 (MQ=16)
 ccGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                                                               <  2:2674948/100‑1 (MQ=255)
   gCTGGCGGCGACGGCCCAACAGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                                                               <  1:921220/98‑1 (MQ=11)
     tCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATc                                                           >  2:763644/1‑100 (MQ=255)
      cgcgGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATcc                                                          >  1:1137575/1‑100 (MQ=255)
      cgcgGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCATAAAGGCGCTTTTTTGCTATTCAGGCATcc                                                          >  1:2523847/1‑100 (MQ=255)
        cggcgACGGCCCAAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc                                                        <  1:186761/100‑1 (MQ=255)
          gctaCGGCCCGAAGGGCGAGCGAAGCGAGTCCTCCTTCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa                                                      >  1:1232415/4‑99 (MQ=255)
            gACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAtt                                                    <  2:145713/100‑1 (MQ=255)
            caCGGCCCGAAGGGCGAGCGCAGCGAGTCATCCTGCACGCCCCACCAATGTAAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATTCAGGCATcc                                                          >  1:441242/2‑94 (MQ=14)
              cGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTc                                                  <  2:2883175/100‑1 (MQ=255)
                gCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                >  1:1342336/1‑100 (MQ=255)
                gCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                <  2:2603510/100‑1 (MQ=255)
                gCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACt                                               <  2:950582/100‑1 (MQ=255)
                 cccAAAGGGCGAGCGCAGCGAGTAATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                <  2:120237/99‑1 (MQ=17)
                  ccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt                                              >  2:596731/1‑100 (MQ=255)
                  ccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                             <  2:299092/100‑1 (MQ=255)
                   cGAAGGGCGAGCGAAGCGAGTCTTCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                             <  2:2517148/100‑1 (MQ=255)
                    gAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg                                            >  2:2574994/1‑100 (MQ=255)
                     aaGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGt                                           <  1:39786/100‑1 (MQ=255)
                        ggCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTaaa                                        >  1:466589/1‑100 (MQ=255)
                          cgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      >  2:2362006/1‑100 (MQ=255)
                          cgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      <  1:191501/100‑1 (MQ=255)
                          cgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      <  1:2303083/100‑1 (MQ=255)
                           gatcgaagctaGTAATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      <  1:731904/96‑1 (MQ=21)
                           gagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      <  1:309737/99‑1 (MQ=255)
                           gagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCt                                     >  2:436816/1‑100 (MQ=255)
                            agcggagcgaGTCATCCTGCACGTCCCACCCATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACGTTGTAAAcc                                      <  1:570551/98‑1 (MQ=18)
                            agcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtg                                    >  2:1900913/1‑100 (MQ=255)
                             gcgaagcgaGTCATCCTGCACGACCCAACAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTaaa                                        <  1:491692/95‑1 (MQ=255)
                              cgcagcgaCTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAcc                                      <  1:2253243/92‑1 (MQ=255)
                              cgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                  >  2:1017853/1‑100 (MQ=255)
                              cgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                  <  1:1468156/100‑1 (MQ=255)
                              cgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtga                                   <  2:2753044/99‑1 (MQ=255)
                               gaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgatg                                 >  1:1741958/1‑100 (MQ=255)
                               gaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgatg                                 <  2:1465408/100‑1 (MQ=255)
                                aagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgatga                                <  1:2302667/100‑1 (MQ=255)
                                  gcgaGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACa                              >  1:2464898/1‑100 (MQ=255)
                                     aGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATc                            <  2:1651778/99‑1 (MQ=255)
                                     aGTCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCg                           <  2:1344231/100‑1 (MQ=255)
                                       tCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTc                         <  1:1043815/100‑1 (MQ=255)
                                       tCATCCTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTc                         >  1:808277/1‑100 (MQ=255)
                                            cTGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGc                    <  2:2655752/100‑1 (MQ=255)
                                             tGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCt                   >  1:2434859/1‑100 (MQ=255)
                                             tGCACGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCt                   >  1:2872570/1‑100 (MQ=255)
                                               cACGGCCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTa                 <  1:2753091/100‑1 (MQ=255)
                                                aCGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTAc                <  2:2279099/100‑1 (MQ=255)
                                                 cGACCCACCAATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACt               >  1:599869/1‑100 (MQ=255)
                                                      caccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgca          >  2:463907/1‑100 (MQ=255)
                                                      caccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgca          >  2:1056770/1‑100 (MQ=255)
                                                      caccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgca          <  2:2106671/100‑1 (MQ=255)
                                                      caccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgca          >  2:1040790/1‑100 (MQ=255)
                                                        ccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag        >  1:514966/1‑100 (MQ=255)
                                                        ccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag        <  2:2218629/100‑1 (MQ=255)
                                                        ccaATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag        <  2:1087214/100‑1 (MQ=255)
                                                        ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagc       <  2:637490/100‑1 (MQ=255)
                                                         caATGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTTCTGTgcaagc       >  2:1692920/1‑100 (MQ=255)
                                                         aattGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagca      <  1:2144804/97‑1 (MQ=255)
                                                          aaTGTAAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagca      <  1:1086205/100‑1 (MQ=255)
                                                             gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAActc  >  1:1448820/1‑100 (MQ=255)
                                                               |                                                                                                  
GCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTC  >  NC_000913/781750‑781910

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