Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,595,570	(A)6→5	coding (282/915 nt)	lgoR ←	putative DNA‑binding transcriptional regulator LgoR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,595,570	0	A	.	100.0%	184.4 / NA	41	coding (282/915 nt)	lgoR	putative DNA‑binding transcriptional regulator LgoR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/18);  total (23/18)

AAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTA  >  NC_000913/4595327‑4595816                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aaaGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTc                                                                                                                                                                                                                                                >  2:307778/1‑251 (MQ=255)      gatGAAGTGAATGGGTTTCAGGCATTTCGCGTATCTCAAACAGTTGTTCGGCGTAGGATTGTACAAATTGCGTGATGCGCCACTGGCCGCGCGTGTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATAGTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAATTGGCGATGGTTTATCATCGTGAAAAA‑GCCTGCTCAAaca                                                                                                                                                                                                                                           <  2:371800/249‑1 (MQ=255)              gAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACtttgtt                                                                                                                                                                                                                                    <  1:90591/250‑1 (MQ=255)                         gAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAACCACTTTGTTTTGTTCACTCAt                                                                                                                                                                                                                        >  2:189850/1‑251 (MQ=255)                          aGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCAtt                                                                                                                                                                                                                       >  1:83384/1‑251 (MQ=255)                                  gcgTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAAACGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCAGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTAACTCATTGAAGCGgt                                                                                                                                                                                                               >  2:466277/1‑251 (MQ=255)                                   cgTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGgt                                                                                                                                                                                                               >  2:429369/1‑250 (MQ=255)                                              cAGCTGCTCGGCGTAGGATTGATAGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAGTAAGATTATAGCGGCCGAATTTTAAAAGATATTCTCTTACTACAAACGGGCTGACACCCGCGGCCCGAGCCAGTTGCCGTTCAGACAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGa                                                                                                                                                                                                    >  2:51475/1‑250 (MQ=255)                                                 ccgctCGGCTTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCTCTGTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCTTCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTTCAGGCGAAAcc                                                                                                                                                                                                <  2:277311/249‑1 (MQ=255)                                                                  gATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGtt                                                                                                                                                                               <  1:307778/251‑1 (MQ=255)                                                                               cATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAAc                                                                                                                                                                  <  1:189850/251‑1 (MQ=255)                                                                                     ccACTGGCCGCGATTTTAGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAT‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGGTTACGCGCAATAACATAGt                                                                                                                                                             <  2:157038/250‑1 (MQ=255)                                                                                      cACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCg                                                                                                                                                           >  2:358538/1‑251 (MQ=255)                                                                                                   gttcgttgagAATAAGAATATAACAGCCTAATTTTAAAAGATATTCTCTAACGACATCCGGGCGTACACTCGCGGCCCGAGCCAGGTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCCCAAAAACATAGTTGATTCCAACCTg                                                                                                                                                <  2:170869/240‑1 (MQ=255)                                                                                                          gtgAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGAc                                                                                                                                       <  2:81837/251‑1 (MQ=255)                                                                                                                   aTTATAACGGCCGAATTTTAAAAGATCTTCTCTTACTACAACCGGACTAACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTGGTACAGACGACACCGTCATCATGGTCCGGTTTCCGCGCAATACCTTAGTCGTTTCCAACCTGCGTCACAACGCCGCATTc                                                                                                                              >  2:296967/1‑251 (MQ=255)                                                                                                                        aaCGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGAGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGt                                                                                                                          >  2:418880/1‑250 (MQ=255)                                                                                                                                aaTTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTg                                                                                                                 <  1:94509/251‑1 (MQ=255)                                                                                                                                   attAAAAGATATTCTCTTACTACAACCGGGCTGACTCCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTAGCATCGTGAAAAA‑GCCTGCTCAAACACGTTGTTTTGTGCACACATTTAAGCTGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgag                                                                                                               <  1:301190/249‑1 (MQ=255)                                                                                                                                    ttAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgaga                                                                                                              >  1:389930/1‑250 (MQ=255)                                                                                                                                    ttAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAACTGGCTgaga                                                                                                              >  2:392879/1‑250 (MQ=255)                                                                                                                                              ttCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACtttgtt                                                                                                                                                                                                                                    <  1:48465/122‑1 (MQ=255)                                                                                                                                              ttCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACtttgtt                                                                                                                                                                                                                                    >  2:48465/1‑122 (MQ=255)                                                                                                                                                          aaCCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAACCGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGCAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCACGACGCCGCATTCGCGTAAGTGGCTGAGAATATGTCGCCCTGTGGTGCGGCt                                                                                       >  2:155004/1‑251 (MQ=255)                                                                                                                                                                   gACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTa                                                                              <  2:83384/251‑1 (MQ=255)                                                                                                                                                                      acCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                                                                           <  1:278722/251‑1 (MQ=255)                                                                                                                                                                            ggCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                                                                     >  1:477868/1‑251 (MQ=255)                                                                                                                                                                                     cAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACt                                                            >  1:168084/1‑251 (MQ=255)                                                                                                                                                                                                         tttCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                                                                 <  2:26341/162‑1 (MQ=255)                                                                                                                                                                                                         tttCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                                                                 >  1:26341/1‑162 (MQ=255)                                                                                                                                                                                                          ttCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCGGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGGGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                        <  1:300549/250‑1 (MQ=255)                                                                                                                                                                                                              ccTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTACGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCCCTTTGCGTTGGAAGTGGGGAAGGGat                                    >  2:272246/1‑250 (MQ=255)                                                                                                                                                                                                              ccTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGata                                   <  2:23328/251‑1 (MQ=255)                                                                                                                                                                                                               cTGGGCGTAACTGGCGCTGGTTTATCATCGTGACAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCg                                                                                                                                                           >  1:362093/1‑130 (MQ=255)                                                                                                                                                                                                               cTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCg                                                                                                                                                           <  2:362093/130‑1 (MQ=255)                                                                                                                                                                                                                 ggtcGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTCGTACAGGCGAAACCGTAATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCAGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGCTGGAAGTGGGGGAGGGATATgg                                >  1:482761/4‑250 (MQ=255)                                                                                                                                                                                                                   gCGTAACTGGCGCTGGTTTATCATAGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGCAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGAGCACTAACATAGTCGTTTCCAACCTGCGTCAAGCCGCCGCATTCGCGTAATAGGCTGCGATTATGACGCACTGTGGAGCGGCTGTTGTTTTACATTTCCGCCAGCGCACTTTGCGGTGGAAGTGGGGagtggatgtgtac                              >  2:470888/1‑239 (MQ=255)                                                                                                                                                                                                                     gTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCAGTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcccc                            <  1:178113/251‑1 (MQ=255)                                                                                                                                                                                                                      tAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcg                           >  1:191161/1‑251 (MQ=255)                                                                                                                                                                                                                               cTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGAACCGGCGAAACCGTCATAATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCGGATGTTGTACATTTCCGCCAGCGCCCTTTGCGGTGGAAGTGGGGAGGGGATTTGATCCAGCGCCATatc                  >  2:219863/1‑251 (MQ=255)                                                                                                                                                                                                                                        catcGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACGCCGCGCCATATCATCAATGAc         >  2:83280/1‑251 (MQ=255)                                                                                                                                                                                                                                            cgtgaaaaaGCCTTCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACAGAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg      <  2:268029/247‑1 (MQ=255)                                                                                                                                                                                                                                                aaaaaGCCTGCTCAAACACTTTGTTTTGGTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACGTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTa  <  1:484159/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       AAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTA  >  NC_000913/4595327‑4595816

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