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	91.5%	123.6 / 2.4	47	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (23/20);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

AAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCA  >  NC_000913/3993292‑3993760                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      aaaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGGGGGGGAAGGCTTTGCCAGCGAAGGGATAATCCAGTTCTTTTTTGAAGAAACGCAAAACGAGAATGGCTTTAATATCTACATTCTCGACGAAAACAACCCGGTTGAGGTATATCCCCACTGCGAAGGCAGCAAAAAGGAGCTGGTACGTGAAGTCAGTCCCTTCTACTCGTCATCGCATGAACGCTTTACCTAccg                                                                                                                                                                                                                            >  1:111438/1‑249 (MQ=255)      gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAATGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCa                                                                                                                                                                                                                        <  1:423436/250‑1 (MQ=255)        aCGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGGCGGCCTTGCCCGCGAAGGGATCCACCCGTTCTTTTTTGGAGAAAAGCAAGACGAGCATGGGTTTAATATCTACAATCTCGGCGAAAGCCACCGGGGTGAGGAATATCACCACTGAGAAAGCAGACAAGAGGAGCTGGTACGTGGCGTCCGTCGCTTCCACTCGTCATCGCATGACCGCTTTAACTACGGCTCacac                                                                                                                                                                                                                     >  1:22370/1‑248 (MQ=255)              tGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGTCGGCTTTGCCAGCGAAGGGGTCATCCAGTTCTTTTTCCAAGAAACGCAAGACGAGAATGGCTTTAATATATACATTCTCGACGAAAGCAACCGGGGTGAGGGATATCACCCATGCGAAGGCAGCAAAGAGGAGCTGGGACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGcatcctc                                                                                                                                                                                                               >  1:311802/1‑245 (MQ=255)                     gCAGGTTTAAACCAATCACGTCAAATTTCCTGCGGTTGTGGACGGCTTGGCCAGCGAAGGGATTATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTc                                                                                                                                                                                                         <  2:38289/250‑1 (MQ=255)                                    ccACGTCAAATTACCGTGGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTGTTCGAAGAAACGCAAGTCGAGAATGGCTGTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCt                                                                                                                                                                                         <  2:120087/250‑1 (MQ=255)                                      aCGTCAAATTACCGGCGGTGGTGGACGGCTTTTCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAAAGCAAGACGAGAATTGCTTTAATATATACATTCTCGACGAAAACAACCGGGGTGAGGTATATCACCCCTGCGGAGGCCGCAAAGAGGAGCTGGTACGTGGCGTCAGTCGCTTCTACTCGTCATCGCATGACCGCCTTACCTACGGCTCAAGCTTCATCCACTTTAACCAGCCGccggtcta                                                                                                                                                                                        >  1:407159/1‑243 (MQ=255)                                                 cgggggggggtggACGGATGTGCCAGCGAAGGGATCATCCAGTTCTTTTGCGAAGAAACGCAAGGCGAGAATGGCTTTAATAGCTACATTATCGCCGAAAGCAGCCGGGTTGAGGTAGAGCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCGTCTACTCGTCAGCGCATGACCGCGTTACCTACGGCGCAAGCTTCAGCAACTTCAACCGGCTGCAGTTCGATCATATTGTGaa                                                                                                                                                                            <  2:222915/243‑1 (MQ=255)                                                 ccGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCCAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGa                                                                                                                                                                             >  2:120837/1‑250 (MQ=255)                                                 ccGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGaa                                                                                                                                                                            <  1:282213/251‑1 (MQ=255)                                                   ggcggTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAgg                                                                                                                                                                          >  1:143463/1‑251 (MQ=255)                                                              aCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCCTAGCATGACCGCTTTACCTACGGCTCAAGATTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                                                                                                                               >  1:281343/1‑251 (MQ=255)                                                                     tGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATAAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAg                                                                                                                                                         <  1:193789/250‑1 (MQ=255)                                                                           cGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                                                                                                                                  >  1:415752/1‑251 (MQ=255)                                                                                    catcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGc                                                                                                                                          >  2:178646/1‑250 (MQ=255)                                                                                         aGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTAAGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGCTTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaa                                                                                                                                     >  2:426820/1‑250 (MQ=255)                                                                                                 ttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                                                                                                                                 <  1:180728/214‑1 (MQ=255)                                                                                                 ttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                                                                                                                                 >  2:180728/1‑214 (MQ=255)                                                                                                      agaaACGCAAGAAGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTACGTTCCGCACAAAATATATCgggaaa                                                                                                                        >  2:320205/1‑246 (MQ=255)                                                                                                                  cGAGAAGGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGAACAAAATCTATCGGTAACCTTCCGCCTGcc                                                                                                            >  2:271304/1‑250 (MQ=255)                                                                                                                         ggCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGGTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAgg                                                                                                     <  1:426820/250‑1 (MQ=255)                                                                                                                            tttAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTAAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATc                                                                                                  >  1:448683/1‑250 (MQ=255)                                                                                                                               aaTATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGa                                                                                              <  1:160457/251‑1 (MQ=255)                                                                                                                                        aTTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgcc                                                                                                                   <  2:234325/221‑1 (MQ=255)                                                                                                                                        aTTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgcc                                                                                                                   >  1:234325/1‑221 (MQ=255)                                                                                                                                              gacgaAAGCAACCGGGTTGAGGTATATCACCCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCCGTCGCTTCTACTCGTCATCGCATTACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATGGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACa                                                                               <  2:151613/251‑1 (MQ=255)                                                                                                                                                   aaGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGACGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACTAAATCTATCGGAAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAat                                                                          >  1:458216/1‑251 (MQ=255)                                                                                                                                                      aaaCCGGGTTGAGGGATAGCACCACTGCGAAGGCAGCAAAGAGGATCTGTCACGGGACGTAAGTCGATTCTACTCGTCTTCGCATGCCCGTCTTTCCTAAGGCTCAAGCTTCATCAACTGCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGGTCCGCACAAAATCTCTCTGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAtt                                                                       <  2:50477/250‑1 (MQ=255)                                                                                                                                                                gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGAGAAGGTTGATGGTCGTGAACCGGTGATTCCGTTCCGCACACAATCTATCGGTAACATGCCGCCTGCCAATCAGGAACACGATAAGCCGCTATTACAGCAATATTTTTCGttatg                                                             >  2:245809/1‑251 (MQ=255)                                                                                                                                                                     atatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACg                                                         <  1:48135/250‑1 (MQ=255)                                                                                                                                                                       cacaccccTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTTCCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATGGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGc                                                      <  2:96849/249‑1 (MQ=255)                                                                                                                                                                          cccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGGCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCgg                                                   <  1:228899/250‑1 (MQ=255)                                                                                                                                                                                gCGAAGGCAGCAAAGAGGCGTGGGGACGTGACGCCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCGCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTTAAAAGGTGATTCCGTTCCGCACAAAAGCTATCGGTAACAGGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATTATCGTGCCGGAAAGCg                                             <  2:284054/251‑1 (MQ=255)                                                                                                                                                                                           aaaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAt                                                                                                                      >  1:454479/1‑167 (MQ=255)                                                                                                                                                                                           aaaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAt                                                                                                                      <  2:454479/167‑1 (MQ=255)                                                                                                                                                                                            aaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGCGGCTTATCCgg                                 >  1:270119/1‑251 (MQ=255)                                                                                                                                                                                            aaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCgg                                 >  1:383367/1‑251 (MQ=255)                                                                                                                                                                                              gaggagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCa                               <  1:307717/251‑1 (MQ=255)                                                                                                                                                                                               gggtgcgggtacggtaCGTCAGTCGCTTCTTCTCGTCATCGCTTGCCCGCTTTACCGACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTTAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCAt                              <  2:273568/236‑1 (MQ=255)                                                                                                                                                                                                gggccTGGTACGGTACGTCATCCGCTTCTACTCGGCATGGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACGTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGGCGTGAACAGGTGATTCCGTTCCGCACAAAAGCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATg                             <  2:457998/247‑1 (MQ=255)                                                                                                                                                                                                           gTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCg                  >  1:308409/1‑251 (MQ=255)                                                                                                                                                                                                             gACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCgg                 >  1:123700/1‑250 (MQ=255)                                                                                                                                                                                                             gACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCgg                 >  1:424596/1‑250 (MQ=255)                                                                                                                                                                                                                    gTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgac         <  1:454664/251‑1 (MQ=255)                                                                                                                                                                                                                         ttCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCGGCCAATCAGGATCACGATACGCCGCTATTCAAGCACTATTTTTCGTGATAAAGTTGCCGGAAAGAGCGCCTTATCCGGCATGCAATCTTAgagagaaacgatgattt    >  2:34605/1‑235 (MQ=255)                                                                                                                                                                                                                           caccTCGGCATCGCATGACCGCTTTACCTACTGCTCAAGCTGCATCAACTTCAACCTGCCGCAGTGCTAGTAGATTGGGAAGGGTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTTCCGGTAACATGCCGCCTTCCAATCAGGATCACGATACGCCGCTATTACATCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAGCCGGCTTGCAATCTTAGCGGACTCTGACTGTTTc   <  2:281343/247‑1 (MQ=255)                                                                                                                                                                                                                             aCTCTTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTTCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCa  <  1:61989/249‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      AAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCA  >  NC_000913/3993292‑3993760

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