BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,907,998	A→T	5.3%	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.3%	1157.0 / 31.6	346	intergenic (‑57/+31)	eno/pyrG	enolase/CTP synthetase
Reads supporting (aligned to +/- strand): ref base A (170/155); new base T (5/15); total (175/171)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

GAAACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACG  >  NC_000913/2907849‑2908147

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

     caccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                                  >  1:4285856/1‑150 (MQ=255)

     caccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGGCGATGATTTCACGACCGATGATTGTTACGATTTTGGACCGTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCACGGtc                                                                                                                                                  >  2:552019/1‑149 (MQ=255)

        ccTCCAGATGTACTTCGGCTTCAACAGTCGGTTTACCTGGGGAGTCGATGATTTCACGACCGTTGATTTTTACGATTTTGTACATTAGGTTTTCCTCATTTCTCTTGTTAAACTGAAACTCCAGTCAAACAACGCGTACCCAGGGTAcgc                                                                                                                                               <  1:5443059/150‑1 (MQ=255)

        ccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                               <  1:1556132/150‑1 (MQ=255)

        ccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                               <  2:5764004/150‑1 (MQ=255)

        ccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                               <  2:2885880/150‑1 (MQ=255)

        ccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                               <  1:795885/150‑1 (MQ=255)

          tCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGt                                                                                                                                             <  1:424387/150‑1 (MQ=255)

           ccAGATGTCCTTCGGCTTCACCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTACACTGAAACTCCAGACAAACACCGCGTACCCAGGGTACGCGtt                                                                                                                                            <  1:552019/150‑1 (MQ=255)

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

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

            cAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGCCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTg                                                                                                                                           <  2:4656953/150‑1 (MQ=255)

              gATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGcc                                                                                                                                         <  2:504707/150‑1 (MQ=255)

                        ggCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCGCTAACtt                                                                                                                               >  2:728993/1‑150 (MQ=255)

                          cTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttt                                                                                                                             <  2:3073636/150‑1 (MQ=255)

                           ttCAACTGGCGGGTTACCAGGGTCGTCTTTTATTTCACGACCGTTTTTTTTTACTATTTTGTACATTAGGTTTTCCTCAAGTCACTTTTTAAACTGAAACTCCAGAAAAACAACGCGTACCCTGGGTACGCGTTGCCGCTCTAACttttt                                                                                                                            <  1:3890240/150‑1 (MQ=255)

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

                            tCAACAGTCGGGTTTCCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACGCAGGGTACGCGTTGCCGCTCTAACtttttt                                                                                                                           <  2:5348462/150‑1 (MQ=255)

                             cAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                                                                          <  2:3980595/150‑1 (MQ=255)

                              aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTAAGATTTTGGACATTAGGTTTTCCTCAAGTCACAAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttttttc                                                                                                                         >  1:5103778/1‑148 (MQ=255)

                               aCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttact                                                                                                                        >  2:2114187/1‑150 (MQ=255)

                                cAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCGCTAACTTTTttactt                                                                                                                       >  2:4225519/1‑150 (MQ=255)

                                 aGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                      >  2:418262/1‑150 (MQ=255)

                                 aGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                      >  1:6473585/1‑150 (MQ=255)

                                 aGTCGGGTTACCACGGGAGTCGATGATTTAACGACCGATGATTTTTACGATTTTTGACATTAGGTTTTCCTCAAGTCACAAGTTAAAAAGAAACTACAGACAAACAAAGCGTACCCACGGTAAGCCTTGCCGCAATAAAATTTttaatta                                                                                                                      >  1:1915866/1‑150 (MQ=255)

                                  gTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttacttac                                                                                                                     >  2:2208106/1‑150 (MQ=255)

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

                                         tACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCt                                                                                                              <  2:4789314/150‑1 (MQ=255)

                                           ccAGGGGTGTCGTGGTTTTCTCGCCCGTTGTTTTTTACGTTTTTGGACTTTTGGTTTTCCTCTAGTCAGTAGTTAAACTGAAACTCCAGACAAACACCGCGTCCCCTGGGTACGGGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGa                                                                                                            <  1:3460982/150‑1 (MQ=255)

                                                 gAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTc                                                                                                      <  2:1825553/150‑1 (MQ=255)

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

                                                   gTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTg                                                                                                    <  1:4199605/150‑1 (MQ=255)

                                                   gTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTg                                                                                                    <  1:3773152/150‑1 (MQ=255)

                                                    tCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTgg                                                                                                   >  2:2785116/1‑150 (MQ=255)

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

                                                       atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAAAGCGTACCCAGGGTACGCGTTGCCGGTCTTACTTTTTTACTTACGTCGCCTGACGTGTCTGGAAc                                                                                                >  2:4621100/1‑150 (MQ=255)

                                                        tgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACt                                                                                               >  2:301494/1‑150 (MQ=255)

                                                          atTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCg                                                                                             >  2:11984/1‑150 (MQ=255)

                                                              cACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTgg                                                                                         >  1:5821015/1‑150 (MQ=255)

                                                              cACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTgg                                                                                         <  1:102299/150‑1 (MQ=255)

                                                                 gACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGGACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcgg                                                                                      <  1:2116756/150‑1 (MQ=255)

                                                                  aCCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGGCGTTTCTGGAACTCGCTggcggc                                                                                     <  1:6359594/150‑1 (MQ=255)

                                                                  aCCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcggc                                                                                     <  1:5059903/150‑1 (MQ=255)

                                                                   ccGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCt                                                                                    >  2:1572657/1‑150 (MQ=255)

                                                                    cGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCtt                                                                                   <  2:105347/150‑1 (MQ=255)

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

                                                                      atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                                 >  1:455928/1‑150 (MQ=255)

                                                                      atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                                 >  2:4321799/1‑150 (MQ=255)

                                                                      atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                                 >  2:6083416/1‑150 (MQ=255)

                                                                       tgatTTTTACGATTTTGGACATTCGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                                                                                <  1:1987541/150‑1 (MQ=255)

                                                                          tttttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaa                                                                             <  1:1335369/150‑1 (MQ=255)

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

                                                                               aCGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCt                                                                        <  1:2535240/150‑1 (MQ=255)

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

                                                                                 gATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAggg                                                                                                                                                    >  2:1220649/1‑72 (MQ=255)

                                                                                 gATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGc                                                                      >  1:5880651/1‑150 (MQ=255)

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

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

                                                                                         aCATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCg                                                              <  2:3813909/150‑1 (MQ=255)

                                                                                           aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTATTTTACTTACTTCGCCTGACGTATCTCGAACTCGCTCCCGGCTTTCACAAAGCCTGCACACAGCggg                                                            >  2:6233946/1‑150 (MQ=255)

                                                                                            ttAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGt                                                           >  2:1796155/1‑150 (MQ=255)

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

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

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

                                                                                                  tttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTgg                                                                                         >  2:1989909/1‑114 (MQ=255)

                                                                                                       tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGt                                                <  2:1022058/150‑1 (MQ=255)

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

                                                                                                       tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGACTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGt                                                >  1:4478619/1‑150 (MQ=255)

                                                                                                        cAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGTCGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTg                                               <  2:4890392/150‑1 (MQ=255)

                                                                                                          aGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGa                                             <  1:479005/150‑1 (MQ=255)

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

                                                                                                            tCACTAGTTAAACTGAAACTCCAGACAAGCAACGCGTACCCTGGGTACGCGt                                                                                                                                             <  2:36297/52‑1 (MQ=37)

                                                                                                            tCACTAGTTAAACTGAAACTCCAGACAAGCAACGCGTACCCTGGGTACGCGt                                                                                                                                             >  1:36297/1‑52 (MQ=37)

                                                                                                             cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                                                                                >  2:56376/1‑112 (MQ=255)

                                                                                                             cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                                                                                <  1:56376/112‑1 (MQ=255)

                                                                                                              aCTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAg                                         >  1:3004600/1‑150 (MQ=255)

                                                                                                               cTAGTTAAACTGAAACTCCAGACAAACCACGCGTACCCAGGGTACGCGTAGCCGCTCTAACTATTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACTATCACGTGGAGTAGa                                        <  2:6523193/150‑1 (MQ=255)

                                                                                                               cTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGa                                        <  1:3717366/150‑1 (MQ=255)

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

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

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

                                                                                                                    tAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTaa                                   >  1:2855251/1‑150 (MQ=255)

                                                                                                                     aaaCTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGaa                                       >  1:2224009/1‑145 (MQ=255)

                                                                                                                        cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAGACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTc                               <  2:599581/150‑1 (MQ=255)

                                                                                                                          gAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCg                             <  2:6566378/150‑1 (MQ=255)

                                                                                                                          gAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCg                             >  2:5165645/1‑150 (MQ=255)

                                                                                                                             aCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGAt                          <  2:3282520/150‑1 (MQ=255)

                                                                                                                               tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAGGTAAACTCCGGATgg                        >  1:6689216/1‑150 (MQ=255)

                                                                                                                                   gACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTGTGAAAGccgc                                                                       >  1:1801646/1‑97 (MQ=255)

                                                                                                                                   gACAAACAACGCGGACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACt                    <  2:901884/150‑1 (MQ=255)

                                                                                                                                          aacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc             >  2:1447897/1‑150 (MQ=255)

                                                                                                                                          aacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCCTCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCACGc             >  1:2672334/1‑150 (MQ=255)

                                                                                                                                           acGCGTACCCTGGGTACGCGTTGTCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTTACAAAGc                                                                          <  1:3227316/88‑1 (MQ=39)

                                                                                                                                               gTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGa        <  2:3521610/150‑1 (MQ=255)

                                                                                                                                               gTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGa        >  2:5447682/1‑150 (MQ=255)

                                                                                                                                                   ccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGc                                                                          >  2:6461379/1‑80 (MQ=255)

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

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

                                                                                                                                                     aGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACg  >  2:3919902/1‑150 (MQ=255)

GAAACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACG  >  NC_000913/2907849‑2908147

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