BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,935	C→A	15.5%	K96N (AAG→AAT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,935	0	C	A	15.5%	342.0 / 40.6	109	K96N (AAG→AAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (43/49); new base A (9/8); total (52/57)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 110 total aligned reads displayed.

ACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTT  >  NC_000913/4116793‑4117083

acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCg                                                                                                                                                >  1:588639/1‑149 (MQ=255)

  ttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGt                                                                                                                                              >  3:448650/1‑149 (MQ=255)

  ttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGt                                                                                                                                              >  3:305045/1‑149 (MQ=255)

  ttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGAGATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGt                                                                                                                                              >  2:803513/1‑149 (MQ=255)

   tcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                                                                                                             <  3:775125/149‑1 (MQ=255)

    cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGTCAATGGCGTTATAGATAGGATTGCCGGttt                                                                                                                                            <  2:737911/149‑1 (MQ=255)

    cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGttt                                                                                                                                            <  1:103188/149‑1 (MQ=255)

    cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGttt                                                                                                                                            <  3:174213/149‑1 (MQ=255)

     acttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTc                                                                                                                                           <  4:131481/149‑1 (MQ=255)

       tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCtt                                                                                                                                         >  2:334948/1‑149 (MQ=255)

       tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCtt                                                                                                                                         <  4:182256/149‑1 (MQ=255)

        ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttt                                                                                                                                        <  1:339100/149‑1 (MQ=255)

        ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttt                                                                                                                                        <  4:560164/149‑1 (MQ=255)

         tGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCtttt                                                                                                                                       <  2:616072/149‑1 (MQ=255)

         tGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCtttt                                                                                                                                       <  4:85536/149‑1 (MQ=255)

          ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttttt                                                                                                                                      <  2:620917/149‑1 (MQ=255)

                aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGa                                                                                                                                >  3:459877/1‑149 (MQ=255)

                   aaaGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACaa                                                                                                                             <  2:593768/149‑1 (MQ=255)

                    aaGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAAt                                                                                                                            <  2:933747/149‑1 (MQ=255)

                     aGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATa                                                                                                                           <  3:211525/149‑1 (MQ=255)

                        aCGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTg                                                                                                                        <  4:843329/149‑1 (MQ=255)

                           ggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGtt                                                                                                                     >  2:893494/1‑149 (MQ=255)

                               aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                                 <  2:250403/149‑1 (MQ=255)

                               aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                                 >  1:459943/1‑149 (MQ=255)

                                aTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACg                                                                                                                >  3:900524/1‑149 (MQ=255)

                                   aCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTg                                                                                                             <  2:462494/149‑1 (MQ=255)

                                   aCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTCCTTTTTCCCAGACAATAGTGGTTTCACGCTg                                                                                                             >  4:75716/1‑149 (MQ=255)

                                     cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTATGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                                                                                                           >  2:412140/1‑149 (MQ=255)

                                       gACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                                                                                                         <  2:145477/149‑1 (MQ=255)

                                          cGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                                                                                                      <  3:347148/149‑1 (MQ=255)

                                           ggTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTaa                                                                                                     <  3:900266/149‑1 (MQ=255)

                                               ttGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAcc                                                                                                 <  4:417005/149‑1 (MQ=255)

                                                 gctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGa                                                                                               <  1:391764/149‑1 (MQ=255)

                                                     cGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                      >  3:142176/1‑122 (MQ=255)

                                                     cGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                      <  4:142176/122‑1 (MQ=255)

                                                              tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                  >  1:81779/1‑149 (MQ=255)

                                                              tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                  >  3:751078/1‑149 (MQ=255)

                                                              tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                  >  1:382729/1‑149 (MQ=255)

                                                               cttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGa                                                                                 >  1:748817/1‑149 (MQ=255)

                                                                ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                >  1:533095/1‑149 (MQ=255)

                                                                ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                >  1:166565/1‑149 (MQ=255)

                                                                ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                <  2:771939/149‑1 (MQ=255)

                                                                 tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                               <  3:825610/149‑1 (MQ=255)

                                                                   tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCgg                                                                             >  4:721906/1‑149 (MQ=255)

                                                                     aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGACAGCTGCAATTTGATCGGaa                                                                           <  2:179846/149‑1 (MQ=255)

                                                                          tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                      >  1:76181/1‑149 (MQ=255)

                                                                          tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                      >  2:453270/1‑149 (MQ=255)

                                                                           cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAGTTGATCGGAACTGata                                                                     <  2:464347/149‑1 (MQ=255)

                                                                            aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGatat                                                                    >  2:439735/1‑149 (MQ=255)

                                                                            aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGatat                                                                    <  4:751078/149‑1 (MQ=255)

                                                                            aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGatat                                                                    >  1:832209/1‑149 (MQ=255)

                                                                              gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                  >  2:66062/1‑149 (MQ=255)

                                                                               ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCTGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                 >  4:307259/1‑149 (MQ=255)

                                                                               ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGCCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                 >  3:568642/1‑149 (MQ=255)

                                                                                  tAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                              <  1:767792/149‑1 (MQ=255)

                                                                                    aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                            >  4:47035/1‑149 (MQ=255)

                                                                                        cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                        <  3:710070/149‑1 (MQ=255)

                                                                                         tCGCNGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAg                                                       >  3:306442/1‑149 (MQ=255)

                                                                                         tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                >  4:735136/1‑124 (MQ=255)

                                                                                         tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                <  3:735136/124‑1 (MQ=255)

                                                                                           gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                     >  1:749911/1‑149 (MQ=255)

                                                                                               aTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTc                                                 <  4:305045/149‑1 (MQ=255)

                                                                                                  tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                              <  2:533095/149‑1 (MQ=255)

                                                                                                     gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAg                                           >  4:145008/1‑149 (MQ=255)

                                                                                                         tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                       <  4:306442/149‑1 (MQ=255)

                                                                                                           cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                     <  1:439735/149‑1 (MQ=255)

                                                                                                           cgacgGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                     <  2:382729/149‑1 (MQ=255)

                                                                                                            gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                    >  4:703087/1‑149 (MQ=255)

                                                                                                                gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                <  4:448650/149‑1 (MQ=255)

                                                                                                                 cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                               <  1:289829/149‑1 (MQ=255)

                                                                                                                      ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                <  1:106191/127‑1 (MQ=255)

                                                                                                                      ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                >  2:106191/1‑127 (MQ=255)

                                                                                                                      ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                          <  3:307259/149‑1 (MQ=255)

                                                                                                                      ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                          >  4:614229/1‑149 (MQ=255)

                                                                                                                      ccAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                          >  2:19576/1‑149 (MQ=255)

                                                                                                                       cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                         <  3:436084/149‑1 (MQ=255)

                                                                                                                       cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                         >  4:49249/1‑149 (MQ=255)

                                                                                                                       cAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                               >  3:339220/1‑95 (MQ=255)

                                                                                                                       cAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                               <  4:339220/95‑1 (MQ=255)

                                                                                                                        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                        >  3:364999/1‑149 (MQ=255)

                                                                                                                        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                         <  3:684009/148‑1 (MQ=255)

                                                                                                                        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                         >  4:684009/1‑148 (MQ=255)

                                                                                                                         gACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACTATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGGTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcca                       >  2:393344/1‑148 (MQ=255)

                                                                                                                          aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccaa                      <  3:361787/149‑3 (MQ=255)

                                                                                                                            aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGa                    <  1:334948/149‑1 (MQ=255)

                                                                                                                            aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGa                    >  4:315630/1‑149 (MQ=255)

                                                                                                                                 cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTcc               <  2:788122/149‑1 (MQ=255)

                                                                                                                                 cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTcc               <  3:315630/149‑1 (MQ=255)

                                                                                                                                  gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCa              <  2:81779/149‑1 (MQ=255)

                                                                                                                                   ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCAt             >  3:770309/1‑149 (MQ=255)

                                                                                                                                    tataGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCAtt            <  4:900524/149‑1 (MQ=255)

                                                                                                                                     atagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTg           <  1:738258/149‑1 (MQ=255)

                                                                                                                                      tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATAgg          <  3:703087/149‑1 (MQ=255)

                                                                                                                                       agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTggg         >  2:322573/1‑149 (MQ=255)

                                                                                                                                       agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTggg         >  3:865561/1‑149 (MQ=255)

                                                                                                                                        gatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGt        >  2:535283/1‑149 (MQ=255)

                                                                                                                                         atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTc       <  3:47035/149‑1 (MQ=255)

                                                                                                                                            ggCTTGCCGGTTTC‑TTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCgtgt   <  2:892405/149‑1 (MQ=255)

                                                                                                                                             gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCgtgt   <  1:893494/149‑1 (MQ=255)

                                                                                                                                              cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGtt  <  2:78082/149‑1 (MQ=255)

ACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTT  >  NC_000913/4116793‑4117083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`