BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	192,579	A→T	5.0%	*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	5.0%	402.7 / 7.7	114	*242L (TAA→TTA)	pyrH	uridylate kinase
Reads supporting (aligned to +/- strand): ref base A (60/48); new base T (2/4); total (62/52)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 117 total aligned reads displayed.

GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCG  >  NC_000913/192433‑192723

ggAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGaataa                                                                                                                                                 >  1:56749/1‑148 (MQ=255)

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

   aaaaGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTcc                                                                                                                                             >  4:758651/1‑149 (MQ=255)

        agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATCCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATTccgtg                                                                                                                                          <  1:282524/148‑4 (MQ=255)

            gAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGAGACTTTAATCAC‑GGAATAATTCCCGTGATgg                                                                                                                                     >  2:176086/1‑148 (MQ=255)

             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATT‑CCGTGAttaaa                                                                                                                                   <  3:732716/149‑5 (MQ=255)

             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGa                                                                                                                                    >  2:898165/1‑148 (MQ=255)

             aaaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGa                                                                                                                                    >  1:661181/1‑148 (MQ=255)

              aaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGata                                                                                                                                  <  3:265712/149‑1 (MQ=255)

                 tCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaat                                                                                                                               >  3:180059/1‑149 (MQ=255)

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

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

                     agACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GAAATAATTCCCGTGATGGATAACTCAgg                                                                                                                           >  4:618702/2‑149 (MQ=255)

                      gACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAgg                                                                                                                           >  1:193746/1‑148 (MQ=255)

                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTa                                                                                                                        <  1:208553/148‑1 (MQ=255)

                         cTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTa                                                                                                                        <  1:176071/148‑1 (MQ=255)

                            gcggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGTGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGa                                                                                                                     >  1:690283/1‑148 (MQ=255)

                               gcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTcc                                                                                                                 <  4:704290/149‑1 (MQ=255)

                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgc                                                                                                               >  3:267882/1‑149 (MQ=255)

                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                              >  2:393436/1‑148 (MQ=255)

                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                              >  2:134482/1‑148 (MQ=255)

                                   tCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGt                                                                                                             <  3:859434/149‑1 (MQ=255)

                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTa                                                                                                            <  2:508782/148‑1 (MQ=255)

                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTa                                                                                                            <  2:280368/148‑1 (MQ=255)

                                            cTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcg                                                                                                     <  2:947207/148‑1 (MQ=255)

                                             tCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgcg                                                                                                   >  4:778428/1‑149 (MQ=255)

                                              cGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgcg                                                                                                   >  1:527144/1‑148 (MQ=255)

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

                                                 gACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCggggg                                                                                               <  4:414437/149‑1 (MQ=255)

                                                    cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGTACTTTAATCAC‑GGAATAATTCCCGGGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                            >  3:745082/1‑149 (MQ=255)

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

                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGTTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTa                                                                                                            >  2:793069/1‑132 (MQ=255)

                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                            >  1:914792/1‑148 (MQ=255)

                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                            >  2:600090/1‑148 (MQ=255)

                                                     aTAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                            >  1:666361/1‑148 (MQ=255)

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

                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTa                                                                                                                        >  4:63407/1‑119 (MQ=255)

                                                       aaaTTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTAGCGCGGGGGCGTaa                                                                                          <  2:844283/148‑1 (MQ=255)

                                                           tACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                      >  2:955790/1‑148 (MQ=255)

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

                                                            aCCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                                    <  4:453729/149‑1 (MQ=255)

                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTTGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGc                                                                                              >  4:513070/1‑138 (MQ=255)

                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGc                                                                                              <  3:513070/138‑1 (MQ=255)

                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTgg                                                                                    <  1:484750/148‑1 (MQ=255)

                                                             ccGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGt                                                                                   >  4:633778/1‑149 (MQ=255)

                                                                  tCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAATTCTGGTTATaa                                                                              <  3:542527/149‑1 (MQ=255)

                                                                   cGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                              >  2:428491/1‑148 (MQ=255)

                                                                       ttttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCg                                                                          <  2:56749/148‑1 (MQ=255)

                                                                        tttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGt                                                                         >  2:292609/1‑148 (MQ=255)

                                                                           cAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                     <  3:877307/149‑1 (MQ=255)

                                                                           cAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAtt                                                                     <  3:203879/149‑1 (MQ=255)

                                                                           cAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAt                                                                      >  1:605707/1‑148 (MQ=255)

                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTa                                                                                                                        >  2:356712/1‑95 (MQ=255)

                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTa                                                                                                                        <  1:356712/95‑1 (MQ=255)

                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                  <  3:624675/149‑1 (MQ=255)

                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                              <  3:738154/149‑1 (MQ=255)

                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                              <  3:336073/149‑1 (MQ=255)

                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                              >  3:385899/1‑149 (MQ=255)

                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                               >  1:217407/1‑148 (MQ=255)

                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                               <  2:23624/148‑1 (MQ=255)

                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                              >  2:849599/1‑148 (MQ=255)

                                                                                     aaaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggca                                                           >  4:240024/1‑149 (MQ=255)

                                                                                      aaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                     >  1:150266/1‑122 (MQ=255)

                                                                                      aaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTg                                                                                     <  2:150266/122‑1 (MQ=255)

                                                                                      aaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggcag                                                          >  3:218122/1‑149 (MQ=255)

                                                                                        ccGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGt                                                         >  1:613499/1‑148 (MQ=255)

                                                                                        ccGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGt                                                         <  1:898165/148‑1 (MQ=255)

                                                                                              gcgctgcgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCATGCAGTTTGGTCa                                                  >  3:407730/1‑149 (MQ=255)

                                                                                                    cgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCa                                            <  4:856636/149‑1 (MQ=255)

                                                                                                    cgcCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCa                                            <  3:778428/149‑1 (MQ=255)

                                                                                                        gtgtGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTa                                                                                                            <  1:808733/81‑1 (MQ=255)

                                                                                                        gtgtGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAgcg                                         >  1:529940/1‑148 (MQ=255)

                                                                                                        gtgtGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGACAAATAAGGGTAAGATTCCGCGTa                                                                                                            >  2:808733/1‑81 (MQ=255)

                                                                                                           tgGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCgca                                     <  3:375771/149‑1 (MQ=255)

                                                                                                                     tGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATGAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAg                           >  3:545146/1‑149 (MQ=255)

                                                                                                                      gAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAg                           <  2:510757/148‑1 (MQ=255)

                                                                                                                        aaaaGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGt                        >  4:803196/1‑149 (MQ=255)

                                                                                                                         aaaGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCgg                             >  1:545235/1‑143 (MQ=255)

                                                                                                                         aaaGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCgg                             <  2:545235/143‑1 (MQ=255)

                                                                                                                           agaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                     >  3:887963/1‑149 (MQ=255)

                                                                                                                            gaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                     >  2:411168/1‑148 (MQ=255)

                                                                                                                              agGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa                   <  1:41412/148‑1 (MQ=255)

                                                                                                                              agGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCATAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaaa                   >  1:859617/1‑148 (MQ=255)

                                                                                                                                   cTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTAc              >  1:769921/1‑148 (MQ=255)

                                                                                                                                    tttAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt            >  4:937316/1‑149 (MQ=255)

                                                                                                                                    tttAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt            >  4:447465/1‑149 (MQ=255)

                                                                                                                                     ttAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt            >  1:458993/1‑148 (MQ=255)

                                                                                                                                      tAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTg           <  2:217407/148‑1 (MQ=255)

                                                                                                                                        aTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTc                                                                                       >  1:72317/1‑70 (MQ=255)

                                                                                                                                        aTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTc                                                                                       <  2:72317/70‑1 (MQ=255)

                                                                                                                                         tcacgggAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGGCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa       >  4:803297/6‑149 (MQ=255)

                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGGAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa       >  4:678493/1‑149 (MQ=255)

                                                                                                                                           acgggAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg     >  3:658786/4‑149 (MQ=255)

                                                                                                                                           aC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGAt      >  1:455051/1‑148 (MQ=255)

                                                                                                                                            c‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                             >  1:749935/1‑76 (MQ=255)

                                                                                                                                            c‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAg                                                                             <  2:749935/76‑1 (MQ=255)

                                                                                                                                            c‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg     >  2:188164/1‑148 (MQ=255)

                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg    <  1:13753/148‑1 (MQ=255)

                                                                                                                                              ggAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTAGAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATgg    <  2:690283/148‑1 (MQ=255)

                                                                                                                                               gAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGc   >  2:665052/1‑148 (MQ=255)

GGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCG  >  NC_000913/192433‑192723

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