BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,757,091	C→A	7.6%	A465E (GCA→GAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,091	0	C	A	7.6%	1034.8 / 52.2	290	A465E (GCA→GAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base C (127/141); new base A (10/12); total (137/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 293 total aligned reads displayed.

TGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAAT  >  NC_000913/1756942‑1757240

gaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACACCGc                                                                                                                                                       <  1:5761189‑M1/10‑1 (MQ=255)

   gCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGAGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                                                                                                                                    >  2:1985380/1‑150 (MQ=255)

      gCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGt                                                                                                                                                 >  2:3790890/1‑150 (MQ=255)

      gCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGt                                                                                                                                                 <  1:5911045/150‑1 (MQ=255)

       cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGtt                                                                                                                                                <  1:6218320/150‑1 (MQ=255)

        aGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                                                                                                                               <  1:2944060/150‑1 (MQ=255)

          tttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCAc                                                                                                                                             >  2:5630146‑M1/131‑150 (MQ=255)

          cTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCAc                                                                                                                                             >  1:3926471/1‑150 (MQ=255)

          cTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCAc                                                                                                                                             >  2:5903412/1‑150 (MQ=255)

              ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTcc                                                                                                                                         >  2:5811805/1‑150 (MQ=255)

              ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTcc                                                                                                                                         >  1:5814134/1‑150 (MQ=255)

               tAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCt                                                                                                                                        >  1:2143076/1‑150 (MQ=255)

                 aaGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGt                                                                                                                                      >  2:3947649/1‑150 (MQ=255)

                        cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTg                                                                                                                               <  1:3633656/150‑1 (MQ=255)

                          cTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCt                                                                                                                             <  1:375290/150‑1 (MQ=255)

                           ttCTACTGATTATTTCTACCGTCTGGATAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCtt                                                                                                                            >  1:2465655/1‑150 (MQ=255)

                           ttCTACTGATTATTTCTACCGTCTGGATAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCtt                                                                                                                            >  1:51719/1‑150 (MQ=255)

                           ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCtt                                                                                                                            >  2:247831/1‑150 (MQ=255)

                                 tgatgaTTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGa                                                                                                                      <  2:3337190/150‑1 (MQ=255)

                                  ctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa                                                                                                                     >  1:5379370‑M1/107‑150 (MQ=255)

                                    tgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatt                                                                                                                   >  2:1000756/1‑150 (MQ=255)

                                        ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaatt                                                                                                               <  1:3790890/150‑1 (MQ=255)

                                         tCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAttt                                                                                                              <  1:4763758/150‑1 (MQ=255)

                                           tACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                                                            >  2:4475515/1‑150 (MQ=255)

                                           tACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                                                            >  1:2235101/1‑150 (MQ=255)

                                                  tGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGa                                                                                                     >  1:1823598/1‑150 (MQ=255)

                                                   gggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGaa                                                                                                    <  2:3269158/150‑1 (MQ=255)

                                                    ggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAg                                                                                                   >  2:187925/1‑150 (MQ=255)

                                                     gTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAgc                                                                                                  >  2:1661436/1‑150 (MQ=255)

                                                       aaaGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAgcgc                                                                                                <  1:2605296/150‑1 (MQ=255)

                                                         aGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                                                                              <  1:4369636/150‑1 (MQ=255)

                                                         aGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                                                                              <  2:2734921/150‑1 (MQ=255)

                                                             cTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGAt                                                                                          >  2:6251793/1‑150 (MQ=255)

                                                               ggCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATgg                                                                                        <  1:1483309/150‑1 (MQ=255)

                                                               ggCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTCTTCACGTCCTGTAAAAAAGCTTTTGTGTAATAATATCAATATCGATGGGCCCTGAtcg                                                                                        >  2:4260834/1‑148 (MQ=255)

                                                                gCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATggg                                                                                       <  1:4630710/150‑1 (MQ=255)

                                                                gCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATggg                                                                                       >  2:1426649/1‑150 (MQ=255)

                                                                 ctctGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgc                                                                                      >  2:1803158/1‑150 (MQ=255)

                                                                    tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgtacgcacctaccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCAGATGGGCGCt                                                                                   <  1:4847567‑M1/78‑1 (MQ=255)

                                                                          gCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAGCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttttt                                                                             <  2:2836914/150‑1 (MQ=255)

                                                                          gCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttttt                                                                             <  1:6232840/150‑1 (MQ=255)

                                                                           cGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTttta                                                                            >  2:4930842/1‑150 (MQ=255)

                                                                           cGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTttta                                                                            <  2:2279767/150‑1 (MQ=255)

                                                                            ggTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttat                                                                           <  2:6169839/150‑1 (MQ=255)

                                                                              tCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattt                                                                         >  1:230239/1‑150 (MQ=255)

                                                                                 ggCACACAACGGTGTCGTTGTTGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                                                                      <  1:3097087/150‑1 (MQ=255)

                                                                                 ggCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                                                                      <  1:187925/150‑1 (MQ=255)

                                                                                 ggCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                                                                      <  1:6305795/150‑1 (MQ=255)

                                                                                 ggCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                                                                      <  1:1299538/150‑1 (MQ=255)

                                                                                 ggCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGGACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATAGTGCTTTTGTGAATTAATTTGTATATCTAAGCGCCCTTATGGGCGCTTTTTTTATTTAAt                                                                      <  1:5894919/150‑1 (MQ=255)

                                                                                    acacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGACTTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                                   >  2:6424817/1‑150 (MQ=255)

                                                                                    acacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                                   >  2:5431188/1‑150 (MQ=255)

                                                                                    acacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                                   >  1:1272112/1‑150 (MQ=255)

                                                                                    acacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                                   >  1:240990/1‑150 (MQ=255)

                                                                                             tGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGaa                                                          <  1:680914/150‑1 (MQ=255)

                                                                                             tGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGaa                                                          <  1:247831/150‑1 (MQ=255)

                                                                                              gACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAg                                                         >  2:6288310/1‑150 (MQ=255)

                                                                                                 gTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCaa                                                      >  2:2195243/1‑150 (MQ=255)

                                                                                                  ttGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAAt                                                     <  2:3830349/150‑1 (MQ=255)

                                                                                                    acaaacatgatctcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaa                                                   >  2:2241462‑M1/41‑150 (MQ=255)

                                                                                                     tAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaaa                                                  <  1:6074641/150‑1 (MQ=255)

                                                                                                      aaacatgatctcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaaaa                                                 <  2:4148253‑M1/112‑1 (MQ=255)

                                                                                                       gTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaaaaa                                                <  2:5982283/150‑1 (MQ=255)

                                                                                                       aacatgatctcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTACATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaaaaa                                                >  1:3128255‑M1/38‑150 (MQ=255)

                                                                                                        ttATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATaaaaaa                                               <  2:5067664/150‑1 (MQ=255)

                                                                                                              tttCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaaat                                         >  2:5135096/1‑150 (MQ=255)

                                                                                                                tcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCg                                       <  2:1537781‑M1/122‑1 (MQ=255)

                                                                                                                tcatatcagggacttgttcgcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCg                                       <  2:2396622‑M1/122‑1 (MQ=255)

                                                                                                                  tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                     <  1:1688330/150‑1 (MQ=255)

                                                                                                                  tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                     >  2:6136295/1‑150 (MQ=255)

                                                                                                                  tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGAGGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                     <  1:4087097/150‑1 (MQ=255)

                                                                                                                       cACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                                <  1:6404623/150‑1 (MQ=255)

                                                                                                                       cACTGGTACCGAGCGGCACTACTAACACCGAATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                                <  1:1115326/150‑1 (MQ=255)

                                                                                                                        aCTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCAc                               >  2:3546726/1‑150 (MQ=255)

                                                                                                                         cTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACt                              <  1:2569371/150‑1 (MQ=255)

                                                                                                                         cTGCTACCGAGCGGCACTACTAACACCGCAACTGTTCACGTCCTGTAATATAGCTTTACTGAATTAATTTGAATAACCAAGCGCCCTGATGGGCGCTTTTTTTAATTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACt                              <  2:2891160/150‑1 (MQ=255)

                                                                                                                              aCCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTAGATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGAGTTCACtatata                         >  2:3457432/1‑150 (MQ=255)

                                                                                                                                 gAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATAGGATTTCACTAAATAAtc                      >  1:5363099/1‑150 (MQ=255)

                                                                                                                                   gcaccttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAGTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAAtctc                    <  2:2375435‑M1/141‑1 (MQ=255)

                                                                                                                                    cGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCa                   <  2:6575602/150‑1 (MQ=255)

                                                                                                                                      gCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACt                 >  2:3517/1‑150 (MQ=255)

                                                                                                                                        ttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAAACGGATTTCACTATATAATCTCACttt               <  2:4870556‑M1/146‑1 (MQ=255)

                                                                                                                                        ttccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATAGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACttt               <  2:2540910‑M1/146‑1 (MQ=255)

                                                                                                                                        actactAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACttt               >  2:5350628/1‑150 (MQ=255)

                                                                                                                                         ctactaACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTa              <  1:530697/150‑1 (MQ=255)

                                                                                                                                          tactaACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTAt             <  2:1900469/150‑1 (MQ=255)

                                                                                                                                          tactaACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTAt             <  2:5030596/150‑1 (MQ=255)

                                                                                                                                          ccCTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACt                              >  1:845900‑M1/3‑133 (MQ=255)

                                                                                                                                            ctaACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAATAATCAAATCGGATTTCACTATATAATCTCTCTTTATCt           >  2:4371698/1‑150 (MQ=255)

                                                                                                                                             taACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTa          >  2:3169455/1‑150 (MQ=255)

                                                                                                                                             taACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATAGGATTTAACTATATAATCTCACTTTATCTa          >  1:2374187/1‑150 (MQ=255)

                                                                                                                                              aaCCCCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTaa         <  2:5284502/150‑1 (MQ=255)

                                                                                                                                              aaCACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAATATCAAATCGGATTTCACTATATAATCTCACTTTATCTaa         <  1:5106500/150‑1 (MQ=255)

                                                                                                                                              aaCACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTaa         <  2:5558301/150‑1 (MQ=255)

                                                                                                                                              aaCACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTaa         <  1:3468061/150‑1 (MQ=255)

                                                                                                                                              aaCACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATAGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATATCACTTTATCTaa         >  1:3978628/1‑150 (MQ=255)

                                                                                                                                                cacCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGa       <  2:71818/150‑1 (MQ=255)

                                                                                                                                                  ccGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATAAAATCGGATTTCACTATATAATCTAACTTTATCTAAGATg     >  1:5554436/1‑150 (MQ=255)

                                                                                                                                                     cATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAAt  <  2:1089588/150‑1 (MQ=255)

                                                                                                                                                     cATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTAACTATATAATCTCACTTTATCTAAGATGaaa  >  1:4981971/1‑149 (MQ=255)

TGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAAT  >  NC_000913/1756942‑1757240

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