BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,907,998	A→T	5.2%	intergenic (‑57/+31)	eno ← / ← pyrG	enolase/CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,907,998	0	A	T	5.2%	837.5 / 20.6	264	intergenic (‑57/+31)	eno/pyrG	enolase/CTP synthetase
Reads supporting (aligned to +/- strand): ref base A (117/127); new base T (6/11); total (124/140)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

AAACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGG  >  NC_000913/2907850‑2908148

aaaCCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAg                                                                                                                                                       <  1:805183/150‑1 (MQ=255)

aaaCCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAg                                                                                                                                                       <  2:2906094/150‑1 (MQ=255)

    caccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGCGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCCGGGTa                                                                                                                                                   >  2:1827684/1‑150 (MQ=255)

      cccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcg                                                                                                                                                 >  2:2000552/1‑150 (MQ=255)

      cccTCCAGATGTACTTCGGCTTCAACAGACGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcg                                                                                                                                                 >  2:3480206/1‑150 (MQ=255)

        cTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                                                                               >  2:3026051/1‑150 (MQ=255)

         tCCAGATGTACTTCGGCTACAAGAGTCGGGTTACCTCGGGAGTCGATGATTTCTCGACAGTTGTTTTTTTAGATTTTGGACTTTAGGTTTTCCTCAAGTCACTAGTTCAACTGAAACTCCAGACAAACAACGATTACCCAGTGTACGCGt                                                                                                                                              <  2:3762485/150‑1 (MQ=255)

          ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGtt                                                                                                                                             <  1:4185398/150‑1 (MQ=255)

           cAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTg                                                                                                                                            <  2:3151003/150‑1 (MQ=255)

            aGATGTACTTCGGCTTCCACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGc                                                                                                                                           <  1:257325/150‑1 (MQ=255)

            aGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgt                                                                                                                                           <  1:678830/150‑2 (MQ=255)

            aGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgt                                                                                                                                           <  1:3981617/150‑2 (MQ=255)

            aGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTgt                                                                                                                                           <  2:139744/150‑2 (MQ=255)

            aGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTgt                                                                                                                                           <  2:2298137/150‑2 (MQ=255)

            aGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGc                                                                                                                                           <  1:1792012/150‑1 (MQ=255)

                gTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCg                                                                                                                                         >  1:3711309/1‑148 (MQ=255)

                gTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCg                                                                                                                                         <  2:3711309/148‑1 (MQ=255)

                    ttCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTaa                                                                                                                                   >  2:2800097/1‑150 (MQ=255)

                        gCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACttt                                                                                                                               <  2:256910/150‑1 (MQ=255)

                          ttCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACttttt                                                                                                                             >  2:3734196/1‑150 (MQ=255)

                            cAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATATTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                                                                           <  2:898263/150‑1 (MQ=255)

                             aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGCTGATTTTTACGCTTTTGGACCTTAGGTTTTCCGCAAGTCACTAGTTAAACTGAAACACCAGAAAAACAACGCCTACCCAGGGTACGCGTTGCCGCTCTAACtttttttc                                                                                                                          >  1:1219710/1‑148 (MQ=255)

                             aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                                                                          >  2:3226237/1‑150 (MQ=255)

                             aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAAATGAAACACAAGACAAACAAAGCGTACCCAGGGTACGCGTTTCCGCTCTAACTTTTttac                                                                                                                          >  1:2676447/1‑150 (MQ=255)

                                cgccGGGTAACCACGGGTGTCCTTTATTTCACGACCTATTTTTTTTACTATTTTGGACATTAGTTTTACCTAAAGACACTATTTAAACTCAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                       <  1:2839435/147‑1 (MQ=255)

                                aGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                       <  2:1268505/150‑1 (MQ=255)

                                      gTTACCAGGGGTGTCTATGATTTCACGACCGATGATTTTTACTTTTTTGTAATTTTGTTTTTCCTCAAGTCAATATTTAAACTTAAACTCCAGACAAACAACGCGTACCCCGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGc                                                                                                                 <  1:540173/150‑1 (MQ=255)

                                      gTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGc                                                                                                                 <  2:2877521/150‑1 (MQ=255)

                                           cACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                                                                            <  2:1745286/150‑1 (MQ=255)

                                           cACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                                                                            <  2:903332/150‑1 (MQ=255)

                                              gggAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTT‑AACTGAAACTCCAGACAAGCAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGttt                                                                                                        >  2:3080444/1‑150 (MQ=255)

                                               ggAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGttt                                                                                                        <  1:464209/150‑1 (MQ=255)

                                                 aGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCt                                                                                                      <  2:1548685/150‑1 (MQ=255)

                                                   tCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTAACTAGTTAAAATGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTAATTACTTCGCCTGACGTTTCTgg                                                                                                    >  2:164223/1‑150 (MQ=255)

                                                    cGATTATTTCACGACCTATGATTTTTACGATTTTGGACATTCGGTTTACCTCAAGTCACTATTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGa                                                                                                   <  1:1980406/150‑1 (MQ=255)

                                                     gctgatTTCACGACCGATGAGTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTCCGCCTGACGTTTCTGGaa                                                                                                  <  2:3811800/148‑1 (MQ=255)

                                                     gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                                  <  1:3363912/150‑1 (MQ=255)

                                                     gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                                  >  1:1471105/1‑150 (MQ=255)

                                                     gatgatTACACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCGCTAACTTTTTTACTTACTTCGCCTGACGGTTCTGGaa                                                                                                  >  1:1665729/1‑150 (MQ=255)

                                                      atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAAc                                                                                                 >  1:3253130/1‑150 (MQ=255)

                                                          tttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGc                                                                                             >  2:948657/1‑150 (MQ=255)

                                                           ttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGATCTCGCt                                                                                            >  2:397276/1‑150 (MQ=255)

                                                               cGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACAAGTTAAACAGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGATGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACACGCAggcg                                                                                        >  1:3005535/1‑150 (MQ=255)

                                                                gACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcgg                                                                                       >  1:3138267/1‑150 (MQ=255)

                                                                 aCCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcggc                                                                                      >  2:2660364/1‑150 (MQ=255)

                                                                    gatgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                                   >  1:2714803/1‑150 (MQ=255)

                                                                      tgatTTTTTCGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                                                                                 <  2:3367328/150‑1 (MQ=255)

                                                                          ttttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                             <  1:157667/150‑1 (MQ=255)

                                                                           tttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAg                                                                            >  2:3611950/1‑150 (MQ=255)

                                                                             tACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGcc                                                                          >  2:1095169/1‑150 (MQ=255)

                                                                             tACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGcc                                                                          >  2:232682/1‑150 (MQ=255)

                                                                               cGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTg                                                                        >  1:4258101/1‑150 (MQ=255)

                                                                                gATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGc                                                                       <  1:3852497/150‑1 (MQ=255)

                                                                                  ttttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaa                                                                     <  1:1350513/150‑1 (MQ=255)

                                                                                   tttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACACCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                                                                    <  1:760846/150‑1 (MQ=255)

                                                                                   tttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                                                                    <  2:1117517/150‑1 (MQ=255)

                                                                                   tttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                                                                    <  2:483960/150‑1 (MQ=255)

                                                                                   tttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                                                                    <  1:2749876/150‑1 (MQ=255)

                                                                                         cATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCCAACAGCgg                                                              >  2:3400218/1‑150 (MQ=255)

                                                                                          aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCggg                                                             <  1:618802/150‑1 (MQ=255)

                                                                                             aGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGa                                                          <  1:2031582/150‑1 (MQ=255)

                                                                                              ggTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc                                                         >  1:395716/1‑150 (MQ=255)

                                                                                              ggTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc                                                         <  2:2144881/150‑1 (MQ=255)

                                                                                                ttttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa                                                       >  2:2830270/1‑150 (MQ=255)

                                                                                                 tttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTCCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCAt                                                      <  2:840598/150‑1 (MQ=255)

                                                                                                   tCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCa                                                    >  1:3950083/1‑150 (MQ=255)

                                                                                                   tCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCa                                                    <  2:667934/150‑1 (MQ=255)

                                                                                                      tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGt                                                 >  2:251883/1‑150 (MQ=255)

                                                                                                       cAAGTCACTAGTTAAACTGAAACTCCCGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGTGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTg                                                <  1:1813451/150‑1 (MQ=255)

                                                                                                       cAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGGACCCAGGGGACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCCCGTg                                                <  2:1768450/150‑1 (MQ=255)

                                                                                                         aGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCCCTAACTTTTTTACTTACTTCGCCAGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGa                                              <  2:82149/150‑1 (MQ=255)

                                                                                                          gTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAg                                             >  1:2451695/1‑150 (MQ=255)

                                                                                                            cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTa                                           <  2:519226/150‑1 (MQ=255)

                                                                                                            cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTa                                           <  1:3772535/150‑1 (MQ=255)

                                                                                                             aCTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTATACTTACTTCGCCTGACGTTTCTGGAAATCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAg                                          >  1:2137718/1‑150 (MQ=255)

                                                                                                               tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGaa                                        <  2:3522482/150‑1 (MQ=255)

                                                                                                               tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGaa                                        <  2:166069/150‑1 (MQ=255)

                                                                                                                  ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTa                                     <  2:1969429/150‑1 (MQ=255)

                                                                                                                   tAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTaa                                    <  1:1057500/150‑1 (MQ=255)

                                                                                                                    aaaCTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTaaa                                   <  2:2739319/150‑1 (MQ=255)

                                                                                                                       cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTc                                >  1:2747833/1‑150 (MQ=255)

                                                                                                                         gAAACTCCAGACAAACAACGAGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCTCAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCg                              >  2:3418359/1‑150 (MQ=255)

                                                                                                                           aaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGa                            <  2:1524520/150‑1 (MQ=255)

                                                                                                                              tCCAGACAACCACCGCGGACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGAggg                         <  1:3604401/150‑1 (MQ=255)

                                                                                                                              tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAACGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATgg                         <  1:1102042/150‑1 (MQ=255)

                                                                                                                              tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATgg                         >  2:468349/1‑150 (MQ=255)

                                                                                                                               ccAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTTCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGa                        >  2:37267/1‑150 (MQ=255)

                                                                                                                                  gACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACt                     >  2:1952334/1‑150 (MQ=255)

                                                                                                                                  gACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACt                     <  1:351671/150‑1 (MQ=255)

                                                                                                                                   acaaacaaCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTg                    <  1:1757231/150‑1 (MQ=255)

                                                                                                                                    gcgacaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTTAACTCCGGATGGAACTgg                   >  1:2332710/4‑150 (MQ=255)

                                                                                                                                    gcgacaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTgg                   >  2:403613/4‑150 (MQ=255)

                                                                                                                                        caacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCCAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAg               >  2:256804/1‑150 (MQ=255)

                                                                                                                                          acGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc             >  2:4213250/1‑150 (MQ=255)

                                                                                                                                          acGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc             >  2:921104/1‑150 (MQ=255)

                                                                                                                                             cgTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACg          <  2:4159213/150‑1 (MQ=255)

                                                                                                                                                  ccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCa     <  2:1366119/150‑1 (MQ=255)

                                                                                                                                                  ccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCa     >  1:1695116/1‑150 (MQ=255)

                                                                                                                                                   cAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCAc    >  2:2620593/1‑150 (MQ=255)

                                                                                                                                                    aGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGC‑AACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACgg  <  2:4148613/150‑1 (MQ=255)

AAACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGG  >  NC_000913/2907850‑2908148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`