BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	192,573:1	+G	7.8%	coding (719/726 nt)	pyrH →	uridylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	192,573	1	.	G	7.8%	373.2 / 28.4	90	coding (719/726 nt)	pyrH	uridylate kinase
Reads supporting (aligned to +/- strand): ref base . (37/46); new base G (4/3); total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

GAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGA  >  NC_000913/192425‑192718

gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  4:664684/149‑1 (MQ=255)

gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  3:178139/149‑1 (MQ=255)

gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  2:295675/149‑1 (MQ=255)

gAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGAGTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGGGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAc                                                                                                                                                    <  1:55990/149‑1 (MQ=255)

    tGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGaa                                                                                                                                               <  3:26823/149‑1 (MQ=255)

      cTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGaata                                                                                                                                             >  4:603263/1‑149 (MQ=255)

        ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑gggaatta                                                                                                                                           >  2:90053/1‑143 (MQ=255)

         gAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAAtt                                                                                                                                          <  1:358538/149‑1 (MQ=255)

           aaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTcc                                                                                                                                        <  1:177449/149‑1 (MQ=255)

              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                     <  3:721699/149‑1 (MQ=255)

              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                     <  4:306099/149‑1 (MQ=255)

              agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                     <  1:298221/149‑1 (MQ=255)

                agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    <  1:90053/149‑5 (MQ=255)

                agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                   >  4:602889/1‑149 (MQ=255)

                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                  <  3:215490/149‑1 (MQ=255)

                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                  <  4:560341/149‑1 (MQ=255)

                 gCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                  >  4:595108/1‑149 (MQ=255)

                  cTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATg                                                                                                                                 <  3:32818/149‑1 (MQ=255)

                   tGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    >  2:239223/1‑142 (MQ=255)

                   tGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                    <  1:239223/146‑5 (MQ=255)

                    gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGa                                                                                                                               <  1:27154/149‑1 (MQ=255)

                    gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGa                                                                                                                               >  2:137707/1‑149 (MQ=255)

                     aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGat                                                                                                                              <  3:53662/149‑1 (MQ=255)

                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                        <  3:153998/149‑1 (MQ=255)

                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                        <  2:31535/149‑1 (MQ=255)

                           aTGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                        <  3:603263/149‑1 (MQ=255)

                            tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAg                                                                                                                       <  3:481725/149‑1 (MQ=255)

                            tGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAg                                                                                                                       <  4:212431/149‑1 (MQ=255)

                             ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      >  3:506512/1‑149 (MQ=255)

                             ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      <  4:691622/149‑1 (MQ=255)

                              gACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTCTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAggg                                                                                                                     <  1:464012/149‑1 (MQ=255)

                                    gcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      >  4:424583/1‑142 (MQ=255)

                                    gcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                      <  3:424583/142‑1 (MQ=255)

                                        ccTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcg                                                                                                           >  1:108849/1‑149 (MQ=255)

                                         cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                          >  2:229759/1‑149 (MQ=255)

                                         cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTAAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                          >  3:466368/1‑149 (MQ=255)

                                          ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                         >  2:34101/1‑149 (MQ=255)

                                          ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                         >  4:417657/1‑149 (MQ=255)

                                                      cGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgg                                                                                             <  3:145195/149‑1 (MQ=255)

                                                       gTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                            <  1:94360/149‑1 (MQ=255)

                                                       gTGACCATAAATCACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggg                                                                                            <  3:701414/149‑1 (MQ=255)

                                                           ccATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCg                                                                                        >  3:76247/1‑149 (MQ=255)

                                                            cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                       <  1:719248/149‑1 (MQ=255)

                                                            cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                       >  3:615008/1‑149 (MQ=255)

                                                             aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTa                                                                                      >  1:581634/1‑149 (MQ=255)

                                                              tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                     <  4:299153/149‑1 (MQ=255)

                                                                   tACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                >  2:511125/1‑149 (MQ=255)

                                                                   tACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                <  2:547657/149‑1 (MQ=255)

                                                                    aCCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  3:261054/149‑1 (MQ=255)

                                                                    aCCGATTCGTGTTTTCAATATGAACAAACCCGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  4:112670/149‑1 (MQ=255)

                                                                      cGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGtt                                                                             <  2:10178/149‑1 (MQ=255)

                                                                         ttCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTtata                                                                          <  4:603469/149‑1 (MQ=255)

                                                                           cGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAATGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                        <  2:453912/149‑1 (MQ=255)

                                                                           cGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                        <  3:67920/149‑1 (MQ=255)

                                                                              gtTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCg                                                                     >  1:431305/1‑149 (MQ=255)

                                                                                tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGtt                                                                   <  1:455486/149‑1 (MQ=255)

                                                                                   cAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                >  3:185874/1‑149 (MQ=255)

                                                                                    aaTATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttg                                                               >  4:602070/1‑149 (MQ=255)

                                                                                     atatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgt                                                              <  3:602889/149‑1 (MQ=255)

                                                                                      tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  2:388115/131‑1 (MQ=255)

                                                                                      tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               >  1:388115/1‑131 (MQ=255)

                                                                                         gAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                          <  3:595108/149‑1 (MQ=255)

                                                                                          aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                >  2:248887/1‑142 (MQ=255)

                                                                                          aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                <  1:248887/142‑1 (MQ=255)

                                                                                           aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               >  2:362491/1‑126 (MQ=255)

                                                                                           aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                               <  1:362491/126‑1 (MQ=255)

                                                                                           aCAAACCGGGTGCGCTGCGCCGTGAGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcagg                                                        >  2:438251/1‑149 (MQ=255)

                                                                                             aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                      <  4:447453/149‑1 (MQ=255)

                                                                                               aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                        >  3:78842/1‑81 (MQ=255)

                                                                                               aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                        <  4:78842/81‑1 (MQ=255)

                                                                                                 cGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATGCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGttt                                                  <  2:374536/149‑1 (MQ=255)

                                                                                                        gctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACg                                           >  3:376674/1‑149 (MQ=255)

                                                                                                        gctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACg                                           >  4:752819/1‑149 (MQ=255)

                                                                                                              ccGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAgc                                     >  1:655747/1‑149 (MQ=255)

                                                                                                                       aaTGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCaa                            >  4:557385/1‑149 (MQ=255)

                                                                                                                            gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  4:193541/1‑149 (MQ=255)

                                                                                                                            gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       <  3:250912/149‑1 (MQ=255)

                                                                                                                            gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  4:634195/1‑149 (MQ=255)

                                                                                                                            gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                       >  2:605502/1‑149 (MQ=255)

                                                                                                                              gAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGc                     <  2:119615/149‑1 (MQ=255)

                                                                                                                                     aagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa              >  1:681667/1‑149 (MQ=255)

                                                                                                                                     aagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa              <  2:681667/149‑1 (MQ=255)

                                                                                                                                      agGGACTTTAATCAC‑GGTATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa             <  2:674203/149‑1 (MQ=255)

                                                                                                                                      agGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaaa             <  2:431305/149‑1 (MQ=255)

                                                                                                                                        ggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGt           >  4:288210/1‑149 (MQ=255)

                                                                                                                                            tttAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt       >  3:371206/1‑149 (MQ=255)

                                                                                                                                                aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTa                  <  2:293961/134‑1 (MQ=255)

                                                                                                                                                aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTa                  >  1:293961/1‑134 (MQ=255)

                                                                                                                                                aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   >  2:13757/1‑149 (MQ=255)

                                                                                                                                                aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   >  4:424874/1‑149 (MQ=255)

                                                                                                                                                aTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg   <  1:591076/149‑1 (MQ=255)

                                                                                                                                                 tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  1:724220/1‑149 (MQ=255)

                                                                                                                                                 tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  3:270009/1‑149 (MQ=255)

                                                                                                                                                 tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCCGAGCGTACAAAAGTACGTGAGGa  >  3:262453/1‑149 (MQ=255)

                                                                                                                                                 tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATGCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa  >  4:439553/1‑149 (MQ=255)

GAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGA  >  NC_000913/192425‑192718

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