Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,722,767 T→C L14P (CTT→CCT)  pssA → phosphatidylserine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,7670TC100.0% 120.6 / NA 35L14P (CTT→CCT) pssAphosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/15);  total (20/15)

GGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGC  >  NC_000913/2722692‑2722844
                                                                           |                                                                             
ggACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGcc                                                                           >  2:486885/1‑80 (MQ=255)
 gACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGccc                                                                          >  1:486894/1‑80 (MQ=255)
   cccTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAAc                                                                       <  1:486848/81‑1 (MQ=255)
          aTGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTAccc                                                                  >  2:486888/1‑79 (MQ=255)
            gAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTAccc                                                                  >  2:486892/1‑77 (MQ=255)
            gAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTAccc                                                                  >  2:486884/1‑77 (MQ=255)
             aaCAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGa                                                              <  1:486849/80‑1 (MQ=255)
                 gagaAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTtct                                                         >  1:486879/1‑81 (MQ=255)
                  agaagaAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTtctc                                                        <  2:486850/81‑1 (MQ=255)
                    aagaaATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCaa                                                      >  1:486890/1‑81 (MQ=255)
                       aaaTGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCa                                                   <  1:486851/81‑1 (MQ=255)
                         aTGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGt                                                 <  1:486852/81‑1 (MQ=255)
                         aTGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGt                                                 >  2:486891/1‑81 (MQ=255)
                            cACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTtga                                              >  2:486883/1‑81 (MQ=255)
                              cTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTtgatg                                            <  2:486853/81‑1 (MQ=255)
                                gtgATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTgatgat                                          >  2:486903/1‑81 (MQ=255)
                                   aTGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGt                                        <  1:486854/80‑1 (MQ=255)
                                        gtCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGAttt                                  >  1:486893/1‑81 (MQ=255)
                                        gtCAAAATTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGAttt                                  >  2:486889/1‑81 (MQ=255)
                                           aaaaTTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCtt                               <  2:486856/81‑1 (MQ=255)
                                           aaaaTTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCtt                               >  2:486909/1‑81 (MQ=255)
                                           aaaaTTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCt                                <  2:486855/80‑1 (MQ=255)
                                            aaaTTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCttt                              <  2:486857/81‑1 (MQ=255)
                                              aTTTAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTa                            >  1:486898/1‑81 (MQ=255)
                                                ttAAGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACg                          >  2:486896/1‑81 (MQ=255)
                                                  aaGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCt                        >  2:486910/1‑81 (MQ=255)
                                                   aGCGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTc                       >  2:486882/1‑81 (MQ=255)
                                                     cGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTccc                     <  1:486859/81‑1 (MQ=255)
                                                     cGTAATAAACATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTccc                     <  1:486858/81‑1 (MQ=255)
                                                              cATCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTc            <  2:486860/81‑1 (MQ=255)
                                                               aTCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTcc           >  1:486887/1‑81 (MQ=255)
                                                                tCAACAACACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCg          <  1:486861/81‑1 (MQ=255)
                                                                  aacaacACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCggg        <  1:486862/81‑1 (MQ=255)
                                                                  aacaacACCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCgg         >  2:486886/1‑80 (MQ=255)
                                                                        acCCTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACgc  >  1:486962/1‑81 (MQ=255)
                                                                           |                                                                             
GGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGC  >  NC_000913/2722692‑2722844

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