Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,622	T→C	D496D (GAT→GAC)	rdcA →	regulator of diguanylate cyclase RdcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,622	0	T	C	100.0%	145.8 / NA	49	D496D (GAT→GAC)	rdcA	regulator of diguanylate cyclase RdcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/29);  total (20/29)

TCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTG  >  NC_000913/4328395‑4328866                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     tCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgccgc                                                                                                                                                                                                                               >  1:302129/1‑251 (MQ=255)    gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCt                                                                                                                                                                                                                             <  1:117139/250‑1 (MQ=255)        cgcATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTCTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGttt                                                                                                                                                                                                                        <  1:9799/251‑1 (MQ=255)             aaaCTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAgg                                                                                                                                                                                                                    >  1:378337/1‑250 (MQ=255)                     gAATTACGCGCAGCAGGCGCGGGAATAACTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCCGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGAGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCCGTTAGCCGCCTTGTTTCCGGATGATTCgg                                                                                                                                                                                                           >  2:281832/1‑251 (MQ=255)                                ggcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCCCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGc                                                                                                                                                                                                <  2:229587/250‑1 (MQ=255)                                   agGCGCGGGAATACCTGGATTGTCTTGCGCACGGGATAAACGTTGCTTGTGAACAATTGCGGCAAAATATCCACCAGAGCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCTCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCgag                                                                                                                                                                                             <  2:420426/251‑1 (MQ=255)                                    ggCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCg                                                                                                                                                                                               >  1:122509/1‑248 (MQ=255)                                       gcgGGAATACCTGGCTTCTCGTGCGCACGGGTTATACGACGCTTGTGACCAATTGCGGCAAAATTTCCAACAGATCGAAGAAAGTTTGCCGCTATTGCAACTCAATCCGGCGCAGGTGAGCGGCGAGATTAAACATGAAATAGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTACGTCAACCGCAAGACGCGCTGAAGGAGCAGTAAGCAGCCTTGTTTCAGGCTGATTCGGAGCCGTTAAGCgagaatc                                                                                                                                                                                         >  2:482637/1‑248 (MQ=255)                                        cgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGATATTCg                                                                                                                                                                                        <  1:98136/251‑1 (MQ=255)                                                ccTGGATTGTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATGGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACGCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCTCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCt                                                                                                                                                                                <  2:4907/251‑1 (MQ=255)                                                   ggATTTTCGTGCGCACGGGTTAAACGTCGCTTGGGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGTCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCtgtg                                                                                                                                                                             <  1:476603/251‑1 (MQ=255)                                                             gcgcACGGGTTAAACGTCGCTTGTGAACAAATGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCGATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGtt                                                                                                                                                                   <  1:115241/251‑1 (MQ=255)                                                               gtaCGGGTTAAACGTCGCTTGTGAACAATGGCGGCAAAATATCCACCAGATCAAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAAGCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTTGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTAc                                                                                                                                                                 <  2:63936/249‑1 (MQ=255)                                                                      ttAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGgc                                                                                                                                                          >  1:104356/1‑251 (MQ=255)                                                                           cGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCACCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGgcgc                                                                                                                                                        >  2:95302/1‑248 (MQ=255)                                                                                    tgAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATc                                                                                                                                             >  1:63300/1‑250 (MQ=255)                                                                                       aCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTc                                                                                                                                         <  1:158417/251‑1 (MQ=255)                                                                                        cAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGTGCAGAACACCATCAGTcg                                                                                                                                        >  1:199768/1‑251 (MQ=255)                                                                                                      tatCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGACGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGCCACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGc                                                                                                                          >  2:134737/1‑251 (MQ=255)                                                                                                       atCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGc                                                                                                                          <  1:408610/250‑1 (MQ=255)                                                                                                             caGATCGGAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTTAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGc                                                                                                                   >  2:305030/1‑251 (MQ=255)                                                                                                                              gCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCTTCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCAccc                                                                                                  <  1:203243/251‑1 (MQ=255)                                                                                                                                aGCTATTGCACCTCAATCAGTCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGATAGCCGCCTTGTTTCAGGATGCTTCGGGGCCGTTAAGCGAGAATCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGGCTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCtt                                                                                                >  2:313456/1‑251 (MQ=255)                                                                                                                                  cTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACAATGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTg                                                                                              >  1:471449/1‑251 (MQ=255)                                                                                                                                          aCTCAATCAGTCGCAGGTGAGCGGCGAGTTTAAACATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgt                                                                                        <  2:158947/249‑1 (MQ=255)                                                                                                                                          aCTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgtgt                                                                                      <  1:259865/251‑1 (MQ=255)                                                                                                                                              aaTCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAAc                                                                                  <  1:281832/251‑1 (MQ=255)                                                                                                                                              aaTCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACACGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAAc                                                                                  >  1:249471/1‑251 (MQ=255)                                                                                                                                                   ggCGCAGGTGGGTGGCGAGATTAAACACGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGGTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATg                                                                             <  2:337617/251‑1 (MQ=255)                                                                                                                                                   ggCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATg                                                                             <  1:386919/251‑1 (MQ=255)                                                                                                                                                            gAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcc                                                                    >  1:366358/1‑251 (MQ=255)                                                                                                                                                            gAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcc                                                                    <  1:94806/251‑1 (MQ=255)                                                                                                                                                                 gcgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCAtt                                                               >  1:220560/1‑251 (MQ=255)                                                                                                                                                                    agaTTAAACATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTGTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGAGGAGTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAg                                                            <  2:489395/251‑1 (MQ=255)                                                                                                                                                                       ttAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGca                                                          >  1:150524/1‑250 (MQ=255)                                                                                                                                                                                    cGAGCTGGCCCTGACCTCCGCCACCCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATg                                             <  2:472249/250‑1 (MQ=255)                                                                                                                                                                                          ggCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGGAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGGCTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGa                                      <  1:325195/251‑1 (MQ=255)                                                                                                                                                                                                     ccGCCAACCTCTTTCTGCGTCAACAACAAGACGCTCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGTCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGGCAggg                           <  1:332659/251‑1 (MQ=255)                                                                                                                                                                                                       gCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAACATGGAATTGAGCGAGTCAGGGtt                         >  1:73781/1‑251 (MQ=255)                                                                                                                                                                                                       gCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt                         <  1:42919/251‑1 (MQ=255)                                                                                                                                                                                                        ccAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGttt                        <  1:95302/251‑1 (MQ=255)                                                                                                                                                                                                        ccAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGttt                        <  1:145896/251‑1 (MQ=255)                                                                                                                                                                                                         cAACCCCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTc                       <  2:104356/251‑1 (MQ=255)                                                                                                                                                                                                                 ttCTGCGTCAACATCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCTTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACTGCTCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCGTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGTAATTGAGCGAGTCAGGGTTTCGTCCTg                 <  1:352050/249‑1 (MQ=255)                                                                                                                                                                                                                    tGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGGCAGGGTTTCGTCTTGGGttt            >  1:149896/1‑251 (MQ=255)                                                                                                                                                                                                                       gTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTGACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAt         <  1:199931/251‑1 (MQ=255)                                                                                                                                                                                                                         cAACAGCAAGACGCGCTGAAGGTGAAGCTAGCCGCCTTGGTTCAGGATGATTCGGCGCCGTTACGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGTTGTTATTCGCCCCATTGAGCCACAAGTAAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAttt       >  2:246314/1‑251 (MQ=255)                                                                                                                                                                                                                              gCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTGTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTg  <  1:41839/251‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTG  >  NC_000913/4328395‑4328866

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