Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,850,545	C→T	A570V (GCA→GTA)	sppA →	protease IV, a signal peptide peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,850,545	0	C	T	97.3%	121.0 / NA	37	A570V (GCA→GTA)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/20);  minor base G (1/0);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCACGCTGGTTGCTGATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAT  >  NC_000913/1850300‑1850790                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      tATCACGCTGGTTGCTGATGAGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAATGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCTGGTAATGt                                                                                                                                                                                                                                                  >  2:282080/1‑251 (MQ=255)   tCACGCTGGTTGCTGATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTg                                                                                                                                                                                                                                                >  1:365928/1‑251 (MQ=255)          ggTTGCTGATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGa                                                                                                                                                                                                                                           <  1:305566/249‑1 (MQ=255)          ggTTGCTGATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGa                                                                                                                                                                                                                                           <  1:408345/249‑1 (MQ=255)                 gATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTc                                                                                                                                                                                                                                    <  2:340462/249‑1 (MQ=255)                           aTTCGACGGCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTa                                                                                                                                                                                                                        <  2:53955/251‑1 (MQ=255)                              cGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACTGGCTGGTCGATCGTCTCGGGGATTTCGATGATGCGGTCGCCAAGGCAGCAGAGCTGGCACAAGTGAAACAGTGGCATCTGGAGTATTACGTTGATGAAACGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGGAATGTGGCCAGATGCGTTCCAGGCCATGTTACCt                                                                                                                                                                                                                     >  2:328795/1‑251 (MQ=255)                                 cGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGAAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCa                                                                                                                                                                                                                  <  2:102097/251‑1 (MQ=255)                                         cAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGc                                                                                                                                                                                                          <  2:146856/251‑1 (MQ=255)                                            aTTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCTTCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGTCAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTTCACCGCTGGCCTc                                                                                                                                                                                                       <  1:311310/251‑1 (MQ=255)                                                            aGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGtt                                                                                                                                                                                       <  2:255984/251‑1 (MQ=255)                                                               gCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTaa                                                                                                                                                                                     >  1:385851/1‑250 (MQ=255)                                                                   cGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGCCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAaaagtg                                                                                                                                                                                >  2:159997/1‑251 (MQ=255)                                                                                       cAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCcgcg                                                                                                                                                            >  2:364605/1‑251 (MQ=255)                                                                                                 cGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGAcc                                                                                                                                                  >  2:92476/1‑251 (MQ=255)                                                                                                  gggCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGAccc                                                                                                                                                 <  1:7761/251‑1 (MQ=255)                                                                                                   ggCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTGCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCca                                                                                                                                                <  1:328795/251‑1 (MQ=255)                                                                                                            aTAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGt                                                                                                                                       <  1:244716/251‑1 (MQ=255)                                                                                                              aGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGtt                                                                                                                                      >  1:273848/1‑250 (MQ=255)                                                                                                                      ggATTTCGATGATGCGGTCGCCAAAGAAGCAGATCTGGAAAAAGTGAAACAGTGGCATCTGGAATATTACGGTGAGGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCGGCACCGCTGGCCTCGGTAGCCTCTTCCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGtttt                                                                                                                             <  2:410378/251‑1 (MQ=255)                                                                                                                             gatgatGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGAc                                                                                                                      >  2:176405/1‑251 (MQ=255)                                                                                                                                   gCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTgcgc                                                                                                                >  1:374910/1‑251 (MQ=255)                                                                                                                                           aaaaGCAGCAGAGCTGGAAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGc                                                                                                        <  2:209190/250‑1 (MQ=255)                                                                                                                                                                   gAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGTACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGACATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCGGGCCGCGTTTAATGACCCACAAAACCGTTATGCGGTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCg                                                                                <  1:212181/251‑1 (MQ=255)                                                                                                                                                                             cATCTGGAATATTACTTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATTTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtc                                                                      <  2:222748/251‑1 (MQ=255)                                                                                                                                                                             cATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtc                                                                      <  2:365928/251‑1 (MQ=255)                                                                                                                                                                                                 tgaACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATg                                                                                                                                                     >  2:399214/1‑152 (MQ=255)                                                                                                                                                                                                 tgaACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATg                                                                                                                                                     <  1:399214/152‑1 (MQ=255)                                                                                                                                                                                                      cGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGACGGCCACATTATTTTTTACGTcgacg                                             >  2:239253/1‑251 (MQ=255)                                                                                                                                                                                                              ttcGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                     >  2:165247/1‑251 (MQ=255)                                                                                                                                                                                                                          aTGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATTCGTTCCAGGCCATGTTACCTGCACCGCTGGCCGCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAACTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                         <  2:313304/251‑1 (MQ=255)                                                                                                                                                                                                                          aTGGACAACATGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                         >  1:6421/1‑251 (MQ=255)                                                                                                                                                                                                                              acaacaTGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                     <  1:250264/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTGTCTGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCatat                <  1:176405/251‑1 (MQ=255)                                                                                                                                                                                                                                        tGGTTCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt           >  1:217114/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttCTGTCCGGGTAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTAt        >  2:266610/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ccGGGTAATGTTGCCAGATGTGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGCGTGGCCGACAGATCGTCGGCCCCATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt  >  2:242119/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      TATCACGCTGGTTGCTGATGCGCGTCATTCGACGCCGGAGCAGATTGATAAAATTGCCCAGGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAT  >  NC_000913/1850300‑1850790

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