Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,125,561	(C)5→6	coding (186/1224 nt)	wcaI ←	putative colanic biosynthesis glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,125,556	1	.	C	100.0%	213.4 / NA	59	A64G (GCC→GGC)	wcaI	putative colanic biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/36);  total (23/36)

AATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑C‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAT  >  NC_000913/2125417‑2125690                                                                                                                                             |                                                                                                                                        aataatGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGACAGCGCCACACCGTggcgg                                                                                                                                          <  1:1119034/140‑1 (MQ=255)   taatGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑c‑c                                                                                                                                      <  2:1301563/140‑1 (MQ=255)      tGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑c‑cccc                                                                                                                                   >  2:179181/1‑140 (MQ=255)         gATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑c‑ccccctc                                                                                                                                <  2:519922/140‑4 (MQ=255)          aTCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑c‑ccccctct                                                                                                                               >  1:813910/1‑136 (MQ=255)           tCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑c‑ccccctctt                                                                                                                              <  1:518794/140‑6 (MQ=255)              agCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTtctc                                                                                                                            >  1:378017/2‑140 (MQ=255)               gCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCg                                                                                                                           <  2:93491/140‑1 (MQ=255)                cTTCCAGCGACGTTGCGCCATCTGTGGTAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAGCGGGCGTGTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑CGCCCCGCTTCTCgg                                                                                                                          >  2:343463/1‑139 (MQ=255)                        gACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAgcg                                                                                                                  <  1:72207/140‑1 (MQ=255)                          cGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGcc                                                                                                                >  1:1107630/1‑140 (MQ=255)                          cGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGcc                                                                                                                <  2:523391/140‑1 (MQ=255)                            ttGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAg                                                                                                              >  1:142961/1‑140 (MQ=255)                                  cATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGaa                                                                                                        <  1:186759/140‑1 (MQ=255)                                       gCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACACCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGtt                                                                                                   >  2:45699/1‑140 (MQ=255)                                              gaaaCTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGccc                                                                                            <  1:722949/140‑1 (MQ=255)                                                aaCTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCAc                                                                                          <  1:913521/140‑1 (MQ=255)                                                   tgctgACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATGCAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTg                                                                                       <  1:343463/140‑1 (MQ=255)                                                   tgctgACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTg                                                                                       <  1:457949/140‑1 (MQ=255)                                                    gctgACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGc                                                                                      >  1:273452/1‑140 (MQ=255)                                                        aCGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCAtt                                                                                  <  2:366748/140‑1 (MQ=255)                                                        aCGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCAtt                                                                                  <  1:1348156/140‑1 (MQ=255)                                                           gCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGTC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCg                                                                               <  2:744865/140‑1 (MQ=255)                                                           gCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCg                                                                               <  2:334804/140‑1 (MQ=255)                                                                    cccAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAg                                                                      >  2:1355457/1‑140 (MQ=255)                                                                    cccAGATGCAACAGGCGTTTCAGGGTGCTCGACTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAg                                                                      >  2:1492370/1‑140 (MQ=255)                                                                     ccAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAgg                                                                     <  2:170699/140‑1 (MQ=255)                                                                     ccAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAgg                                                                     <  2:1030771/140‑1 (MQ=255)                                                                      cAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgc                                                                    <  1:45699/140‑1 (MQ=255)                                                                       aGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcg                                                                   <  2:1006458/140‑1 (MQ=255)                                                                          tGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcggtg                                                                <  1:670668/140‑1 (MQ=255)                                                                          tGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcggtg                                                                <  1:965480/140‑1 (MQ=255)                                                                          tGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcggtg                                                                <  2:142961/140‑1 (MQ=255)                                                                               cAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTa                                                           <  1:834304/140‑1 (MQ=255)                                                                                aGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTaa                                                          <  1:1464253/140‑1 (MQ=255)                                                                                 ggCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAAt                                                         <  1:457336/140‑1 (MQ=255)                                                                                   cGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGTGGTAATGa                                                       >  2:654710/1‑140 (MQ=255)                                                                                   cGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGa                                                       >  1:611695/1‑140 (MQ=255)                                                                                    gTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGAc                                                      >  1:1245184/1‑140 (MQ=255)                                                                                     tttCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGAcc                                                     <  1:262181/140‑1 (MQ=255)                                                                                     tttCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGAcc                                                     >  1:1293648/1‑140 (MQ=255)                                                                                     tttCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGAcc                                                     <  1:179181/140‑1 (MQ=255)                                                                                       tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCGCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCg                                                   <  1:11145/140‑1 (MQ=255)                                                                                       tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCg                                                   <  1:1017014/140‑1 (MQ=255)                                                                                       tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCg                                                   <  1:505450/140‑1 (MQ=255)                                                                                       tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCg                                                   <  2:1051612/140‑1 (MQ=255)                                                                                        cAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGc                                                  <  1:725606/140‑1 (MQ=255)                                                                                        cAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGc                                                  <  2:1148980/140‑1 (MQ=255)                                                                                        cAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGc                                                  <  1:328386/140‑1 (MQ=255)                                                                                             tGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTc                                             >  1:443116/1‑140 (MQ=255)                                                                                              gCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCa                                            >  1:1311120/1‑140 (MQ=255)                                                                                                tCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATg                                          >  2:51555/1‑140 (MQ=255)                                                                                                       ttttGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTtgt                                   >  1:952957/1‑140 (MQ=255)                                                                                                             cacaTACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGCCCCGCACCTCATGACCTTGTgccgcc                             <  2:1245184/140‑1 (MQ=255)                                                                                                                aTACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc                          <  2:490446/140‑1 (MQ=255)                                                                                                                aTACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGAGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc                          <  2:1107630/140‑1 (MQ=255)                                                                                                                 tACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagcc                         <  1:973155/140‑1 (MQ=255)                                                                                                                 tACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagcc                         >  1:285530/1‑140 (MQ=255)                                                                                                                 tACAGCGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagcc                         <  1:811045/140‑1 (MQ=255)                                                                                                                      cGGGCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTcc                    >  1:826102/1‑140 (MQ=255)                                                                                                                       gggCAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTcca                   >  1:1425880/1‑140 (MQ=255)                                                                                                                          cAGCGCCACACCGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTccacca                >  2:897482/1‑140 (MQ=255)                                                                                                                                    ccGTGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGgt      <  2:611695/140‑1 (MQ=255)                                                                                                                                       tGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTa   >  2:327596/1‑140 (MQ=255)                                                                                                                                       tGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTa   >  2:828265/1‑140 (MQ=255)                                                                                                                                       tGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTa   >  2:734675/1‑140 (MQ=255)                                                                                                                                       tGGCGGCC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTa   >  1:995966/1‑140 (MQ=255)                                                                                                                                        ggcggcC‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCCCCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAt  >  1:1135592/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                        AATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGG‑C‑CCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAT  >  NC_000913/2125417‑2125690

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