Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,741	C→T	R109C (CGC→TGC)	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,741	0	C	T	98.3%	201.4 / ‑4.3	60	R109C (CGC→TGC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (31/27);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAC  >  NC_000913/3966523‑3967007                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                             ttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGGTTCGCCCGCCGACAGCAGGTGAATGCTATTTTTCGCTGCCTAACGTTAACGAAGTTAACTTTGCCAACACCGGAAAACCCCCCcacaaacaa                                                                                                                                                                                                          >  2:743643/1‑280 (MQ=255)                ntaccaaagaaggTGGCGAGGTTATCTTTGTGGATGGCGTACTGGGGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGca                                                                                                                                                                                                                 <  2:149815/252‑1 (MQ=255)                 gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATcc                                                                                                                                                                                                        >  1:149815/1‑272 (MQ=255)                  caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAACGAAGGTGAAGGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAACCGCCCTCAACACAATCCAC‑TTTTAGGACTCaa                                                                                                                                                                                        >  2:720026/1‑286 (MQ=255)                  caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAATTGGTTTCCTTCGTTCCGCCGACAGCTCCTACATCCTCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAACAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTAGTCAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAAc                                                                                                                                                                                       >  1:735833/1‑287 (MQ=255)                      cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCATAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGACAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAATTTAAccc                                                                                                                                                                                     >  1:236503/1‑285 (MQ=255)                               gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCaa                                                                                                                                                                         >  1:365606/1‑288 (MQ=255)                                     atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAACGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAATTTAACTTCGACAA‑ACCTTAAAAGCCTCGCCTCAAAATCCCCATTTTAAAATTTAAACCCCCTGctccccaacct                                                                                                                                                                     >  2:436779/1‑277 (MQ=255)                                                       tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTgc                                                                                                                                                                                                                                                                    >  2:427633/1‑175 (MQ=255)                                                       tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTgc                                                                                                                                                                                                                                                                    <  1:427633/175‑1 (MQ=255)                                                       tACTGGAGATATTGCAGGATGGATTTGGTTCCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATATCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTga                                                                                                                                                                                                <  1:106288/241‑1 (MQ=255)                                                          tgaagaTATTGCAGGATGGAGTTTGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGTAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACt                                                                                                                                                                                           <  1:565368/240‑1 (MQ=255)                                                                 atTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTc                                                                                                                                       >  1:56087/1‑288 (MQ=255)                                                                   tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGACCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTc                                                                                                                                       >  1:722508/1‑286 (MQ=255)                                                                              tttGGTTTCCTCCTTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAGAACGCCCGCAACATAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                                                                                                                                        <  1:548400/274‑1 (MQ=255)                                                                                                gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGAt                                                                                                        <  1:836570/288‑1 (MQ=255)                                                                                                       gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTGTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                   <  1:852097/254‑1 (MQ=255)                                                                                                          ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                                                    <  1:743643/282‑1 (MQ=255)                                                                                                              ccTCGCCTGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAATCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAt                                                                                          <  1:417507/288‑1 (MQ=255)                                                                                                                  gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                                                                                                                                                <  1:110640/136‑1 (MQ=255)                                                                                                                  gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                                                                                                                                                >  2:110640/1‑136 (MQ=255)                                                                                                                  gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGt                                                                                      <  1:436779/288‑1 (MQ=255)                                                                                                                  gCCGGTCCTGATGACATCTAAGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGCACGAAGGTGAATGCTATTTTgc                                                                                                                                                                                                                                                                    >  2:499564/1‑116 (MQ=255)                                                                                                                   ccgtgccTGATGACAGCTACGTTTCCCCTAGCCAAAGCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                     <  1:466511/283‑1 (MQ=255)                                                                                                                            atgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCCTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTgg                                                                                                          >  1:118531/1‑258 (MQ=255)                                                                                                                                 aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa                                                                                                                                                     >  1:86296/1‑210 (MQ=255)                                                                                                                                 aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa                                                                                                                                                     <  2:86296/210‑1 (MQ=255)                                                                                                                                 aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTACCGGTTCTACTGAAGATTTAACTGCTAGCGTACTGGATCTGGCATCACCTACCGGTCGTGGTCAGCGTg                                                                         >  1:404658/1‑286 (MQ=255)                                                                                                                                           ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                                                                                                                                                                                                    >  1:471815/1‑123 (MQ=255)                                                                                                                                           ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                                                                                                                                                                                                    <  2:471815/123‑1 (MQ=255)                                                                                                                                           ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTGCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTACCGGCCCTGGTCAGCGTGGTCTGGTtttg                                                              >  1:234904/1‑287 (MQ=255)                                                                                                                                                           ccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                                                       >  1:416102/1‑262 (MQ=255)                                                                                                                                                                         gCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGt                                                                                <  1:787449/239‑1 (MQ=255)                                                                                                                                                                                tGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTAGTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGAAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCACGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                        <  1:319420/192‑1 (MQ=255)                                                                                                                                                                                tGATAACATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTCACGAAGTTCACTTCGACAA‑ACCTGAAAACGCACGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGTGTATGGAACGTGGTAACGGTTCTAGTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                                                       <  1:120971/241‑1 (MQ=255)                                                                                                                                                                                 gATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                                                           <  1:220417/236‑1 (MQ=255)                                                                                                                                                                                  aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGGAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                      <  1:670359/288‑1 (MQ=255)                                                                                                                                                                                  aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                      >  1:823430/1‑288 (MQ=255)                                                                                                                                                                                       aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCg                                                                                        >  1:4699/1‑217 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                                                                                                                                                                                                    <  2:625929/77‑1 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                                                                                                                                                                                                    >  1:625929/1‑77 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTga                                                                                                                                                                                                >  1:506686/1‑111 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTga                                                                                                                                                                                                <  2:506686/111‑1 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTg                                                                                                                                                                                <  2:838034/127‑1 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTg                                                                                                                                                                                >  1:838034/1‑127 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGTAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa                                                  >  1:659181/1‑253 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                   <  1:336197/172‑1 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                   >  2:336197/1‑172 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCg                                                                                        >  1:701568/1‑215 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                                                                  >  1:188094/1‑237 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCa               >  1:284292/1‑288 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCCCa               >  1:749804/1‑288 (MQ=255)                                                                                                                                                                                         ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCAGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                 >  1:237317/1‑222 (MQ=255)                                                                                                                                                                                              gTAAGATTCGCCCGCCGAAAGAATGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa          <  1:720026/288‑1 (MQ=255)                                                                                                                                                                                              gTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAAc                       >  1:331601/1‑275 (MQ=255)                                                                                                                                                                                              gTAAGATTCGCCCGCCGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa          <  1:709061/288‑1 (MQ=255)                                                                                                                                                                                              gTAAGAGTCGCCCGCCGCAAGAAGGTGCATGCTATTTTGCGCTGCCGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTACAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGTTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCCGGCATCACCTATCGGTCGTGGTCAGCGGGGTCTGATTGTGGCACCGCCGAAAGGCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa          <  1:316806/288‑1 (MQ=255)                                                                                                                                                                                                      cgcccgccGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGCTTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAc  >  1:147578/1‑288 (MQ=255)                                                                                                                                                                                                        cccgccGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGt                                                                                      >  2:335954/1‑202 (MQ=255)                                                                                                                                                                                                        cccgccGAAAGAAGGTGAATGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGt                                                                                      <  1:335954/202‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                             TTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAA‑ACCTGAAAACGCCCGCAACAAAATCCTC‑TTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAC  >  NC_000913/3966523‑3967007

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