Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,332	T→A	E229V (GAA→GTA)	dgcQ ←	diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,332	0	T	A	95.6%	217.2 / ‑0.2	70	E229V (GAA→GTA)	dgcQ	diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (37/28);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

CGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAT  >  NC_000913/2027085‑2027567                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            cgcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCg                                                                                                                                                                                                                                          >  2:251453/1‑251 (MQ=255)  gcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCgg                                                                                                                                                                                                                                         >  1:93394/1‑251 (MQ=255)              aGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGaa                                                                                                                                                                                                                             <  2:249926/251‑1 (MQ=255)              aGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACGGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGaa                                                                                                                                                                                                                             <  2:410351/251‑1 (MQ=255)               gTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCGCAGCTTACTGTCATAGAGCTGAGACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATGCCAGTAGTGATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTGCAGTAATGCTGTCGTGGACCGGTTGCGAAGTGAAc                                                                                                                                                                                                                            <  2:224408/251‑1 (MQ=255)                 tCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGCTTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAAcc                                                                                                                                                                                                                           >  1:229273/1‑250 (MQ=255)                   aCGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGc                                                                                                                                                                                                                        >  1:266126/1‑251 (MQ=255)                        gATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCCATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTAcc                                                                                                                                                                                                                   >  2:500248/1‑251 (MQ=255)                         aTCAATAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGTTCAAAAATCTCAGCTTCCTGTCTTAGCGATGATACTCACCATCGAGGTTTTTATCGATGCCGTTTCTTAAAAATTGCTGCCTAGTACGCACGGGAATACTCATCCCCAGCACGCCATACAAGTCGTTATTACTGTCTACCGGACAACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTCCCGGTTGCTAAGTGAACCAGCGTACcg                                                                                                                                                                                                                  >  2:540977/1‑251 (MQ=255)                               aaTATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCggg                                                                                                                                                                                                            >  1:171276/1‑250 (MQ=255)                                atatTCCCTGTTGGAGGGTCAGGATTGGAAGAGGACAAAAATCGCAGCTTACTGTCATAGAGCTGATTCTCACCATCGAGGTTTTTATCGAGGCCGTTTCTTAAAAATTGCTGCACAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGGAGTTATTACTGTCTCCCGTACCACTGACGGTAACCTGCGGTTCAATATTGCGGGCGTGTACCGGTTGCGAAGTGAACCAGGGTGCCGCTCGGTg                                                                                                                                                                                                           <  2:362967/251‑1 (MQ=255)                                    tCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACg                                                                                                                                                                                                        >  2:483440/1‑250 (MQ=255)                                            ggATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCg                                                                                                                                                                                               >  1:366128/1‑251 (MQ=255)                                              aTGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGtt                                                                                                                                                                                             >  1:271312/1‑251 (MQ=255)                                              aTGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAAGTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGtt                                                                                                                                                                                             <  1:415577/251‑1 (MQ=255)                                                    aGGATTGGAAGAGGTCAAAAAGCTCAGCTTACTGTCATTGAGCTGATACTCACCATCAAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATATTTCGCACGGGAATACTCATCCCCAGCACGCCAGACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAAt                                                                                                                                                                                       <  1:339326/251‑1 (MQ=255)                                                     ggATTGGAAGAGGTCAACAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATg                                                                                                                                                                                      >  2:41134/1‑251 (MQ=255)                                                      gATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGCTTTTCTCGTTGTGAATgg                                                                                                                                                                                     >  1:33602/1‑251 (MQ=255)                                                         tGGAAGAGGTCAACAATCTCAGCTTACTGCCATAGAGCTGATACTCACCATCGAGGTTTTTATCGCTGGCGTTTCTTAACAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGACCACTGACGGTAACCTGCGGTTCAGACTTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGAACCGCGCGGTGACGATTTTCTCGTTGTGAATGGcc                                                                                                                                                                                   >  2:85777/1‑250 (MQ=255)                                                         tGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTTCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGcc                                                                                                                                                                                   <  1:277734/250‑1 (MQ=255)                                                              gaggtcaaaactctaaaccaaccgTCATCGAGCTGACACTCACCCTCGAGGATTTTATCGATGGCGTTTCTTAAAAATTCCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACAGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATaa                                                                                                                                                                              >  2:337705/24‑250 (MQ=255)                                                                ggTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAAcc                                                                                                                                                                           >  2:142715/1‑251 (MQ=255)                                                                            cAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTg                                                                                                                                                                >  1:327828/1‑250 (MQ=255)                                                                                       cATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGta                                                                                                                                                     >  2:504147/1‑250 (MQ=255)                                                                                       cATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGCAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGAGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGta                                                                                                                                                     >  1:75446/1‑250 (MQ=255)                                                                                         tAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGCTTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGGGACATAGCCGtaat                                                                                                                                                   >  2:312243/1‑250 (MQ=255)                                                                                                      cACCATAGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCCACGTACCGCGCGGTGACGGTTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGCCATAGCCGCAATAACGCGTTGGAAc                                                                                                                                      >  2:254030/1‑250 (MQ=255)                                                                                                            cGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTgcg                                                                                                                                <  1:142715/250‑1 (MQ=255)                                                                                                                gTTTTTATCGATGGCGTTGCTTAAAAAGTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATATCCGTAATAACGCGTTGGAACATTGCTCGTaa                                                                                                                           <  2:181898/251‑1 (MQ=255)                                                                                                                 tttttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTaa                                                                                                                           >  1:280029/1‑250 (MQ=255)                                                                                                                   tttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCTTTGGAACATTGCGCGTAAAc                                                                                                                         <  1:112059/250‑1 (MQ=255)                                                                                                                     tATCGATGGCGTTTCTTTAAAATTGCTGCTTAGTACGCACGGGAATACTCATCCCCAGCACGCCATACAAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGAGACCCAGCGTACCGCGCGGAGACGATTTTCTCGTTGTGATTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAAAGCGTTGGACCATCGCGCGTAAACa                                                                                                                        >  1:456757/1‑249 (MQ=255)                                                                                                                     tATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGCGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACaa                                                                                                                       <  1:365189/250‑1 (MQ=255)                                                                                                                         gATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATCCTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACAGAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCg                                                                                                                  <  2:266126/251‑1 (MQ=255)                                                                                                                                  tCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTc                                                                                                          <  1:41134/250‑1 (MQ=255)                                                                                                                                             gctgcATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTGCCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAt                                                                                              <  1:526204/251‑1 (MQ=255)                                                                                                                                             gctgcATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAACCGTAAAAt                                                                                              >  2:427398/1‑251 (MQ=255)                                                                                                                                                   tAGTACGCACGGGAATACTCATCCCCGGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCCATAAACCAg                                                                                                                                                                         >  2:49869/1‑170 (MQ=255)                                                                                                                                                   tAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAg                                                                                                                                                                         <  1:49869/170‑1 (MQ=255)                                                                                                                                                        cGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAAACAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAcgcg                                                                                    >  2:366012/1‑250 (MQ=255)                                                                                                                                                            cgggaacaccaaTCCCCAGCACGCCATAACAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGGTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATGTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGGAATAACTCGTTGGAACATTGCGGGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCTCGCGAGaca                                                                               <  2:43383/240‑1 (MQ=255)                                                                                                                                                                   aCTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATc                                                                         <  1:251453/250‑1 (MQ=255)                                                                                                                                                                    cTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGGGACATAGCCGTAATAACGCGTTGGAACATTGCGCGGAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCCCGCGAGCAATACATCGc                                                                       >  2:304477/1‑251 (MQ=255)                                                                                                                                                                     tCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCGCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATc                                                                         <  1:552813/248‑1 (MQ=255)                                                                                                                                                                     tCATCCCAAGCACGCCCTACCAGTAGTTATTACTGTCTACAGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAAAGTACAGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACCTAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacatacctc                                                                         >  2:552813/1‑245 (MQ=255)                                                                                                                                                                     cctcccccATCACGCCATACCAGTTGTTATCACTGTCTACCTGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAATTGAGCCAGCGTACCGCGCGTTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCTTTGCAACATTGCGCGTCAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAcgcggttcattcatcgcc                                                                      <  2:37480/247‑15 (MQ=255)                                                                                                                                                                            catcACGCCATACCAGAAGTTATTACTGTCTACCGGATCACGTACTGGAACCTGCGGTTCAGTATTGCTGGCGGGTTCCGTTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGTAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCaa                                                               <  1:233102/248‑1 (MQ=255)                                                                                                                                                                                 cGCCATACCATTAGTTATGTCTGTCGACCGGAACACCGCAGGTAACCTGCGGCTCAGTATTGCTGGCGTGTACCGGTTTCGATGTGAACCTGCGTACGGGTTGGTGACGATTTTCTCGTTGCGAATGGCCAATAAACCAGGGGTGGGAGCCATGGCCGTAATAACGCGTTGGAACATTGCGCGGAAATAAGGTCGGCTGCGGAGAAACTTAAAATCCGGCCCGCGAGacatacatatcatggtcaaacatc                                                          <  2:385767/251‑17 (MQ=255)                                                                                                                                                                                     aTACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacg                                                      <  2:380476/251‑1 (MQ=255)                                                                                                                                                                                             agtTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTtg                                              >  1:261114/1‑251 (MQ=255)                                                                                                                                                                                                tattaCTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGCGGAGTTGTgc                                           >  1:349254/1‑251 (MQ=255)                                                                                                                                                                                                 attaCTGTCTACCTGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTgcg                                          <  1:535706/251‑1 (MQ=255)                                                                                                                                                                                                    aCTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCt                                        >  1:194634/1‑250 (MQ=255)                                                                                                                                                                                                      ggTCTGCCGGAACACTGACGGTACCCTGCGGTTCAGGATTGCTGGCGTGTGCCGGTTGCGAATTGACCCAGCGTACCGCGCGGTGACGATTTGCTCGTGGTGAATGTCCAAGAAACCAGGGTTGGGTGACATGGCCGTAATAACGCGGTGGAGCGTTGCGCGTAAACAAGGTTGGTTGGGTCGTAACATAAAATCCGGGTCGCAAGACATGCATCGCCGGTTCAACCATCGACGAGGAGTTGTGCGCTAGt                                     <  1:337705/250‑1 (MQ=255)                                                                                                                                                                                                           aCCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCgc                                  >  1:543524/1‑249 (MQ=255)                                                                                                                                                                                                             cGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGCGGAGTTGTGCGCTAGTCgcagc                               >  1:53824/1‑250 (MQ=255)                                                                                                                                                                                                               gAACACTGACGGTAACCTGCGGTTCAGTATGGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAg                            <  1:254030/251‑1 (MQ=255)                                                                                                                                                                                                                  cacTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTaa                         <  1:483440/251‑1 (MQ=255)                                                                                                                                                                                                                     tGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGGCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTCGTCGCAGCAAGTAAcc                       >  1:222420/1‑250 (MQ=255)                                                                                                                                                                                                                     tGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGGTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTAGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAAcc                       >  1:483184/1‑250 (MQ=255)                                                                                                                                                                                                                      gACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAAcc                       <  1:504147/249‑1 (MQ=255)                                                                                                                                                                                                                      gACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCa                      >  1:133554/1‑250 (MQ=255)                                                                                                                                                                                                                       aCGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAAc                    >  1:456194/1‑251 (MQ=255)                                                                                                                                                                                                                        cGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACt                   <  2:327828/251‑1 (MQ=255)                                                                                                                                                                                                                         ggTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGCGGAGTTGTGCGCTAGTCGCCGCAAGTAACCAACtt                  >  2:51037/1‑251 (MQ=255)                                                                                                                                                                                                                            aaCCTGCGGTTGAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTcc                <  1:77827/250‑1 (MQ=255)                                                                                                                                                                                                                               cTGCGGTTCAGCATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTg            >  1:359770/1‑251 (MQ=255)                                                                                                                                                                                                                                tGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGAACCGCGCGGGGACGATTTTCTCGTTGTGAATGGCCAATACACCAGGGTTGGGTGAAATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGGAGGCTGCGTCGAAACGTAAAATCCGGACCGCGAGACATACATCCCCTGTTCAACCAACGACGAGGAATTGTGAGCTAGTCGAAGCAAGAAACCAACATCCAGTGc           >  1:346084/1‑251 (MQ=255)                                                                                                                                                                                                                                ggCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGc           <  1:7891/250‑1 (MQ=255)                                                                                                                                                                                                                                         gTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAt  <  1:427398/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            CGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAT  >  NC_000913/2027085‑2027567

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