Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,009,349 C→T C60C (TGC→TGT pldB → lysophospholipase L2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,009,3490CT100.0% 194.8 / NA 50C60C (TGC→TGTpldBlysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/20);  total (30/20)

AAAACGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGCCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATC  >  NC_000913/4009204‑4009494
                                                                                                                                                 |                                                                                                                                                 
aaaaCGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCggg                                                                                                                                              >  1:830410/1‑151 (MQ=255)
  aaCGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGc                                                                                                                                            >  2:168149/1‑151 (MQ=255)
     gcgTTTGCTGCTTTTACCATGGGAACGCTGACTGATTTCTGGCGGCAGAGTGATGAAGCAGAGTTTACTGGGGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACAGGGTGGTAGTCATCTGTCCGGGGCg                                                                                                                                           >  2:244866/1‑149 (MQ=255)
         ttGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTg                                                                                                                                      >  2:637970/1‑150 (MQ=255)
          tgctgcTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTg                                                                                                                                      >  2:593389/1‑149 (MQ=255)
            ctgctTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTgag                                                                                                                                    >  2:228350/1‑149 (MQ=255)
               ctTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCt                                                                                                                                >  1:569275/1‑150 (MQ=255)
                ttttACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTAc                                                                                                                              >  2:689945/1‑151 (MQ=255)
                  ttACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGt                                                                                                                            >  1:798783/1‑151 (MQ=255)
                   tACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGt                                                                                                                            <  1:895786/150‑1 (MQ=255)
                          ggACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATg                                                                                                                    <  2:666763/151‑1 (MQ=255)
                           gACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGc                                                                                                                   >  2:892104/1‑151 (MQ=255)
                               gCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGaa                                                                                                               <  2:891631/151‑1 (MQ=255)
                                     tgaTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGcc                                                                                                         >  1:220018/1‑151 (MQ=255)
                                       aTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTa                                                                                                       <  2:474856/151‑1 (MQ=255)
                                         ttCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTa                                                                                                       <  2:538963/149‑1 (MQ=255)
                                         ttCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATg                                                                                                     >  2:766377/1‑151 (MQ=255)
                                          tCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATg                                                                                                     >  2:315682/1‑150 (MQ=255)
                                                gTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGtt                                                                                              >  2:710381/1‑151 (MQ=255)
                                                  cAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTcc                                                                                            >  1:127331/1‑151 (MQ=255)
                                                   aGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTcc                                                                                            >  2:676560/1‑150 (MQ=255)
                                                   aGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTcc                                                                                            >  2:653511/1‑150 (MQ=255)
                                                        gatgaAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTg                                                                                       <  1:10881/150‑1 (MQ=255)
                                                            aaGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGtt                                                                                  <  1:69060/151‑1 (MQ=255)
                                                              gCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTg                                                                                <  2:134500/151‑1 (MQ=255)
                                                               cAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGa                                                                               <  2:332081/151‑1 (MQ=255)
                                                                             gtgGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAatcatc                                                                  <  1:266510/150‑1 (MQ=255)
                                                                             gtgGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCg                                                                 <  2:335942/151‑1 (MQ=255)
                                                                                    aCATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATcg                                                          <  2:829734/151‑1 (MQ=255)
                                                                                         ccGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGc                                                     <  2:852745/151‑1 (MQ=255)
                                                                                          cGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGc                                                     >  2:80859/1‑150 (MQ=255)
                                                                                           ggTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGc                                                     >  2:624594/1‑149 (MQ=255)
                                                                                             tGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAggg                                                 >  2:599735/1‑151 (MQ=255)
                                                                                                gcTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACg                                              >  2:714653/1‑151 (MQ=255)
                                                                                                   ttGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGt                                             <  2:339030/149‑1 (MQ=255)
                                                                                                         gTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCg                                     >  1:75100/1‑151 (MQ=255)
                                                                                                         gTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCg                                     <  2:173249/151‑1 (MQ=255)
                                                                                                                cgcACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGtt                               >  2:146211/1‑150 (MQ=255)
                                                                                                                      gCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGAt                        <  1:168173/151‑1 (MQ=255)
                                                                                                                        aCCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATcc                      <  2:123405/151‑1 (MQ=255)
                                                                                                                          cATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGc                    >  1:790717/1‑151 (MQ=255)
                                                                                                                               ccGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCAtctc               <  1:426389/151‑1 (MQ=255)
                                                                                                                                   gtggtAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGc           >  2:555621/1‑151 (MQ=255)
                                                                                                                                   gtggtAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGc           >  1:858786/1‑151 (MQ=255)
                                                                                                                                     ggtAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCa          <  1:188482/150‑1 (MQ=255)
                                                                                                                                        aGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGtt      >  2:638768/1‑151 (MQ=255)
                                                                                                                                        aGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGtt      <  2:219992/151‑1 (MQ=255)
                                                                                                                                         gTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTAAGACATGTACCATTTGGGGTTTAATGTCTTAATAATCGACCATCGCGGGCTGGGACAGTCCGGTCGCATGTTAGCCGATCCGCATCTCGGGCATGtt      >  2:872193/1‑150 (MQ=255)
                                                                                                                                            aTCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATc  >  1:241114/1‑151 (MQ=255)
                                                                                                                                            aTCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATc  >  2:377904/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AAAACGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGCCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATC  >  NC_000913/4009204‑4009494

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