Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,676,056 C→T G102D (GGC→GAC)  pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,676,0560CT98.2% 182.4 / ‑5.3 54G102D (GGC→GAC) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/22);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CACACAGTTTGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCA  >  NC_000913/1675919‑1676186
                                                                                                                                         |                                                                                                                                  
cacacaGTTTGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCa                                                                                                                                  <  1:693111/140‑1 (MQ=255)
        ttGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCAc                                                                                                                          <  1:1404538/140‑1 (MQ=255)
         tGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACt                                                                                                                         <  1:845007/140‑1 (MQ=255)
                           aaCCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAg                                                                                                       <  2:949815/140‑1 (MQ=255)
                           aaCCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAg                                                                                                       <  1:2738/140‑1 (MQ=255)
                           aaCCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAg                                                                                                       <  1:1188528/140‑1 (MQ=255)
                                 tGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGc                                                                                                 <  2:500063/140‑1 (MQ=255)
                                  gAACGTTACCGCCCCGATGAAGATACCGAGGAACCCTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATGCCCGCGTCATGATGCAGATAGCTGTTAAAACCACCCAGCACTGCCGCCAGACCCACGAAGCTATGCa                                                                                                <  1:247007/140‑1 (MQ=255)
                                        tACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCg                                                                                          >  1:1446932/1‑140 (MQ=255)
                                                  aTGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTc                                                                                >  1:1351612/1‑140 (MQ=255)
                                                          aCCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCAttt                                                                        >  1:795696/1‑140 (MQ=255)
                                                               ggAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTc                                                                   >  2:1290992/1‑140 (MQ=255)
                                                                 aaCACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCAt                                                                 >  1:1261369/1‑140 (MQ=255)
                                                                 aaCACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCAt                                                                 >  2:834336/1‑140 (MQ=255)
                                                                  acacTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCAtt                                                                >  2:331705/1‑140 (MQ=255)
                                                                      ttCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCGCTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCaa                                                            <  1:1290992/140‑1 (MQ=255)
                                                                      ttCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCaa                                                            >  1:1263996/1‑140 (MQ=255)
                                                                      ttCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCaa                                                            <  2:1082298/140‑1 (MQ=255)
                                                                         cGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACtt                                                         >  1:565523/1‑140 (MQ=255)
                                                                           tCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTttc                                                       >  2:481489/1‑140 (MQ=255)
                                                                           tCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTttc                                                       >  2:82167/1‑140 (MQ=255)
                                                                            cAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTttct                                                      <  1:1264210/140‑1 (MQ=255)
                                                                              ggTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTttcttc                                                    >  2:912525/1‑140 (MQ=255)
                                                                                tGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGc                                                  >  2:1427141/1‑140 (MQ=255)
                                                                                tGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGc                                                  <  1:331705/140‑1 (MQ=255)
                                                                                  aaTATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCa                                                <  1:481489/140‑1 (MQ=255)
                                                                                  aaTATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCa                                                >  1:1223464/1‑140 (MQ=255)
                                                                                    tatTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGa                                              <  1:1066801/140‑1 (MQ=255)
                                                                                         aCCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGAt                                         >  2:717159/1‑140 (MQ=255)
                                                                                         aCCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGAt                                         >  2:683003/1‑140 (MQ=255)
                                                                                         aCCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGAt                                         >  2:599213/1‑140 (MQ=255)
                                                                                         aCCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGAt                                         >  2:1169282/1‑140 (MQ=255)
                                                                                         aCCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGAt                                         >  2:335237/1‑140 (MQ=255)
                                                                                              aaTCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCaa                                    >  2:1105765/1‑140 (MQ=255)
                                                                                              aaTCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCaa                                    >  1:1417361/1‑140 (MQ=255)
                                                                                                   gTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGcc                                                 >  2:35164/1‑122 (MQ=255)
                                                                                                   gTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGcc                                                 <  1:35164/122‑1 (MQ=255)
                                                                                                   gTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGcc                               <  1:1105765/140‑1 (MQ=255)
                                                                                                   gTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACAATTGccc                              <  2:1223464/140‑1 (MQ=255)
                                                                                                    tGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGccc                              <  2:1446932/140‑1 (MQ=255)
                                                                                                          tCCCGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCaa                        <  1:730161/140‑1 (MQ=255)
                                                                                                            ccGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGcc                               >  1:1083801/1‑131 (MQ=255)
                                                                                                            ccGCGTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGcc                               <  2:1083801/131‑1 (MQ=255)
                                                                                                              gcgTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGAc                    >  2:1186291/1‑140 (MQ=255)
                                                                                                                gTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGccc                              >  1:403614/1‑128 (MQ=255)
                                                                                                                gTCATGATGCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGccc                              <  2:403614/128‑1 (MQ=255)
                                                                                                                   atgatgCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGccc                              >  2:57645/1‑125 (MQ=255)
                                                                                                                   atgatgCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGccc                              <  1:57645/125‑1 (MQ=255)
                                                                                                                   atgatgCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCg               >  1:1497107/1‑140 (MQ=255)
                                                                                                                    tgatgCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGc              <  1:1186291/140‑1 (MQ=255)
                                                                                                                     gatgCAGATAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGcc             >  1:1075616/1‑140 (MQ=255)
                                                                                                                            aTAGCTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGa      >  1:1254938/1‑140 (MQ=255)
                                                                                                                                cTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCa  >  1:1153978/1‑140 (MQ=255)
                                                                                                                                cTGTTAAAGTCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCa  >  1:117958/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CACACAGTTTGCCGAACGCCACCACCGAACCCGTGAACGTTACCGCCCCGATGAAGATACCGAGGAACACTTCCGTCAGGTGAATATTGACCAGAATCGGTGCCATTCCCGCGTCATGATGCAGATAGCTGTTAAAGCCAACCAGCACTGCCGCCAGACCCACGAAGCTATGCAGGATCGCCACCAGTTCTGGCATTTCGGTCATTTCAACTTTCTTCGCCAGACGGATACCAATTGCCCCACCAATGACCATCGCCAGCAAGATCCA  >  NC_000913/1675919‑1676186

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