Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,613,359 C→T P317S (CCA→TCA)  yjjU → putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,613,3590CT91.9% 118.7 / 2.3 37P317S (CCA→TCA) yjjUputative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (19/15);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

CGAAGCCATTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGCCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGAC  >  NC_000913/4613213‑4613501
                                                                                                                                                  |                                                                                                                                              
cgAAGCCATTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAgt                                                                                                                                            >  1:354329/1‑151 (MQ=255)
cgAAGCCATTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAgt                                                                                                                                            >  1:600043/1‑151 (MQ=255)
  aaGCCATTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAgtgg                                                                                                                                          >  2:132482/1‑151 (MQ=255)
       aTTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACg                                                                                                                                     >  1:62776/1‑151 (MQ=255)
           cATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCgg                                                                                                                                 >  1:564696/1‑151 (MQ=255)
               gTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATc                                                                                                                             <  2:364739/151‑1 (MQ=255)
                    gcgcTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTg                                                                                                                        >  2:590627/1‑151 (MQ=255)
                          ggCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCAtt                                                                                                                   <  2:590271/150‑1 (MQ=255)
                           gCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGCCAGTGGTGACGCCGGAATCGATTGTCATTcc                                                                                                                 >  2:606695/1‑151 (MQ=255)
                            cAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCg                                                                                                                <  1:181988/151‑1 (MQ=255)
                               tCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCt                                                                                                             >  2:354738/1‑151 (MQ=255)
                                           gcTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGcc                                                                                                 <  2:583453/151‑1 (MQ=255)
                                              gCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCa                                                                                                <  2:606252/149‑1 (MQ=255)
                                               cGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACg                                                                                             >  1:403978/1‑151 (MQ=255)
                                               cGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACg                                                                                             <  2:185004/151‑1 (MQ=255)
                                                  gAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGCCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGAt                                                                                           >  2:250499/1‑150 (MQ=255)
                                                   aaGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATAc                                                                                         >  2:391372/1‑151 (MQ=255)
                                                     gACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGc                                                                                       >  1:381102/1‑151 (MQ=255)
                                                           aaaCTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGCCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTg                                                                                 >  1:252806/1‑151 (MQ=255)
                                                            aaCTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGc                                                                                <  2:578216/151‑1 (MQ=255)
                                                                   gCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAATCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGt                                                                          <  2:89293/150‑1 (MQ=255)
                                                                   gCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTg                                                                         >  1:359674/1‑151 (MQ=255)
                                                                    cGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTg                                                                         >  2:237378/1‑150 (MQ=255)
                                                                    cGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGa                                                                        <  2:134712/151‑1 (MQ=255)
                                                                       ttATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCg                                                                     >  1:10387/1‑151 (MQ=255)
                                                                                    aTTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAACAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGc                                                        >  2:240399/1‑151 (MQ=255)
                                                                                                        aaGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAAc                                     >  2:416292/1‑150 (MQ=255)
                                                                                                          gCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAAt                                  <  2:553669/151‑1 (MQ=255)
                                                                                                           cTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATg                                 >  2:617594/1‑151 (MQ=255)
                                                                                                               tAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAgg                              <  2:158864/150‑1 (MQ=255)
                                                                                                               tAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAgg                              >  1:84982/1‑150 (MQ=255)
                                                                                                                aaCTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGAt                            <  2:479742/151‑1 (MQ=255)
                                                                                                                 aCTGAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATc                           >  2:537352/1‑151 (MQ=255)
                                                                                                                    gAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCt                          <  1:537363/149‑1 (MQ=255)
                                                                                                                    gAAAAAGCGCCGCTTACCCGCCATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCt                          <  1:152322/149‑1 (MQ=255)
                                                                                                                                         cATCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGa   <  2:603148/151‑1 (MQ=255)
                                                                                                                                          aTCTGGTGTCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGac  >  2:244087/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGAAGCCATTACATAGTATCGCGCTTGGCAGTCGGATTCCGGCGCTGCGTGAAGACTATAAACTTGGGCGTTTATGCGGTCGTTATTTCCTCGCCACGGTTGGCAAGCTATTAACTGAAAAAGCGCCGCTTACCCGCCATCTGGTGCCAGTGGTGACGCCGGAATCGATTGTCATTCCGCCTGCGCCAGTCGCCAACGATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGAC  >  NC_000913/4613213‑4613501

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