Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,836	T→C	E61G (GAA→GGA)	dgcQ ←	diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,836	0	T	C	97.9%	165.8 / NA	48	E61G (GAA→GGA)	dgcQ	diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/32);  minor base A (1/0);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCG  >  NC_000913/2027595‑2028076                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ggaCAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAgatga‑t                                                                                                                                                                                                                                         <  1:161586/250‑1 (MQ=255)     aGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAAATGGCAGAGTGCGTCGTCGGTTGCGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGCGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACg                                                                                                                                                                                                                                     >  2:178943/1‑251 (MQ=255)         gaTTCCCATCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACgctgg                                                                                                                                                                                                                                 >  1:109803/1‑251 (MQ=255)         caTTCCCCTCGTTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCTGCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTGGTCGCGATGCTGTTCGAACTCGGTTACAGCATCCCGCAGAGCGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACgctgg                                                                                                                                                                                                                                 <  2:228605/250‑1 (MQ=255)          agtCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACGGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACgctggc                                                                                                                                                                                                                                <  1:181606/249‑1 (MQ=255)               ccTCGCTTACGAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGgag                                                                                                                                                                                                                           <  2:327120/251‑1 (MQ=255)                  cGCTGACTAATGCAGCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGAGTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACGCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGAGCAGTTTGGCGACGTTATACTGCAAAAGCTTATCGAGCGCGTTGGCAACGTTGCCCAGATGA‑TTACGCTGGCTGgaga                                                                                                                                                                                                                          <  2:225403/249‑1 (MQ=255)                      ggcTAATTCATCCGACACACCGTTACCTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGATCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCATATTA‑TTACGCTGGCTGGAGATATAg                                                                                                                                                                                                                     <  2:311246/248‑1 (MQ=255)                             gCATCCGACACAACATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATGTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGGCGACGTTATACTGCAAATGCTTATCGAGCGCGGTGGCAACGTTTCCCAGATGA‑TTACGCGGGCTGGAGATATAGGCATCtt                                                                                                                                                                                                              <  2:207388/250‑1 (MQ=255)                                         ccATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTTCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGcaccaccac                                                                                                                                                                                                 <  1:178943/251‑1 (MQ=255)                                           aTTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTACAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACtt                                                                                                                                                                                               >  1:269812/1‑251 (MQ=255)                                                  ggCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAAGGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGTCTGGAGATATAGGCATCTTCCAGCACCTCCACTTCACGc                                                                                                                                                                                          <  1:47946/249‑1 (MQ=255)                                                       agTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTg                                                                                                                                                                                    >  1:416190/1‑250 (MQ=255)                                                            gtggtcgGTTGAGTTCGATTTTACAGGCGGGCTCGTAGCGAGGCTGTTCGAACTCGGTTACAGCATCACGCAGAGCGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTAAGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCGGGTGATATAGGCATCTTACATCACTCCCACTTCACGCCATGTGAGCa                                                                                                                                                                                >  2:319493/4‑249 (MQ=255)                                                            gtcgtcGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCa                                                                                                                                                                                >  1:392506/1‑249 (MQ=255)                                                                   gTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTGCGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaa                                                                                                                                                                        >  1:295730/1‑250 (MQ=255)                                                                      gAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaaaaa                                                                                                                                                                     >  1:177130/1‑250 (MQ=255)                                                                       ggtgctattgtcacggCGTGCGCGGCGCGATGCTGTTCGAACTCGGCGATCGTATCACGCAGAGAGGTGAAATCCAGTGCCGCTACGAGAGCTTCGCGCATGCCATTACGTGAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATAGAGCTCGTTGGCAACGTTGCCCAGATGA‑GTACGCTGGATGTAGATATAGGCATCTTCCAGCACCCCCACTTCCCGCCAGGTGAGCAAGGTGGAAAAAAg                                                                                                                                                                   <  2:10812/237‑1 (MQ=255)                                                                       aGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAg                                                                                                                                                                   <  1:332156/251‑1 (MQ=255)                                                                                                 gcgATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACTATAAAGCAGAGATTaa                                                                                                                                         <  1:149463/251‑1 (MQ=255)                                                                                                  cgATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAAc                                                                                                                                        <  2:209325/251‑1 (MQ=255)                                                                                                    aTGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGa                                                                                                                                      <  1:74933/251‑1 (MQ=255)                                                                                                                       acagtagcACGCAGAGAGGTGAAATCCAGGGGCGCTACGAGAGCTTCGCGCATGTCATTACGTTAAATGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGTGTGTTGGCAACGTTTCCCATATGA‑TTTCGCTGGCGGGAGATATAGGCATCTTCCATCTCGACCACTTCACTCCAGGTGAGCAAGGTGGATGAAAGCAATACCACGAGAAAGCAGAGATTAACGACATTACCTGGCGCCAggc                                                                                                                    <  2:165829/243‑1 (MQ=255)                                                                                                                              cACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCg                                                                                                            >  1:91306/1‑251 (MQ=255)                                                                                                                                  cAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGAGCACGCCATTACGTAAAAAGATCAGTTTGTCGACGTAATACTGCAACTGCTTATCGAGCGCGTTGGCAACGATTACCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGACAGGTGAGCAAGGTGAACAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtt                                                                                                        >  1:421917/1‑251 (MQ=255)                                                                                                                                      gagGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttttt                                                                                                    >  2:309828/1‑251 (MQ=255)                                                                                                                                                     ggCGCTACGAGATCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt                                                                                      <  1:313983/250‑1 (MQ=255)                                                                                                                                                      gcgcTACGAGAGCTTCGCGCACGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTGTCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTGCGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGttt                                                                                    <  2:208174/251‑1 (MQ=255)                                                                                                                                                           aCGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCAt                                                                               >  1:36084/1‑251 (MQ=255)                                                                                                                                                                     cgcgCATGCCATTACGGAAAAAGACCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAGCGTTTCCCAGAGGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGTCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGtctc                                                                      <  2:277679/250‑1 (MQ=255)                                                                                                                                                                         cATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGtgct                                                                 >  2:301242/1‑251 (MQ=255)                                                                                                                                                                             ccATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTAGCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGTCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCAc                                                             <  2:159369/251‑1 (MQ=255)                                                                                                                                                                             ccATTACGTAAAAAGATCAGTTTGTCGACGGTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCAc                                                             <  2:320229/251‑1 (MQ=255)                                                                                                                                                                                      aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTgacgac                                                    >  2:179657/1‑251 (MQ=255)                                                                                                                                                                                                      cGACGTTACACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGAGTAACGTCATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTGTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCtttt                                     <  2:283403/250‑1 (MQ=255)                                                                                                                                                                                                        aCGTTCTACTGCAAATGCTTATAGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGCTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCCCCCTGACGACTCCCAACCTGCttt                                      >  2:388709/1‑247 (MQ=255)                                                                                                                                                                                                        aCGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCttt                                      <  1:388709/247‑1 (MQ=255)                                                                                                                                                                                                               actacAAATGCTTATCGAGCGTGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCGTCCAGCACCACCACTTCACGCAAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGTCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTGTTCCATTTTGGTCTCGTGCTGCACCCGGACGACTCCCAACCTGCTTTTTTATGATTc                            <  2:68798/246‑1 (MQ=255)                                                                                                                                                                                                               aCTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCCGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCt                           >  2:399611/1‑251 (MQ=255)                                                                                                                                                                                                                  gCAAAGGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCCCCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTGTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTGATGATTCTGGc                        <  2:425340/251‑1 (MQ=255)                                                                                                                                                                                                                   cAAATGCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCa                       <  2:177130/251‑1 (MQ=255)                                                                                                                                                                                                                     agtGCTTAACGGGCGCGTTGGCAACGTTTCCCAGGTGA‑TTGCGCTGGCGGGAGATATAGGCATCTTCCAGCACCACCACATCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATGAAGCAGAGACTAACGACATGACCAGGCACCAGGCGGCGTGCGAGGTTTTTCAACCCGCTCTAGTTTTCTATTTTTGTCGCGTGCTGCTCCCTGACGACTCCCAACCGGCTTTGTTATGATTCTGGCATa                     <  2:4210/249‑1 (MQ=255)                                                                                                                                                                                                                       agCTTATCGAGCGCGTTGGGAACGGTTCCCAGATGAGTTA‑GCTGGCTGGGGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAAc                   <  2:4447/250‑1 (MQ=255)                                                                                                                                                                                                                        gCTTATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑GTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAAcc                  <  2:91306/251‑1 (MQ=255)                                                                                                                                                                                                                           tATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGTTTAACGACATGACCAGGCCCCATGCGGCGTGCTAGTGTTTTCAACCAGCTCTGGTTTTCCATTGTTGGCTCGTGCTGCGCCCTGACGACTCCCAGCCTGCTTTTGTAGGATTCTGTCATAACCt                 <  2:262132/249‑1 (MQ=255)                                                                                                                                                                                                                           tATCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTg                <  1:81225/250‑1 (MQ=255)                                                                                                                                                                                                                            aTCGAGCGCGTTGGCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTGTTGTCTCGTGCTGCTCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTgg               <  2:157283/250‑1 (MQ=255)                                                                                                                                                                                                                                        ggCAACGTTTCCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCg  <  1:399611/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  AGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGA‑TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCG  >  NC_000913/2027595‑2028076

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