Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,724,058 C→T D444D (GAC→GAT pssA → phosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,724,0580CT97.2% 115.8 / NA 36D444D (GAC→GATpssAphosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/17);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAG  >  NC_000913/2723923‑2724192
                                                                                                                                       |                                                                                                                                      
acaGCGAGAAAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCAtcgatcga                                                                                                                                     <  1:794697/139‑1 (MQ=255)
 caGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATcgatcgat                                                                                                                                    >  2:837730/1‑139 (MQ=255)
                        gATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATca                                                                                                             >  1:1092356/1‑139 (MQ=255)
                         aTCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAGGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATGGATTAATTAGCCGCATCCTGTAATcac                                                                                                            >  2:713822/1‑139 (MQ=255)
                                  cATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTc                                                                                                   <  1:619983/139‑1 (MQ=255)
                                           aTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                          >  2:748667/1‑139 (MQ=255)
                                                 aaGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTc                                                                                                   >  1:518047/1‑124 (MQ=255)
                                                 aaGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTc                                                                                                   <  2:518047/124‑1 (MQ=255)
                                                    cACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATca                                                                                                             <  2:942883/111‑1 (MQ=255)
                                                    cACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATca                                                                                                             >  1:942883/1‑111 (MQ=255)
                                                       tATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTg                                                                              >  1:143797/1‑139 (MQ=255)
                                                         tCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGattaat                                                                                                                                <  2:1119236/87‑1 (MQ=255)
                                                         tCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGattaat                                                                                                                                >  1:1119236/1‑87 (MQ=255)
                                                                 tGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGt                                                                    >  1:775380/1‑139 (MQ=255)
                                                                         aTTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTa                                                            <  2:143797/139‑1 (MQ=255)
                                                                           tGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTAtt                                                          >  2:360541/1‑139 (MQ=255)
                                                                             ccGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCt                                                        >  1:225661/1‑139 (MQ=255)
                                                                                 ttATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                    <  2:775380/139‑1 (MQ=255)
                                                                                                 cGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGt                                    >  1:847822/1‑139 (MQ=255)
                                                                                                 cGTAAACTCATCCGCCGTTTGCGCCGTATACGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGCTTGTTTTTTGGAGGCCACGTGTTGCGTATTCTTTTTGTCTGTTCGCTATTGt                                    <  1:837730/139‑1 (MQ=255)
                                                                                                       cTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTgg                                                                             >  2:568914/1‑92 (MQ=255)
                                                                                                       cTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTgg                                                                             <  1:568914/92‑1 (MQ=255)
                                                                                                       cTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCg                                           <  2:951924/126‑1 (MQ=255)
                                                                                                       cTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCg                                           >  1:951924/1‑126 (MQ=255)
                                                                                                         cATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCa                                                                                                                        <  2:198335/47‑1 (MQ=255)
                                                                                                         cATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCa                                                                                                                        >  1:198335/1‑47 (MQ=255)
                                                                                                         cATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTa                                                            >  2:1081433/1‑107 (MQ=255)
                                                                                                         cATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTa                                                            <  1:1081433/107‑1 (MQ=255)
                                                                                                         cATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCt                            >  2:630793/1‑139 (MQ=255)
                                                                                                          aTCCGCCGTTTGCGCAGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGGGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTg                           >  2:1093413/1‑139 (MQ=255)
                                                                                                           tCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAAcccc                                                                                                      <  2:117546/63‑1 (MQ=255)
                                                                                                           tCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAAcccc                                                                                                      >  1:117546/1‑63 (MQ=255)
                                                                                                            ccgccgTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGTGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                         <  2:844613/139‑1 (MQ=255)
                                                                                                            ccgccgTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                         <  1:748667/139‑1 (MQ=255)
                                                                                                                cgTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCa                     <  2:1092356/139‑1 (MQ=255)
                                                                                                                cgTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCa                     <  1:360541/139‑1 (MQ=255)
                                                                                                                     gcgcCGTATCCGCATCGATCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                <  1:881609/139‑1 (MQ=255)
                                                                                                                            aTCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCAGCCTGTACGGGGTTTGTTGTTTGGAGGCCACGTTTTGCGTATTCTCTATGTCTGTTCGCTATTGTTACTTTCTGGCAGCAACAATATGGCTa         >  2:396492/1‑139 (MQ=255)
                                                                                                                                   tcgatcgaTTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAg  >  2:370990/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAG  >  NC_000913/2723923‑2724192

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