Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,932 A→C 11.1% D325E (GAT→GAG pntAA ← pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,9320AC11.1% 111.3 / 8.0 45D325E (GAT→GAGpntAApyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base C (2/3);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTC  >  NC_002947/165783‑166070
                                                                                                                                                     |                                                                                                                                          
ccTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCa                                                                                                                                            >  1:1200497/1‑150 (MQ=255)
         tttttGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACtg                                                                                                                                   <  1:576643/150‑2 (MQ=255)
         attttGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTc                                                                                                                                   >  2:977419/2‑150 (MQ=255)
             tGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCt                                                                                                                               <  2:859917/150‑1 (MQ=255)
                       cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGcccgacc                                                                                                                     >  1:1057491/1‑145 (MQ=255)
                       cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                     >  2:1412162/1‑150 (MQ=255)
                           ccATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTg                                                                                                                 <  1:477530/150‑1 (MQ=255)
                             aTCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGt                                                                                                               >  2:209494/1‑150 (MQ=255)
                             aTCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGt                                                                                                               >  1:376215/1‑150 (MQ=255)
                             aTCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGt                                                                                                               >  1:413310/1‑150 (MQ=255)
                                       aTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCAc                                                                                                                                     >  1:741883/1‑118 (MQ=255)
                                       aTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCAc                                                                                                                                     <  2:741883/118‑1 (MQ=255)
                                            gcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatgat                                                                                                <  1:1412162/150‑1 (MQ=255)
                                               ggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAATGACAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAc                                                                                             <  1:1085738/150‑1 (MQ=255)
                                               ggcCGCGACGATGTCGTCTTCGAGGTTGATGATCACCACGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAc                                                                                             <  1:436051/150‑1 (MQ=255)
                                                     gacgaTGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGtt                                                                                       <  2:413310/150‑1 (MQ=255)
                                                            tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCTGACGATGATCACGCCGTTCTCCTGc                                                                                <  1:977419/150‑1 (MQ=255)
                                                                 ttCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCCGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGcaccac                                                                           >  1:1189096/1‑150 (MQ=255)
                                                                    gAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTg                                                                        >  2:1296575/1‑150 (MQ=255)
                                                                       gTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACt                                                                                                                                    >  1:1119/1‑87 (MQ=255)
                                                                       gTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACt                                                                                                                                    <  2:1119/87‑1 (MQ=255)
                                                                             gatCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGAc                                                               >  1:378897/1‑150 (MQ=255)
                                                                                      gcCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCgg                                                      >  1:203233/1‑150 (MQ=255)
                                                                                       ccATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgat                                                                                                   <  2:365832/104‑1 (MQ=255)
                                                                                       ccATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgat                                                                                                   >  1:365832/1‑104 (MQ=255)
                                                                                        cATCCTTGTCGAACCGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGt                                                            <  2:330338/142‑1 (MQ=255)
                                                                                        cATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGt                                                            >  1:330338/1‑142 (MQ=255)
                                                                                            cTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                >  1:451293/1‑150 (MQ=255)
                                                                                            cTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGCCTGCGACGGTCAGCGGGCAGtt                                                <  1:623563/150‑1 (MQ=255)
                                                                                              tGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCCACGGTCACCGGGCAGTTgc                                              >  2:285797/1‑150 (MQ=255)
                                                                                                      agcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTgg                                      >  1:1168045/1‑150 (MQ=255)
                                                                                                             tCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAg                               >  1:669678/1‑150 (MQ=255)
                                                                                                              cATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGc                              >  1:393055/1‑150 (MQ=255)
                                                                                                                tGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCa                            >  1:1387947/1‑150 (MQ=255)
                                                                                                                 gAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAg                           >  1:16616/1‑150 (MQ=255)
                                                                                                                  aaGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAgg                          >  1:1317628/1‑150 (MQ=255)
                                                                                                                       cagcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                     >  2:350155/1‑54 (MQ=255)
                                                                                                                       cagcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                     <  1:350155/54‑1 (MQ=255)
                                                                                                                         gcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCa                                                   <  2:462002/118‑1 (MQ=255)
                                                                                                                         gcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTGGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCa                                                   >  1:462002/1‑118 (MQ=255)
                                                                                                                           agGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGacc                 <  2:1200497/150‑1 (MQ=255)
                                                                                                                                  gtgcgtACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGc          >  1:1191585/1‑150 (MQ=255)
                                                                                                                                       tACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGc     <  2:669678/150‑1 (MQ=255)
                                                                                                                                          aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGa                                                                >  1:414084/1‑88 (MQ=255)
                                                                                                                                          aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGa                                                                <  2:414084/88‑1 (MQ=255)
                                                                                                                                          aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGATGCAGCGAGGTCGATGAGCACCGAGCCGGGCTTc  <  2:376215/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
CCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTC  >  NC_002947/165783‑166070

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: