BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,908,486	G→C	5.7%	A394G (GCC→GGC)	pyrG ←	CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,908,486	0	G	C	5.7%	349.5 / 9.0	101	A394G (GCC→GGC)	pyrG	CTP synthetase
Reads supporting (aligned to +/- strand): ref base G (54/41); new base C (3/3); total (57/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 101 total aligned reads displayed.

GCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGC  >  NC_000913/2908341‑2908631

gcCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGCGTACGTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCaa                                                                                                                                                <  2:71295/149‑1 (MQ=255)

   gagaCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACAt                                                                                                                                             >  1:522907/1‑149 (MQ=255)

      aCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAGTCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGCCAACATgga                                                                                                                                          <  4:618488/149‑2 (MQ=255)

        gCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGCCAACATGgaga                                                                                                                                        <  1:803883/149‑1 (MQ=255)

          aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGc                                                                                                                                      >  4:204895/1‑149 (MQ=255)

          aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCGCCATGTTGCCAACATGgagaac                                                                                                                                      <  2:523828/149‑3 (MQ=255)

          aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCAGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGc                                                                                                                                      >  3:542812/1‑149 (MQ=255)

              tACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTTCAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCg                                                                                                                                         >  4:245188/1‑142 (MQ=255)

              tACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCg                                                                                                                                         <  3:245188/142‑1 (MQ=255)

               aCCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAAt                                                                                                                                 <  1:524526/149‑1 (MQ=255)

                   ccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAt                                                                                                                             <  4:263179/149‑1 (MQ=255)

                      agaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTaa                                                                                                                          <  3:401295/149‑1 (MQ=255)

                         tCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAACCAATTAACGc                                                                                                                       <  3:190184/149‑1 (MQ=255)

                           gCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCa                                                                                                                     >  1:420877/1‑149 (MQ=255)

                           gCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCa                                                                                                                     >  2:525457/1‑149 (MQ=255)

                               ttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGTCACTCGGTAATCAGCGCCACAACCGGGAACTTACAGTCTGGCACAAATTCCGTAGAGTTCGCGTTCTCCGTGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTg                                                                                                                 >  3:548515/1‑149 (MQ=255)

                               ttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTg                                                                                                                 >  4:135905/1‑149 (MQ=255)

                                  tcGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGGGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCAt                                                                                                              <  2:766691/149‑1 (MQ=255)

                                  tcGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCAt                                                                                                              >  3:756558/1‑149 (MQ=255)

                                      tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                          >  2:806633/1‑149 (MQ=255)

                                        cGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAg                                                                                                        <  4:853707/149‑1 (MQ=255)

                                             ttCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                   <  1:806633/149‑1 (MQ=255)

                                             ttCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                   >  2:109161/1‑149 (MQ=255)

                                                    ttGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAg                                                                                            >  4:698647/1‑149 (MQ=255)

                                                      gCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                          >  3:616201/1‑149 (MQ=255)

                                                       ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                          >  1:66914/1‑148 (MQ=255)

                                                       ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                          <  2:66914/148‑1 (MQ=255)

                                                        cGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                        >  1:856602/1‑149 (MQ=255)

                                                         gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                       <  4:458588/149‑1 (MQ=255)

                                                         gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                       >  3:365802/1‑149 (MQ=255)

                                                         gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                       <  3:832932/149‑1 (MQ=255)

                                                         gTTTTCATCGCGCCACTCGGTAATCAACGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                       >  4:186240/1‑149 (MQ=255)

                                                              cATCGCGCGACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCTCAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATAACTTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAata                                                                                  <  4:542812/149‑1 (MQ=255)

                                                               aTCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAatat                                                                                 >  4:12581/1‑149 (MQ=255)

                                                                   cgcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGtt                                                                             <  4:756329/149‑1 (MQ=255)

                                                                   cgcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGtt                                                                             <  2:671969/149‑1 (MQ=255)

                                                                      cACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTtctc                                                                          >  3:400614/1‑149 (MQ=255)

                                                                        cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAc                                                                        <  4:629314/149‑1 (MQ=255)

                                                                         tCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAcg                                                                       >  2:529810/1‑149 (MQ=255)

                                                                         tCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAcg                                                                       >  4:52561/1‑149 (MQ=255)

                                                                           ggTAATCAGCGCCACAACCGGGTACTTACAGTCTGGGACAAATTTCGTAGAGTTGGCGTTCTCGATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCGAATGCCCAGATAAGGAATATTGTCCTCAcgcg                                                                     >  3:550697/1‑149 (MQ=255)

                                                                           ggTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAcgcg                                                                     <  4:365802/149‑1 (MQ=255)

                                                                           ggTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAcgcg                                                                     >  1:80988/1‑149 (MQ=255)

                                                                            gTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACgcgc                                                                    >  3:530591/1‑149 (MQ=255)

                                                                            gTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACgcgc                                                                    >  2:422638/1‑149 (MQ=255)

                                                                            gTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACgcgc                                                                    >  2:143567/1‑149 (MQ=255)

                                                                                  aGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                              >  2:56728/1‑149 (MQ=255)

                                                                                     gcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgc                                                           >  2:512203/1‑149 (MQ=255)

                                                                                     gcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgc                                                           <  2:475543/149‑1 (MQ=255)

                                                                                     gcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgc                                                           >  4:142408/1‑149 (MQ=255)

                                                                                     gcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgc                                                           <  1:519480/149‑1 (MQ=255)

                                                                                       cacaACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGcc                                                                                                                      <  2:538470/88‑1 (MQ=255)

                                                                                       cacaACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGcc                                                                                                                      >  1:538470/1‑88 (MQ=255)

                                                                                         caACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTc                                                       <  3:445190/149‑1 (MQ=255)

                                                                                              gggTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAAt                                                  >  1:359848/1‑149 (MQ=255)

                                                                                                gTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCa                                                <  2:723559/149‑1 (MQ=255)

                                                                                                gTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCa                                                <  4:411562/149‑1 (MQ=255)

                                                                                                  aCTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATg                                              <  3:142408/149‑1 (MQ=255)

                                                                                                  aCTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATg                                              <  1:512203/149‑1 (MQ=255)

                                                                                                    ttACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGcc                                            >  2:269447/1‑149 (MQ=255)

                                                                                                    ttACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGcc                                            <  2:420877/149‑1 (MQ=255)

                                                                                                     tACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCt                                           <  3:176764/149‑1 (MQ=255)

                                                                                                       cAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTc                                         <  1:280254/149‑1 (MQ=255)

                                                                                                       cAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTc                                         >  1:519149/1‑149 (MQ=255)

                                                                                                         gTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTa                                       <  3:801075/149‑1 (MQ=255)

                                                                                                         gTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTa                                       >  2:23002/1‑149 (MQ=255)

                                                                                                            ttgCACAAATTCCGTAGAGTACGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCGGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAAGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGc                                    <  4:400614/147‑1 (MQ=255)

                                                                                                            tGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGc                                    >  1:513656/1‑149 (MQ=255)

                                                                                                              gCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCa                                  >  2:197546/1‑149 (MQ=255)

                                                                                                               cacaAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCAc                                 >  3:518326/1‑149 (MQ=255)

                                                                                                                  aaaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGAt                              >  4:388520/1‑149 (MQ=255)

                                                                                                                  aaaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGAt                              >  2:457828/1‑149 (MQ=255)

                                                                                                                   aaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATa                             >  1:779418/1‑149 (MQ=255)

                                                                                                                   aaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATa                             >  1:85012/1‑149 (MQ=255)

                                                                                                                         gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAATGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaa                       >  4:525183/1‑149 (MQ=255)

                                                                                                                               ttGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGcc                 >  2:199545/1‑149 (MQ=255)

                                                                                                                                tGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCa                <  1:197546/149‑1 (MQ=255)

                                                                                                                                tGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCa                <  2:522907/149‑1 (MQ=255)

                                                                                                                                tGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCa                <  3:795042/149‑1 (MQ=255)

                                                                                                                                    gTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa            >  4:594858/1‑149 (MQ=255)

                                                                                                                                    gTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa            >  2:775894/1‑149 (MQ=255)

                                                                                                                                    gTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa            >  1:338104/1‑149 (MQ=255)

                                                                                                                                    gTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa            >  1:552679/1‑149 (MQ=255)

                                                                                                                                     ttCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTAc           >  2:630430/1‑149 (MQ=255)

                                                                                                                                      tctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg          >  3:649038/1‑149 (MQ=255)

                                                                                                                                      tcgccgtGTTGCCAACATGGCGTGCGTAATAAAGTTACGCCACATGCCTACCCAGCCAAATACCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGATCGTACTCATGCCTTCTACGCCCCGATAGCCGAAACCGCCAGGTACg          >  4:309314/7‑149 (MQ=255)

                                                                                                                                      tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg          >  2:530852/4‑149 (MQ=255)

                                                                                                                                       ctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGa         <  1:525457/149‑1 (MQ=255)

                                                                                                                                       cgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGa         >  3:736155/3‑149 (MQ=255)

                                                                                                                                        tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg        <  1:584544/149‑1 (MQ=255)

                                                                                                                                        tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg        >  3:728059/1‑149 (MQ=255)

                                                                                                                                          cATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAAcc                    <  1:487656/135‑1 (MQ=255)

                                                                                                                                          cATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAAcc                    >  2:487656/1‑135 (MQ=255)

                                                                                                                                          cATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg        <  3:627359/147‑1 (MQ=255)

                                                                                                                                          cATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg        >  4:627359/1‑147 (MQ=255)

                                                                                                                                           aTGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGAt     <  1:529810/149‑1 (MQ=255)

                                                                                                                                            tGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGAtt    <  4:649038/149‑1 (MQ=255)

                                                                                                                                              ttGGTAACATGGCGAGCGTAATCAATTAACGCCACCTGCCTACCCAGCCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGAGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc  <  4:550697/149‑1 (MQ=255)

                                                                                                                                              ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc  >  2:486197/1‑149 (MQ=255)

                                                                                                                                              ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGCAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc  <  2:80988/149‑1 (MQ=255)

GCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGC  >  NC_000913/2908341‑2908631

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