Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,028,170	G→A	R98* (CGA→TGA)	ECB_02828 ←	Polysialic acid transport protein kpsM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,028,170	0	G	A	100.0%	198.7 / NA	52	R98* (CGA→TGA)	ECB_02828	Polysialic acid transport protein kpsM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/25);  total (27/25)

TATAAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAA  >  REL606/3028098‑3028238                                                                         |                                                                     tataAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATa                                                                   <  1:798404/76‑1 (MQ=255)   taAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGCTCGTATCGATGGGTTTTACTGGTCaataat                                                                 <  1:1391481/76‑1 (MQ=255)   taAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCaataat                                                                 <  1:557891/76‑1 (MQ=255)   taAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCaataat                                                                 <  2:582562/76‑1 (MQ=255)   taAACACCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGTTCGTATCGATGGGTTTTACTGGTCaataat                                                                 <  1:1403264/76‑1 (MQ=255)     aaCAGCAACGTAAATCAGTGTCTCGAGCGGTNCGCGTGCACTGATCNTATCGNTGGGTTTTACTGGGCAATAATTa                                                               <  2:201188/76‑1 (MQ=255)       cAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTaaa                                                             <  1:610082/76‑1 (MQ=255)       cAGCAACGTAAATCAGTGTCGCNNGCAGTGCGCGNNCAATGATCGTATCGTTGGGTTTTACTGGTCAATAATTaaa                                                             <  1:1378528/76‑1 (MQ=255)        aGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaac                                                            >  1:153988/1‑76 (MQ=255)        aGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaac                                                            >  1:154010/1‑76 (MQ=255)          cAACGTAAGTCAGTTTCTCAAGCAGGGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaacaa                                                          <  2:315201/76‑1 (MQ=255)          cAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaacaa                                                          <  2:1083235/76‑1 (MQ=255)          cAACGTAAAGCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaacaa                                                          <  1:49913/76‑1 (MQ=255)           aaCGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAaacaac                                                         <  2:830983/76‑1 (MQ=255)               tAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCt                                                     >  2:839409/1‑76 (MQ=255)                aaaTCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTg                                                    >  2:578055/1‑76 (MQ=255)                aaaTCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTg                                                    >  2:1495775/1‑76 (MQ=255)                aaaTCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTCACTGGTCAATAATTAAACAACCCCTg                                                    <  1:305171/76‑1 (MQ=255)                  aTCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGAt                                                  <  1:465723/76‑1 (MQ=255)                   tCAGTGTCTTAAGCAGTGCGCGTNCAATGNTCGTAGCGATGGGTNTTCCTGCTCAATAATTAAACAACCCCTGAtt                                                 >  2:206174/1‑76 (MQ=255)                   tCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGAtt                                                 >  1:947028/1‑76 (MQ=255)                     aGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTCCTGGTCAATAATTAAACAACCCCTGATTCg                                               >  1:1152572/1‑76 (MQ=255)                        gtCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCtt                                            <  2:1361829/76‑1 (MQ=255)                         tctcAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATCCCCTTc                                           <  1:94688/76‑1 (MQ=255)                          ctcAAGCAGTGCGCGTGTAATGAGCGTATCGATGGGTTTTACTGGTGAATAATTAAACAACCCCTGATACGCTTCa                                          >  1:81955/1‑76 (MQ=255)                          ctcAAGCAGCGCGCGTGCAATGATCGTATCGATTGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCa                                          >  2:1310611/1‑76 (MQ=255)                            cAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAAt                                        >  1:373112/1‑76 (MQ=255)                            cAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAAt                                        <  2:121083/76‑1 (MQ=255)                            cAAGCAGTGCGCGTGCAATGATCGTATCGATGGGGTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAAt                                        >  1:373117/1‑76 (MQ=255)                             aaGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATa                                       >  2:1142848/1‑76 (MQ=255)                               gCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGc                                     <  1:438537/76‑1 (MQ=255)                                    gcgcGTGCAATGATNGTATCGATTGGTTTTACTGNTCAATCATNAAACAACCCCTGATTCGCTTCAATAGCACCTa                                <  2:1271630/76‑1 (MQ=255)                                    gcgcGTGCAATGATNGTATCGATGGGTTGTACTGNTCAATCATNAAACAACCCCTGATTCGCTCCAATAGCACCTa                                <  2:1271627/76‑1 (MQ=255)                                       cgTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAg                             <  2:1410093/76‑1 (MQ=255)                                        gTGCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGa                            <  1:1068786/76‑1 (MQ=255)                                          gCAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAAc                          >  2:700897/1‑76 (MQ=255)                                           cAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACg                         >  1:927257/1‑76 (MQ=255)                                           cAATGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACg                         >  2:1162908/1‑76 (MQ=255)                                             aTGATCGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGATTCAATAGCCCCTACAGAACGAt                       >  1:1315013/1‑76 (MQ=255)                                                  cGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTg                  >  1:551942/1‑76 (MQ=255)                                                  cGTATCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTg                  <  2:394317/76‑1 (MQ=255)                                                     aTCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATa               <  2:464272/76‑1 (MQ=255)                                                     aTCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATa               >  2:770074/1‑76 (MQ=255)                                                      tCGATGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATAc              <  2:633080/76‑1 (MQ=255)                                                        gATCGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTg            >  2:563490/1‑76 (MQ=255)                                                         aTGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGc           <  2:723862/76‑1 (MQ=255)                                                          tGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCt          <  1:991029/76‑1 (MQ=255)                                                          tGGGTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCt          >  1:1335166/1‑76 (MQ=255)                                                            ggTTTTACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGNTACTGCTaa        <  1:39564/76‑1 (MQ=255)                                                              ntttACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGNTGATACTGCTAAAg      >  2:217297/2‑76 (MQ=255)                                                               tttACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGa     >  1:889497/1‑76 (MQ=255)                                                                ttACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGAt    >  2:328019/1‑76 (MQ=255)                                                                ttACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGAt    >  1:737286/1‑76 (MQ=255)                                                                 tACTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATa   >  1:1474625/1‑76 (MQ=255)                                                                 tACTGGTCAAGAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATa   >  1:1417100/1‑76 (MQ=255)                                                                  aCTGGTCAATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATaa  <  2:843897/76‑1 (MQ=255)                                                                         |                                                                     TATAAACAGCAACGTAAATCAGTGTCTCAAGCAGTGCGCGTGCAATGATCGTATCGATGGGTTTTACTGGTCGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAA  >  REL606/3028098‑3028238

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