BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,717,169	T→G	5.5%	L139R (CTC→CGC) ‡	ung →	uracil‑DNA‑glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,717,169	0	T	G	5.5%	439.2 / 15.0	128	L139R (CTC→CGC) ‡	ung	uracil‑DNA‑glycosylase
Reads supporting (aligned to +/- strand): ref base T (58/63); new base G (4/3); total (62/66)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

Only 100 of 129 total aligned reads displayed.

TTGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCAT  >  NC_000913/2717021‑2717314

ttGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCt                                                                                                                                                   >  1:233991/1‑149 (MQ=255)

    atatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCGGCCTCGGc                                                                                                                                               <  1:271270/149‑1 (MQ=255)

    atatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCGGCCTCGGc                                                                                                                                               >  2:271270/1‑149 (MQ=255)

       tGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTgg                                                                                                                                            <  1:671404/149‑1 (MQ=255)

         tataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGa                                                                                                                                          <  2:563228/149‑1 (MQ=255)

         tataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGa                                                                                                                                          >  4:852722/1‑149 (MQ=255)

          ataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGaa                                                                                                                                         >  3:11387/1‑149 (MQ=255)

           taAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                        <  4:511685/149‑5 (MQ=255)

            aaaGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtgg                                                                                                                                       <  1:664546/149‑6 (MQ=255)

              agagCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCt                                                                                                                                     >  1:469162/1‑149 (MQ=255)

               gagCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCtt                                                                                                                                    <  2:465985/149‑1 (MQ=255)

                   tGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGCTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCAcc                                                                                                                                >  3:16606/1‑149 (MQ=255)

                     gAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTAACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGa                                                                                                                              >  3:458775/1‑149 (MQ=255)

                         aTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCGGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                          <  3:612176/149‑1 (MQ=255)

                         aTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                          <  3:223041/149‑1 (MQ=255)

                          tACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAg                                                                                                                         >  3:895551/1‑149 (MQ=255)

                             tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCa                                                                                                                                  <  1:280295/137‑1 (MQ=255)

                             tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCa                                                                                                                                  >  2:280295/1‑137 (MQ=255)

                             tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGa                                                                                                                      >  2:661117/1‑149 (MQ=255)

                              aTTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGAt                                                                                                                     <  4:612279/149‑1 (MQ=255)

                                   gggCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGcc                                                                                                                <  2:297166/149‑1 (MQ=255)

                                   gggCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGcc                                                                                                                <  1:34660/149‑1 (MQ=255)

                                    ggCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCt                                                                                                               <  2:645317/149‑1 (MQ=255)

                                      cTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCGCGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGa                                                                                                             <  3:631830/149‑1 (MQ=255)

                                        tCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAtt                                                                                                           <  1:843686/149‑1 (MQ=255)

                                          aCCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGcc                                                                                                                >  1:373029/1‑142 (MQ=255)

                                          aCCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGcc                                                                                                                <  2:373029/142‑1 (MQ=255)

                                           cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAc                                                                                                        <  1:371471/149‑1 (MQ=255)

                                           cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCAACGATAAAGTGATCAGCCTGATTAAc                                                                                                        <  1:97850/149‑1 (MQ=255)

                                                 cTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCAt                                                                                                  >  4:405856/1‑149 (MQ=255)

                                                  tAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCt                                                                                                                                              >  1:393424/1‑104 (MQ=255)

                                                  tAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCt                                                                                                                                              <  2:393424/104‑1 (MQ=255)

                                                   aaTCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcg                                                                                                <  3:852722/149‑1 (MQ=255)

                                                     tCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgcg                                                                                              >  4:514727/1‑149 (MQ=255)

                                                      cATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGc                                                                                                                 <  2:889303/129‑1 (MQ=255)

                                                      cATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGa                                                                                             >  2:800386/1‑149 (MQ=255)

                                                       aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                          >  3:811695/1‑119 (MQ=255)

                                                       aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                          <  4:811695/119‑1 (MQ=255)

                                                       aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGaa                                                                                            >  2:847087/1‑149 (MQ=255)

                                                       aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGaa                                                                                            <  1:715820/149‑1 (MQ=255)

                                                             aTCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtg                                                                                      <  3:372708/149‑1 (MQ=255)

                                                                ttGAAAGCTGGGCGCGTCATGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAGCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                   >  4:581878/1‑149 (MQ=255)

                                                                ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                   >  3:168438/1‑149 (MQ=255)

                                                                  gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                 >  4:15419/1‑149 (MQ=255)

                                                                  gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                 >  3:118097/1‑149 (MQ=255)

                                                                  gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTAACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                 <  4:458775/149‑1 (MQ=255)

                                                                   aaaGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGttt                                                                                <  2:390541/149‑1 (MQ=255)

                                                                   aaaGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGttt                                                                                <  2:599096/149‑1 (MQ=255)

                                                                    aaGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                               <  4:222130/149‑1 (MQ=255)

                                                                         gggCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttgttg                                                                          >  1:392987/1‑149 (MQ=255)

                                                                          ggCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTtgttgt                                                                         >  2:726208/1‑149 (MQ=255)

                                                                            cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                       <  1:847087/149‑1 (MQ=255)

                                                                               gTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGCTAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGa                                                                    >  3:43167/1‑149 (MQ=255)

                                                                                  aGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCg                                                                 >  3:491916/1‑149 (MQ=255)

                                                                                   gggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGc                                                                <  1:479923/149‑1 (MQ=255)

                                                                                       gTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATgc                                                            <  1:519061/149‑1 (MQ=255)

                                                                                        ttCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATgcg                                                           >  1:533668/1‑149 (MQ=255)

                                                                                             cTACTCAATCCTGTGTTGGCTGCGCGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGCCTGGGAAACAGTCCCCTATACAGTGTCCAGGCTAATGATCCAGCATGGCGTAGGCGTGGTGTTTTTGTTGTTGGTATCCCATGCGCaaaa                                                      <  2:119381/149‑1 (MQ=255)

                                                                                             cTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                      >  2:463221/1‑149 (MQ=255)

                                                                                             cTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                      <  1:833307/149‑1 (MQ=255)

                                                                                              tACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaag                                                     >  4:190092/1‑149 (MQ=255)

                                                                                              tACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaag                                                     >  2:589574/1‑149 (MQ=255)

                                                                                                cTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaa                                                   >  2:62978/1‑149 (MQ=255)

                                                                                                 tCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaa                                                  >  3:824788/1‑149 (MQ=255)

                                                                                                 tCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaa                                                  <  3:743162/149‑1 (MQ=255)

                                                                                                   aaTACTGGGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGAGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAgg                                                >  2:799149/1‑149 (MQ=255)

                                                                                                     tACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAgggg                                              <  1:661117/149‑1 (MQ=255)

                                                                                                       cTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCg                                            >  3:20165/1‑149 (MQ=255)

                                                                                                        tgtgTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGa                                           <  3:850183/149‑1 (MQ=255)

                                                                                                        tgtgTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGa                                           >  2:7485/1‑149 (MQ=255)

                                                                                                          tgtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGAtt                                         >  4:627122/1‑149 (MQ=255)

                                                                                                                 ggTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATa                                  <  2:321461/149‑1 (MQ=255)

                                                                                                                  gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATaa                                 <  1:633994/149‑1 (MQ=255)

                                                                                                                  gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATaa                                 <  4:118097/149‑1 (MQ=255)

                                                                                                                   tACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAg                                <  3:190092/149‑1 (MQ=255)

                                                                                                                      gcgcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaa                                                  <  3:512863/128‑1 (MQ=255)

                                                                                                                      gcgcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaa                                                  >  4:512863/1‑128 (MQ=255)

                                                                                                                        gcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATa                                                                                                                            <  4:873806/52‑1 (MQ=255)

                                                                                                                        gcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATa                                                                                                                            >  3:873806/1‑52 (MQ=255)

                                                                                                                        gcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAggg                                               <  2:53022/129‑1 (MQ=255)

                                                                                                                        gcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAggg                                               >  1:53022/1‑129 (MQ=255)

                                                                                                                        gcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACg                           >  1:719598/1‑149 (MQ=255)

                                                                                                                            gTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCcat                       <  4:20165/149‑1 (MQ=255)

                                                                                                                             tCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCcatc                      <  4:16606/149‑1 (MQ=255)

                                                                                                                             tCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCcatc                      >  4:474811/1‑149 (MQ=255)

                                                                                                                                 gcgcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGt                  >  2:149141/1‑149 (MQ=255)

                                                                                                                                  cgcATTCCCTCGACAGCCTCGGCTGGGAAGCCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGGGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTa                 <  3:627133/149‑1 (MQ=255)

                                                                                                                                  cgcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTa                 <  3:581878/149‑1 (MQ=255)

                                                                                                                                   gcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTAc                <  4:188935/149‑1 (MQ=255)

                                                                                                                                        cccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGaaa           <  1:7485/149‑1 (MQ=255)

                                                                                                                                         ccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAg          >  3:497530/1‑149 (MQ=255)

                                                                                                                                         ccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAg          >  1:17651/1‑149 (MQ=255)

                                                                                                                                          cACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGc         >  2:614354/1‑149 (MQ=255)

                                                                                                                                           ttcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCa        >  3:138889/4‑149 (MQ=255)

                                                                                                                                            tcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAc       >  1:835174/3‑149 (MQ=255)

                                                                                                                                            tcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAc       >  2:476281/3‑149 (MQ=255)

                                                                                                                                             gCCAGCCTCGGCTGGGAAACCTTCACCGATACAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAcc      <  4:491916/149‑1 (MQ=255)

                                                                                                                                             gCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAcc      <  2:826765/149‑1 (MQ=255)

                                                                                                                                             gCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAcc      >  2:439103/1‑149 (MQ=255)

                                                                                                                                                 gCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCAt  >  3:218240/1‑149 (MQ=255)

TTGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCAT  >  NC_000913/2717021‑2717314

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