Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,685,542	(C)6→5	coding (1047/1404 nt)	fumC ←	fumarase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,685,537	0	C	.	100.0%	222.9 / NA	54	coding (1052/1404 nt)	fumC	fumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (24/30);  total (24/30)

TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCA  >  NC_000913/1685398‑1685671                                                                                                                                            |                                                                                                                                       tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcgc                                                                                                                                        >  2:1221809/1‑140 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                       <  1:1210437/140‑1 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                       <  2:143664/140‑1 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                       <  2:830861/140‑1 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                       <  1:471963/140‑1 (MQ=255)   gCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                      <  1:846444/140‑1 (MQ=255)   gCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTCTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                      <  1:1322429/140‑1 (MQ=255)     tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGTAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                                    >  2:170303/1‑140 (MQ=255)     tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                                    >  1:331197/1‑140 (MQ=255)          aTTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                      <  2:594551/133‑1 (MQ=255)          aTTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                      >  1:594551/1‑133 (MQ=255)                ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGc                                                                                                                        <  2:985044/140‑1 (MQ=255)                ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGc                                                                                                                        <  2:651313/140‑1 (MQ=255)                 tCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTAGATCGcc                                                                                                                       >  1:922961/1‑140 (MQ=255)                     tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgt                                                                                                                   >  1:386419/1‑140 (MQ=255)                          aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTc                                                                                                              >  1:1076630/1‑140 (MQ=255)                            tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                            >  1:258901/1‑140 (MQ=255)                             gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGAt                                                                                                                           >  2:1224985/1‑124 (MQ=255)                             gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGAt                                                                                                                           <  1:1224985/124‑1 (MQ=255)                             gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                           >  2:70194/1‑140 (MQ=255)                             gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                           >  1:883967/1‑140 (MQ=255)                               gcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                         >  2:432922/1‑140 (MQ=255)                                cAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgt                                                                                                                   <  1:466887/129‑1 (MQ=255)                                cAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgt                                                                                                                   >  2:466887/1‑129 (MQ=255)                                  gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                                                                      >  2:696619/1‑140 (MQ=255)                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                       <  1:20288/136‑1 (MQ=255)                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                       >  2:20288/1‑136 (MQ=255)                                      ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                         <  2:639003/133‑1 (MQ=255)                                      ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                         >  1:639003/1‑133 (MQ=255)                                      ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                                  <  2:310062/140‑1 (MQ=255)                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                               <  2:1076630/140‑1 (MQ=255)                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                               <  2:806380/140‑1 (MQ=255)                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGAACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                                             <  1:935406/140‑1 (MQ=255)                                             aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                    >  2:1279448/1‑131 (MQ=255)                                             aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                    <  1:1279448/131‑1 (MQ=255)                                             aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACagaaga                                                                                           >  2:632855/1‑140 (MQ=255)                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                         <  1:1016280/140‑1 (MQ=255)                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                         <  2:331197/140‑1 (MQ=255)                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGTAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                        <  1:170303/140‑1 (MQ=255)                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                        <  2:386419/140‑1 (MQ=255)                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                        <  2:883967/140‑1 (MQ=255)                                                   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                     >  1:1022978/1‑140 (MQ=255)                                                      tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                                  <  1:1475366/140‑1 (MQ=255)                                                       gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                                 >  2:654536/1‑140 (MQ=255)                                                       gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTAGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                                 <  2:922961/140‑1 (MQ=255)                                                            cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                            >  1:617662/1‑140 (MQ=255)                                                             tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                           <  2:22681/140‑1 (MQ=255)                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        <  2:1022978/140‑1 (MQ=255)                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        <  2:617662/140‑1 (MQ=255)                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGAAGAGCATGGTTAATGCCTcaca                                                                        <  1:632855/140‑1 (MQ=255)                                                                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        >  2:513550/1‑139 (MQ=255)                                                                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        <  1:513550/139‑1 (MQ=255)                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        <  1:1257013/133‑1 (MQ=255)                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        >  2:1257013/1‑133 (MQ=255)                                                                                tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        >  2:499541/1‑124 (MQ=255)                                                                                tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                        <  1:499541/124‑1 (MQ=255)                                                                                        aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                                <  2:258901/140‑1 (MQ=255)                                                                                             ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                           <  2:1415652/140‑1 (MQ=255)                                                                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                          >  1:1300698/1‑140 (MQ=255)                                                                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                          >  2:298921/1‑140 (MQ=255)                                                                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                          >  2:998197/1‑140 (MQ=255)                                                                                                                   gttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                     <  1:654536/140‑1 (MQ=255)                                                                                                                        gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGa                >  2:1488683/1‑140 (MQ=255)                                                                                                                               gTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt         <  2:869198/140‑1 (MQ=255)                                                                                                                                      gAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCa  <  1:1221809/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                       TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCA  >  NC_000913/1685398‑1685671

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