Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,530	2 bp→TT	coding (114‑115/1260 nt)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,530	0	C	T	100.0%	136.8 / NA	37	I38I (ATC→ATT)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/14);  total (23/14)
*	NC_000913	3,966,531	0	C	T	100.0%	129.2 / NA	37	L39L (CTG→TTG)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/14);  total (23/14)

CGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGA  >  NC_000913/3966289‑3966769                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                               cGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagcag                                                                                                                                                                                                                                        <  2:150076/251‑1 (MQ=255) cGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagc                                                                                                                                                                                                                                          >  1:110261/1‑249 (MQ=255)  gTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGGATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagccgc                                                                                                                                                                                                                                       >  2:13490/1‑248 (MQ=255)    aCTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAc                                                                                                                                                                                                                                     >  1:2373/1‑251 (MQ=255)                 cAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGAATGTATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGct                                                                                                                                                                                                                        <  2:303179/251‑2 (MQ=255)                      aGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGa                                                                                                                                                                                                                       >  1:385857/1‑247 (MQ=255)                                                               gatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAtt                                                                                                                                                                          >  1:135291/1‑251 (MQ=255)                                                                        aTACCATTCACAACCTTACGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGCTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcct                                                                                                                                                                 >  2:59627/1‑251 (MQ=255)                                                                        aTACCATTCACAACATTAAGTTCGAGATTTCCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcct                                                                                                                                                                 >  1:559666/1‑251 (MQ=255)                                                                                cacaACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTGTTTTTGCCATTTTGAAGCCGCCCGCCAAGAGTGGCGAGGATATCTTTGGTCATAGCGTACTGGCAATATTGCAGGTTAGTTTTGGCTTCCTCCGTTCCg                                                                                                                                                         >  2:254180/1‑251 (MQ=255)                                                                                  caacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGc                                                                                                                                                        >  1:226359/1‑250 (MQ=255)                                                                                               gagaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                                                          >  1:212866/1‑251 (MQ=255)                                                                                                agaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACCGCTCCTACCt                                                                                                                                         >  2:428314/1‑251 (MQ=255)                                                                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCATACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCg                                                                                                                                    >  1:140962/1‑251 (MQ=255)                                                                                                      aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTCTTTTTGCCATTTTGAAGCAGCACGCACAGAGCGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCTGCAGAAAGCTCCTACCTCGCCgg                                                                                                                                   >  1:294321/1‑251 (MQ=255)                                                                                                                tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgattac                                                                                                                         >  2:314002/1‑248 (MQ=255)                                                                                                                 aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                         >  1:12310/1‑250 (MQ=255)                                                                                                                                 aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                                                                          >  2:524601/1‑249 (MQ=255)                                                                                                                                         tACCGAATTAAAGAATTCGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGAAATTATGTTTGCCATTTTGAAGCAGCACGGAAAGAGTGGCGAAGATATCTTTGGGGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCTCCGGTCGTGATGTCATCTTCTTTTCCCCTTGCCaaa                                                                                                   <  2:344601/248‑1 (MQ=255)                                                                                                                                                         aCGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                                                                                >  2:380626/1‑251 (MQ=255)                                                                                                                                                           gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                                                                               >  1:52167/1‑250 (MQ=255)                                                                                                                                                            ccGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCAc                                                                             <  1:67428/251‑1 (MQ=255)                                                                                                                                                               gTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                                                                          <  1:251302/251‑1 (MQ=255)                                                                                                                                                                         tGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                >  1:132956/1‑251 (MQ=255)                                                                                                                                                                          gaacACTCCCGGCGAAAAGATGGGGCTGGAAAACCTGG‑CGCGTAGGCGTAAGCAGGACAGTATTTTGGCCATTTTGAAGCAGCTCGCAAAGATTGGGGAAGATATCTTTGTTGATGCCTTACTGGATATATTGCAGGATGGATTGGGTTTCCTCCGTTCCTCAAACTGCTCCTACCGCGCCGGTCCTGATGACATTTACGTTTCCCCGAGCTAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCAt                                                                 <  1:558633/246‑1 (MQ=255)                                                                                                                                                                             cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATctct                                                             <  1:101930/250‑1 (MQ=255)                                                                                                                                                                                 ctcGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt                                                          <  1:278418/249‑1 (MQ=255)                                                                                                                                                                                          aaTATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                               <  1:518464/251‑1 (MQ=255)                                                                                                                                                                                               gggggcggAAAACCCGGGCTCGTAGGCGTAAGCAGGAAATTATTTTGGCCTTTTTGAAGCAGCACGAAAAGAGGGGCGAAGATATCTTTGGGGTGGGCGTACTGGAGATATGGCAGGAGGGATTGGGTTTCCTCCGTCCCGCAGACAGCCCCTCCCTCGCGGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTGCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                           <  1:428314/245‑1 (MQ=255)                                                                                                                                                                                                 ggCTGGAAAACCGGG‑CTCGTATGGGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATTTCTTTGGTGAGGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCCCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTTCGTTTCCCCTAGCCAAATCCGCCTTTTCAACCTCCGCGCTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                        <  1:365157/251‑1 (MQ=255)                                                                                                                                                                                                     gaaaaaCCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGAAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTCCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                    <  1:402886/249‑1 (MQ=255)                                                                                                                                                                                                                      aTGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                     <  1:161679/250‑1 (MQ=255)                                                                                                                                                                                                                        gCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                  >  2:413706/1‑251 (MQ=255)                                                                                                                                                                                                                          gTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                >  2:37952/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaGCGGGACATTATTTTTGCCATTTTGAAGCAGACCGAAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACGGGAGATATTGCAGGAGGGATTTGGTTTCCTCCGTTCCGCAAACAGCTCTTCCCTCGCCGTTCCTGATGACTTCTACGTTCCCCCTAGCAAATTCCGCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc               <  1:508891/250‑1 (MQ=255)                                                                                                                                                                                                                               cAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa            >  1:310021/1‑250 (MQ=255)                                                                                                                                                                                                                                         tttttGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa  >  2:3242/1‑250 (MQ=255)                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                               CGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGG‑CTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGA  >  NC_000913/3966289‑3966769

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