BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	192,579	A→T	7.7%	*242L (TAA→TTA)	pyrH →	uridylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	192,579	0	A	T	7.7%	309.0 / 13.5	91	*242L (TAA→TTA)	pyrH	uridylate kinase
Reads supporting (aligned to +/- strand): ref base A (38/46); new base T (4/3); total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

CTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCG  >  NC_000913/192431‑192723

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

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

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

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

        agagCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGt                                                                                                                                          <  3:721699/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                                                                                                                                       >  4:595108/1‑149 (MQ=255)

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

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

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

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

             tGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtga                                                                                                                                         >  2:239223/1‑142 (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                                                                                                                             <  3:603263/149‑1 (MQ=255)

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

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

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

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

                       ggACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                           >  3:506512/1‑149 (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                                                                                                              >  4:417657/1‑149 (MQ=255)

                                    ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                              >  2:34101/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                                                                                            >  3:615008/1‑149 (MQ=255)

                                                      cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                            <  1:719248/149‑1 (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                                                                                    >  1:388115/1‑131 (MQ=255)

                                                                                tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                                    <  2:388115/131‑1 (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                                                                                    <  1:362491/126‑1 (MQ=255)

                                                                                     aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                                    >  2:362491/1‑126 (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                            <  3:250912/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                               aagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa                   >  1:681667/1‑149 (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       >  3:270009/1‑149 (MQ=255)

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

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

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

                                                                                                                                                ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGc   >  1:622309/1‑149 (MQ=255)

                                                                                                                                                ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGc   <  2:622309/149‑1 (MQ=255)

                                                                                                                                                 gAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCg  >  1:300965/1‑149 (MQ=255)

CTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCG  >  NC_000913/192431‑192723

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