Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,320,763 T→C S149P (TCG→CCG)  proP → glycine betaine/L‑proline transporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,320,7630TC100.0% 148.3 / NA 43S149P (TCG→CCG) proPglycine betaine/L‑proline transporter ProP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/18);  total (25/18)

GCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGC  >  NZ_CP009273/4320659‑4320896
                                                                                                        |                                                                                                                                     
gCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              <  2:69087/114‑1 (MQ=255)
gCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              >  1:69087/1‑114 (MQ=255)
     aTACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAAt                                                                                                                     <  1:147510/118‑1 (MQ=255)
     aTACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAAt                                                                                                                     >  2:147510/1‑118 (MQ=255)
     aTACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTg                                                                                                <  1:484838/139‑1 (MQ=255)
     aTACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTg                                                                                                >  1:117476/1‑139 (MQ=255)
          gTCCTACGACACGCTTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTCAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTCAACGTGGCttt                                                                                           >  2:294458/1‑139 (MQ=255)
          gTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              <  2:31125/104‑1 (MQ=255)
          gTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              >  1:31125/1‑104 (MQ=255)
          gTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCg                                                                                                        <  2:458529/126‑1 (MQ=255)
          gTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCg                                                                                                        >  1:458529/1‑126 (MQ=255)
          gTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCtt                                                                                            >  1:372656/1‑138 (MQ=255)
           tCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCg                                                                                                                         <  1:96837/108‑1 (MQ=255)
           tCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCg                                                                                                                         >  2:96837/1‑108 (MQ=255)
           tCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCt                                                                                                             >  2:553716/1‑120 (MQ=255)
           tCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCt                                                                                                             <  1:553716/120‑1 (MQ=255)
              tACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATgg                                                                                       >  2:546678/1‑139 (MQ=255)
               aCGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATggg                                                                                      >  1:51446/1‑139 (MQ=255)
                     cGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAgctg                                                                                >  1:491203/1‑139 (MQ=255)
                     cGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAgctg                                                                                >  1:544752/1‑139 (MQ=255)
                              tttGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACtt                                                                       <  2:235038/139‑1 (MQ=255)
                                 gggCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              >  1:494655/1‑81 (MQ=255)
                                 gggCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATttg                                                                                                                              <  2:494655/81‑1 (MQ=255)
                                  ggCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGt                                                                   >  1:222099/1‑139 (MQ=255)
                                                          tAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCt                                               <  2:265509/135‑1 (MQ=255)
                                                            aGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGc                                         <  2:222099/139‑1 (MQ=255)
                                                               tGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGAc                                                                         <  2:137284/104‑1 (MQ=255)
                                                               tGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGAc                                                                         >  1:137284/1‑104 (MQ=255)
                                                                 gCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCTGAATACTCCCCTGACCGTAAACg                                                                                                  <  2:557065/77‑1 (MQ=255)
                                                                 gCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACg                                                                                                  >  1:557065/1‑77 (MQ=255)
                                                                 gCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCg                                    >  2:468759/1‑139 (MQ=255)
                                                                  cacaAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTg                                          <  2:62806/132‑1 (MQ=255)
                                                                  cacaAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTg                                          >  1:62806/1‑132 (MQ=255)
                                                                  cacaAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCgt                                   >  2:389805/1‑139 (MQ=255)
                                                                  cacaAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCgt                                   >  2:350467/1‑139 (MQ=255)
                                                                   acaAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCgtg                                  <  1:25314/139‑1 (MQ=255)
                                                                           tctcGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCt                                                                                                             <  2:265416/56‑1 (MQ=255)
                                                                           tctcGGTCGGCGGTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCt                                                                                                             >  1:265416/1‑56 (MQ=255)
                                                                                     cggTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCgtggtggtg                            <  2:415711/127‑1 (MQ=255)
                                                                                     cggTGAATATACCGGGGCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCgtggtggtg                            >  1:415711/1‑127 (MQ=255)
                                                                                                    ggCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAg                                                                                   >  2:374946/1‑57 (MQ=255)
                                                                                                    ggCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAg                                                                                   <  1:374946/57‑1 (MQ=255)
                                                                                                    ggCGCCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGgcgaagc  >  1:309708/1‑138 (MQ=255)
                                                                                                        tCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGTGCAGCTGGCTGGACTTCgg                                                                    <  2:564871/68‑1 (MQ=255)
                                                                                                        |                                                                                                                                     
GCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGC  >  NZ_CP009273/4320659‑4320896

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