Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	98.4%	197.7 / ‑4.3	61	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/31);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993286‑3993767                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             cATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTAc                                                                                                                                                                                                                                         <  2:303716/251‑1 (MQ=255)       aaaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGGTCTTTTTCGAAGAAACGCAAAACGAGAATGGCTTTAATATCTAAATTCTCGACGAAAGCAACCGGGGTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGGGACGTCAGTCGCTTCTAATCGTCATTGCATGACCGCTTTACCTACgg                                                                                                                                                                                                                                   >  1:330433/1‑251 (MQ=255)                          tGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGGATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGGGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTTCCCACGGCTCAAACTTCATCAACTTc                                                                                                                                                                                                                >  1:88840/1‑251 (MQ=255)                            cAGGTTGAAACCAATCACGTCAAATTACCGGCTGTGTTGGACGGCTTTGCCAGCGAAGGGAGCATCCAGTTCTTTTTCGAATAAACGCAAGACGAGAAGTGCTTTAATATCTACATTCTCGACGTAAGCAACCGGGTTGAGGTATATCTCCACTGCGAAGGCAGCAAAGAGGGGCTGGTACGTGGCGTCAGTCTCTTCGACTCGTCATCGCATGGCCGCGTTACCTACGGCTCAAGCTTCATCAGCTTCaa                                                                                                                                                                                                              <  1:349289/251‑1 (MQ=255)                                 tgacaCCAATCACGTCAAATTACCTGCGGAGGTGGACGGATTTGCCAGCGAAGTTATCAGCCAGTTCGTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCGACATTCTCGACGAAAGCAACAGTGTTGAGGTATATCACCACTGCGAAGGCAGCATAGAGGATCTGGTACGTGACGTCAGTCGCTTCTACTCTTCATCGCAGTACCGCTTTACCTTCGGCTCAAGCTTCATCTACTCCAACCTg                                                                                                                                                                                                          <  2:69943/246‑1 (MQ=255)                                            aCGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCATTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATAGCACCACTGCGAAGGCAGCAAAGAGGAGGTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAGGACCTCTTTACCTACGGCTCAGGCTTCATCAACTTCAACCTGCCGCAGTTCTAt                                                                                                                                                                                              <  1:411253/251‑1 (MQ=255)                                                            ggtggtGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGCAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGCCGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                                                                                                                                              >  2:49805/1‑251 (MQ=255)                                                            ggtggtGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                                                                                                                                              <  1:331563/251‑1 (MQ=255)                                                                gtggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATgg                                                                                                                                                                          >  2:13426/1‑251 (MQ=255)                                                                  ggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTGCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                                                                                                                                        <  2:333990/251‑1 (MQ=255)                                                                           tGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAg                                                                                                                                                                <  1:297669/250‑1 (MQ=255)                                                                                 cGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGCTGGTCGTGAACAGGAGAttc                                                                                                                                                         >  1:15313/1‑251 (MQ=255)                                                                                  taaGGGATCATCCAGTTATTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGACAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTc                                                                                                                                                                                                 >  2:69326/2‑210 (MQ=255)                                                                                  gAAGGTATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTCATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCCTCAACTTCCACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACCGGTGAttcc                                                                                                                                                        >  2:239292/1‑251 (MQ=255)                                                                                  gAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTc                                                                                                                                                                                                 <  1:69326/210‑1 (MQ=255)                                                                                     ggtatcatcCAGTTCTTTTTCGAAGAAACTCAAGACGAGAATGGCTTTAATATCTACATTCTCTACGAAAGCAACCGGGTTGAGGTATAACACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGATACGTCAGTCTCTTCTACTCTTCATCGCATGACCGCTTTACCGATGTCTCAATCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCTTGAACAGGTTAttccgtt                                                                                                                                                     <  2:317544/248‑1 (MQ=255)                                                                                          catcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAGCAGGTGATTCCGTTCCGca                                                                                                                                                <  1:106286/251‑1 (MQ=255)                                                                                             ccAGTTCTTTTTCGCAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAACGTTGCTGGTCGTGAACAGGTGATTCCGTTCCGccca                                                                                                                                              >  2:310999/1‑247 (MQ=255)                                                                                                      tttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCaa                                                                                                                                                                                                                    <  1:137147/171‑1 (MQ=255)                                                                                                      tttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCaa                                                                                                                                                                                                                    >  2:137147/1‑171 (MQ=255)                                                                                                      tttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCg                                                                                                                                    >  1:363797/1‑251 (MQ=255)                                                                                                          gaagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTGCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTaa                                                                                                                                <  1:374089/251‑1 (MQ=255)                                                                                                           aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCACCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGTTCGTGAACAGGAGATTCCGTTCCGCACAAAATCTCTCGGTAAc                                                                                                                               >  2:222580/1‑251 (MQ=255)                                                                                                                aCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgcc                                                                                                                          <  1:62153/251‑1 (MQ=255)                                                                                                                        cGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGcc                                                                                                                   >  1:91490/1‑250 (MQ=255)                                                                                                                           aaaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCTCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTTCTCGTCATCGCATGACCGCTTTGCCTACGGCGCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAATTTTGATGGTCGTGAACAGGTGATGCCGTTCCGTACAAAATCTATCGGTAACAAGCCTCCTTCCAATc                                                                                                               <  2:182510/250‑1 (MQ=255)                                                                                                                              tGGCTTTAATATCTACATTCTCAACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAgg                                                                                                            >  2:52283/1‑251 (MQ=255)                                                                                                                                      atatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTTACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGAt                                                                                                    <  1:223355/251‑1 (MQ=255)                                                                                                                                      atatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAGTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGAt                                                                                                    <  2:214312/251‑1 (MQ=255)                                                                                                                                           tACATTCTCGACGAAAACAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAACTCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgcc                                                                                               >  2:40986/1‑251 (MQ=255)                                                                                                                                             cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgccgc                                                                                             >  1:236959/1‑251 (MQ=255)                                                                                                                                               ttCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGCTACGCCGCt                                                                                            >  2:118788/1‑250 (MQ=255)                                                                                                                                                tctcAACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTa                                                                                           <  1:52283/250‑1 (MQ=255)                                                                                                                                                    gacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTAc                                                                                       >  1:50930/1‑250 (MQ=255)                                                                                                                                                                 gggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAAGCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCg                                                                         <  1:336314/251‑1 (MQ=255)                                                                                                                                                                 gggTTGAGGTATATATCCACGGCGAAGGCAGCAAAGAGGAGCTGGTCCGTGACGTCAGTCGTTTCTACGCGTCATCTCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTTTCGGTAACATGCCGCCTGCCAATCAGGATCACGTTACGCCGCTATTACAGCAATATTTTTCg                                                                         <  2:302863/251‑1 (MQ=255)                                                                                                                                                                    ttGAGGTATATCACCACTGCGACGACAGCGACGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACCGGTGCTTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtg                                                                       >  2:375336/1‑250 (MQ=255)                                                                                                                                                                     tGAGGTATTTCACCTCTGCGACGGCAGCACAGAGGAGCTGGTTCGTGACGTCAGTCGCTTCTACTCGCCATCGCATGCCCGCTTTACCTACGGCTCAAGTTTCATCAACTTCAACCTGCCGCAGTTCTATCCGATTGTTAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTCACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTGTCGtga                                                                      <  2:376898/250‑1 (MQ=255)                                                                                                                                                                     tGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgct                                                                     >  2:299878/1‑249 (MQ=255)                                                                                                                                                                      gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatg                                                                    <  1:101842/251‑1 (MQ=255)                                                                                                                                                                      gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatg                                                                    >  1:218930/1‑251 (MQ=255)                                                                                                                                                                       aGGTATATCACCACTGCGAAGGCAGCAAAGAGGACCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCACATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCCATATTTTTCGtgatg                                                                    >  1:41610/1‑250 (MQ=255)                                                                                                                                                                          tataTCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAAc                                                                 <  1:256877/250‑1 (MQ=255)                                                                                                                                                                               caccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCg                                                           <  1:49805/251‑1 (MQ=255)                                                                                                                                                                                   acTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGaa                                                        >  1:333035/1‑250 (MQ=255)                                                                                                                                                                                        gAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAGTCAGGATCACGATACGCCGCTATTACAGCAAGATTTTTCGTGATGAACGTGCCGGAAAGCGAg                                                  <  1:400667/251‑1 (MQ=255)                                                                                                                                                                                        gAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCACTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGCTGGTCGTGAACAGGTGATTACGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCg                                                    >  1:252037/1‑249 (MQ=255)                                                                                                                                                                                          aGGCAGCAAAGAGGAGCTGGTACGGGACGACCGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCACGCTCCATCAACTTCAACCTGCCGCAGTTCTATCCGCTTGTGCAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCCCTATTACCGCAATATTTTTCGGGATGAACGTGCCGGAAAGCGAGGc                                                >  2:414477/1‑251 (MQ=255)                                                                                                                                                                                              cgcaAAGAGGAGCTGGTACGGGACGTCATTCGCTTCTACTCGTCATCGCAGGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAt                                            <  2:330433/250‑1 (MQ=255)                                                                                                                                                                                              agcaAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAt                                            >  1:61982/1‑251 (MQ=255)                                                                                                                                                                                                caAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATcc                                          >  2:346327/1‑251 (MQ=255)                                                                                                                                                                                                caAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATCACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATcc                                          >  1:85782/1‑251 (MQ=255)                                                                                                                                                                                                       gagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGc                                   >  1:4623/1‑251 (MQ=255)                                                                                                                                                                                                               aCGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                           <  1:181110/251‑1 (MQ=255)                                                                                                                                                                                                                            cGCTCCTACTCGTCATCGCATGACCGCTTTTCCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACATGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgact               <  2:447/250‑1 (MQ=255)                                                                                                                                                                                                                             gCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttt                                                                             >  2:317273/1‑187 (MQ=255)                                                                                                                                                                                                                             gCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttt                                                                             <  1:317273/187‑1 (MQ=255)                                                                                                                                                                                                                                   aCTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAc        <  1:230626/250‑1 (MQ=255)                                                                                                                                                                                                                                   aCTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAAGCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTGCGTGATGAACGGGCCGGAAAGCGAGGCTTAGCCGGCATGCAATCTTAGCGGAAACTGACTTTGTCAcc       <  2:20292/251‑1 (MQ=255)                                                                                                                                                                                                                                    cTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATCTTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGCCTGTTTCAccc      >  2:119414/1‑251 (MQ=255)                                                                                                                                                                                                                                         cATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGAGTGTGAAGGTTGATGGTCGTGAACAGGGGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTGATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct  <  1:239292/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             CATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993286‑3993767

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: ‑