Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,202	A→G	N118S (AAC→AGC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,202	0	A	G	100.0%	264.9 / NA	75	N118S (AAC→AGC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/36);  total (39/36)

TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCA  >  NC_000913/3061073‑3061335                                                                                                                                  |                                                                                                                                      tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                              >  2:1217961/1‑139 (MQ=255)            gCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACg                                                                                                                   >  2:517252/1‑139 (MQ=255)              tGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTc                                                                                                                 >  1:1513102/1‑139 (MQ=255)               ggTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCg                                                                                                                >  2:1478022/1‑139 (MQ=255)                   tttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCaa                                                                                                            >  2:1353606/1‑139 (MQ=255)                        aCAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCgg                                                                                                       <  1:517252/139‑1 (MQ=255)                           gcagaaGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    >  1:574524/5‑139 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    <  1:1478022/139‑1 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    <  2:781646/139‑1 (MQ=255)                                     ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                          >  2:94980/1‑139 (MQ=255)                                      cAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGac                                                                                         >  1:635938/1‑139 (MQ=255)                                         cGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACAcc                                                                                      >  2:1346655/1‑139 (MQ=255)                                         cGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACAcc                                                                                      >  1:990109/1‑139 (MQ=255)                                            tttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   <  2:1513102/139‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGACTACGACTCTGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                                 >  1:774610/1‑139 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGACTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                                 <  2:774610/139‑1 (MQ=255)                                               ttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGaa                                                                                <  2:738686/139‑1 (MQ=255)                                                 aTCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGa                                                                              <  1:1353606/139‑1 (MQ=255)                                                   cGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATTGACACCGTg                                                                                   <  1:1363665/132‑1 (MQ=255)                                                   cGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATTGACACCGTg                                                                                   >  2:1363665/1‑132 (MQ=255)                                                          aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGt                                                                     <  2:356687/139‑1 (MQ=255)                                                              ggCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCATGGCGGGTGATGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGATAGTCCTg                                                                 <  1:1056492/139‑1 (MQ=255)                                                               gccgccGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   >  2:1088243/1‑120 (MQ=255)                                                               gccgccGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   <  1:1088243/120‑1 (MQ=255)                                                               gccgccGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGt                                                                >  2:16119/1‑139 (MQ=255)                                                               gccgccGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGt                                                                <  2:990109/139‑1 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                               <  1:1213226/139‑1 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                               >  1:916731/1‑139 (MQ=255)                                                                   ccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                                 >  2:721819/1‑102 (MQ=255)                                                                   ccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                                 <  1:721819/102‑1 (MQ=255)                                                                   ccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGa                                                            >  2:542973/1‑139 (MQ=255)                                                                    cgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   <  1:1596162/115‑1 (MQ=255)                                                                    cgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   >  2:1596162/1‑115 (MQ=255)                                                                    cgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                           >  1:586116/1‑139 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                  <  2:916731/139‑1 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                  <  1:394819/139‑1 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                  <  1:94980/139‑1 (MQ=255)                                                                                  tttCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTg                                                                 >  2:592088/1‑119 (MQ=255)                                                                                  tttCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTg                                                                 <  1:592088/119‑1 (MQ=255)                                                                                          gCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGt                                     <  1:1079077/139‑1 (MQ=255)                                                                                             tttGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCggg                                                                                                      >  2:350530/1‑71 (MQ=255)                                                                                             tttGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCggg                                                                                                      <  1:350530/71‑1 (MQ=255)                                                                                             tttGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCgg                                  >  2:1431874/1‑139 (MQ=255)                                                                                              ttGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGt                                 >  2:222861/1‑139 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                            <  1:16119/139‑1 (MQ=255)                                                                                                     tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATg                          >  1:933271/1‑139 (MQ=255)                                                                                                      gCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    <  2:37614/64‑1 (MQ=255)                                                                                                      gCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    >  1:37614/1‑64 (MQ=255)                                                                                                        cacccaccGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGAcc                       <  1:222861/139‑1 (MQ=255)                                                                                                          cccaccGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                               <  2:609562/97‑1 (MQ=255)                                                                                                          cccaccGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                               >  1:609562/1‑97 (MQ=255)                                                                                                           ccaccGTGGCTACGACTCCGATGGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATg                    >  1:743832/1‑139 (MQ=255)                                                                                                           ccaccGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATg                    >  1:125087/1‑139 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGt                                                                     >  2:633107/1‑86 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGt                                                                     <  1:633107/86‑1 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATg                >  2:1216654/1‑139 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc              <  1:566609/139‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc              >  1:659600/1‑139 (MQ=255)                                                                                                                  ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    >  2:145935/1‑52 (MQ=255)                                                                                                                  ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                    <  1:145935/52‑1 (MQ=255)                                                                                                                   gCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   <  2:634254/68‑1 (MQ=255)                                                                                                                   gCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                   >  1:634254/1‑68 (MQ=255)                                                                                                                       cGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg             <  1:488326/134‑1 (MQ=255)                                                                                                                       cGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg             >  2:488326/1‑134 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                        >  1:1039426/1‑57 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                        <  2:1039426/57‑1 (MQ=255)                                                                                                                           tCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAc    <  2:125087/139‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                                 >  1:1233872/1‑44 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                                 <  2:1233872/44‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTAt                                                                                             >  2:77077/1‑48 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTAt                                                                                             <  1:77077/48‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGt         >  2:307204/1‑132 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGt         <  1:307204/132‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCa  >  2:1605100/1‑139 (MQ=255)                                                                                                                               aTAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGa                                                       <  1:256970/84‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCA  >  NC_000913/3061073‑3061335

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