Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,723,976 A→G H417R (CAC→CGC)  pssA → phosphatidylserine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,723,9760AG97.5% 124.9 / NA 41H417R (CAC→CGC) pssAphosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/22);  minor base C (1/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACCTGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGT  >  NC_000913/2723845‑2724106
                                                                                                                                   |                                                                                                                                  
taaCCTGAACCCGCGCTCCTGGTGTCTGGAACTTGGAAAGGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATcgcg                                                                                                                            <  2:320133/140‑1 (MQ=255)
taaCCTGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATcgcg                                                                                                                            <  2:529751/140‑1 (MQ=255)
taaCCTGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATcgcg                                                                                                                            <  2:378874/140‑1 (MQ=255)
taaCCTGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATcgcg                                                                                                                            <  2:84830/140‑1 (MQ=255)
     tGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTg                                                                                                                       <  1:410675/140‑1 (MQ=255)
       aaCCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCTTCGTTAAGCGCTATCGCGATCTGCa                                                                                                                     >  1:206202/1‑140 (MQ=255)
       aaCCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCa                                                                                                                     >  1:142078/1‑140 (MQ=255)
       aaCCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCCAGCATACAACCCTCGTTAAGCCCTATCGCGATCTGCa                                                                                                                     >  1:300160/1‑140 (MQ=255)
            gcgcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTa                                                                                                                >  1:656479/1‑140 (MQ=255)
            gcgcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTa                                                                                                                >  1:401259/1‑140 (MQ=255)
             cgcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTAt                                                                                                               >  2:523141/1‑140 (MQ=255)
              gcgcCTGGCGTCTGGATGTGGATAACGCCATGTCGGTCCGCGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTAtt                                                                                                              <  2:294822/140‑1 (MQ=255)
              gcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATGCGCGAGCAAACCCCCAGCGTTAAGCGCTATGGCGATCTGCAAAGTAtt                                                                                                              >  2:181244/1‑140 (MQ=255)
              gcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATAGCACCATCGTTAGGCTCTATCGGGATCTGCAAAGTAtt                                                                                                              >  1:199348/1‑140 (MQ=255)
              gcgcCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCCCGATCGTTAAGCGCTATCGCGATCTGCAAAGTAtt                                                                                                              >  2:296208/1‑140 (MQ=255)
                       gTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTAt                                                                                                     <  2:142078/140‑1 (MQ=255)
                             atctggaAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTg                                                                                               <  1:110230/140‑1 (MQ=255)
                                tggaAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAg                                                                                            <  1:523141/140‑1 (MQ=255)
                                                  tCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATccgc                                                                          >  1:161604/1‑140 (MQ=255)
                                                        aTCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTg                                                                    <  2:161604/140‑1 (MQ=255)
                                                                   gAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCg                                                         >  1:630965/1‑140 (MQ=255)
                                                                   gAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGAATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCg                                                         >  1:425462/1‑140 (MQ=255)
                                                                    aGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGc                                                        >  2:656241/1‑140 (MQ=255)
                                                                        ggCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCg                                                    <  2:497797/140‑1 (MQ=255)
                                                                         gcgcCACAGCGAGAGCAAGAACTGGTGCTTATCCGCGCGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGa                                                   <  1:181244/140‑1 (MQ=255)
                                                                          cgcCACAGCGAGCGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTACTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGAc                                                  >  2:362809/1‑140 (MQ=255)
                                                                          cgcCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCTTAAACTCTTCCACCCTTTGCGCCGTATCCGCATCGAc                                                  >  1:238030/1‑140 (MQ=255)
                                                                           gcCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGAcc                                                 >  2:504259/1‑140 (MQ=255)
                                                                                  ggagagaAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTTAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGattaat                                          <  2:628182/139‑1 (MQ=255)
                                                                                  cGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGattaat                                          <  2:51238/140‑1 (MQ=255)
                                                                                         aagaaCTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGc                                   <  1:362809/140‑1 (MQ=255)
                                                                                          agaaCTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCa                                  <  2:113148/140‑1 (MQ=255)
                                                                                           gaaCTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCAt                                 <  2:515043/140‑1 (MQ=255)
                                                                                               tGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTg                             <  2:630965/140‑1 (MQ=255)
                                                                                               tGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCAGTATCCGCATCGACCGATTAATTAGCCGCATCCTg                             <  2:425462/140‑1 (MQ=255)
                                                                                                ggAGCTGATCCGCGAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGt                            <  1:656241/140‑1 (MQ=255)
                                                                                                     tGATCCGCGAGCATACCTCCATGGTTAAGCGCTATGGCGATCTGCAAAATATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATca                       <  2:238030/140‑1 (MQ=255)
                                                                                                          cgcgAGCATACCACCATAGTTATGCGCTATCGCGACCTGTAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAAcc                  <  2:267070/140‑1 (MQ=255)
                                                                                                          cgcgAGCATACCACCACGGTTATGCGCTATCCCGATCTGGAAAGTATTCCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAAcc                  <  2:300160/140‑1 (MQ=255)
                                                                                                            cgAGCATACCACCATCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAAcccc                >  1:488694/1‑140 (MQ=255)
                                                                                                                         aTCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACgggg   >  2:485869/1‑140 (MQ=255)
                                                                                                                          tCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgt  <  2:92016/140‑1 (MQ=255)
                                                                                                                          tCGTTAAGCGCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGAACGATTAATTAGCGGGATTCTGTAATAACAACCCCGTCCTGTACGGGgt  >  2:317837/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                  
TAACCTGAACCCGCGCGCCTGGCGTCTGGATCTGGAAAACGCCATTTTGATCCACGATCCGCAACTTGAGCTGGCGCCACAGCGAGAGAAAGAACTGGAGCTGATCCGCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGT  >  NC_000913/2723845‑2724106

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