BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	129,798	A→T	9.4%	I488N (ATC→AAC)	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	129,798	0	A	T	9.4%	311.6 / 3.4	138	I488N (ATC→AAC)	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base A (58/59); new base T (2/12); total (67/71)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GCTCTGGTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGT  >  NC_000913/129649‑129945

gcgctGGGGTTGCGGACTGTCTACGCTCACGCGATGCGGCTCGTCGCTGAGGGTGCTGCGGGTTTGTCGGCATTTTCTCGCGAACCGCCTGGCGCTGCGCAGGCGAGGCGGACGGAATGCGCTCACGCGCTTCCTCTCTTCGCTGCTGGa                                                                                                                                                     <  1:942578/147‑1 (MQ=255)

 ctctGGTGATGCGGACTGAATACGCTCACGCGCTGCCTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCCGGCGAGGCGGACTGAAAGCGCTCACGCGCTGCCGCTCGTCGCTGCTGGAt                                                                                                                                                    >  2:1389192/1‑150 (MQ=18)

  tctGGTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACAGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATACGCTCACGAGCTGCCTCTATTAGGTGCTGGTTc                                                                                                                                                   >  2:1906126/1‑150 (MQ=18)

  tctGGTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTGCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATc                                                                                                                                                   <  1:1421237/150‑1 (MQ=255)

      gTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTGCGCTGACGCGCTGCCTCGATTCGATGCTGGTTCTTCg                                                                                                                                               >  2:62813/1‑150 (MQ=18)

       tGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGATGCTGGATCTTCgg                                                                                                                                              >  2:775998/1‑150 (MQ=255)

        gATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCggg                                                                                                                                             <  2:491683/150‑1 (MQ=255)

        gATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGGCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCggg                                                                                                                                             <  1:1707846/150‑1 (MQ=255)

          tGCGGAGTGAATACGCTCTCGCGCTGCGTCTCTTCGCTGCTGCTTCTGCGGGATAGTCTGCATTTTCTCGAGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGtt                                                                                                                                           <  1:167559/150‑1 (MQ=25)

          tGCGGAATGAAAACGCTCCCGCGCTGCTTCTCTTCGCTGCTGGTTCTGCGGGTGATTCTGCAGTTTCTCGCGAACCTCCTGGCGCTGCTCTGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGtt                                                                                                                                           <  2:325380/150‑1 (MQ=255)

            ccgACTGAACTCGCTCCCGGGCGGCCTCTCTTGGCTGGGGGGTCTGGGGGCGCGGCTGCACTATCTCGCTAAACTCCTGGCGCTGCTCCGGGGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAg                                                                                                                                         <  2:274447/148‑1 (MQ=255)

             ggACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGGTCTTCGGGTTAGt                                                                                                                                        >  1:863503/1‑150 (MQ=255)

                  gAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACAGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAGCGCTCACGCGCTGCCTCTCTTCGGTGCTGGATCTTCGGGTTAGTCTGcg                                                                                                                                   >  2:925664/1‑149 (MQ=255)

                     tACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTGTTCGCTGCTGGATCTTCGGGTGAGTCTGCCTtc                                                                                                                                >  1:1016475/1‑150 (MQ=255)

                         cTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCAGGCTtctctc                                                                                                                            >  1:1506566/1‑150 (MQ=255)

                            ccgcgcgGCGTCTCTTGGCTGCTGGTGCTGCGGGTTCGGCTGCCTTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTcgcg                                                                                                                         <  1:351295/149‑1 (MQ=37)

                            ccgcgcgGCGTCTCTTCGCTGCTGTTTCTGCGGGGTAGTCTGCATTTTCTCGCTAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGTATTCTCTcgcg                                                                                                                         <  2:1310278/149‑1 (MQ=37)

                            aCGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTcgcg                                                                                                                         <  1:408310/150‑1 (MQ=255)

                                     gTCTCTTCGCTGCTGGTTCTGCGGGCTAGGCTGCATTTTCTCCCTAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTg                                                                                                                <  2:387518/150‑1 (MQ=255)

                                      gctctTCGCTGCTGGTTCTGCGGGGTAGGCTGAACTTTCTCGCGAACCGCCTGGCGCTGCGCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGGTAGTCTGCATTCTCTCGCGAACTGCCTgg                                                                                                               <  2:843433/149‑1 (MQ=255)

                                       ctctTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTAGGGTTAGTCAGCATTCTCTCGCGAACTGCCTGgc                                                                                                              >  2:650456/1‑150 (MQ=255)

                                       ctctTAGCTGCTGGTTCTGAGGGGTAGGCTGCATTTTCTCGCTAACCGCCTGGCGCTGCTCAGTCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGgc                                                                                                              <  2:1106792/150‑1 (MQ=255)

                                         ctTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGCGGCGGACTGATTGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGGTCTTCGGGTTGGTCTGCATGCTCTCGGGACCTTCCTGgctc                                                                                                            >  1:368802/1‑148 (MQ=255)

                                            cGCTGCTGTTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCgctgc                                                                                                         <  2:429421/150‑1 (MQ=255)

                                             gccgctgGTGCTGCGGGCGAGTCTGCACTGTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCgctgct                                                                                                        <  2:131937/147‑1 (MQ=38)

                                              ctgctgTTTCTGCGGGCTAGGATGCAATTTCTCCTTAACCTCCTGGCGCTGCTCAGGCGATGCGGACTGAATCCGCTGACGGGCTGCCTCTCTTCGCTGCTGTTTCTTCGGGTTAGTCTGCATTGTCTCGCTAACTGCCTGGCGCTGCTc                                                                                                       <  2:1811697/150‑1 (MQ=255)

                                              ctgctgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTc                                                                                                       >  1:657710/1‑150 (MQ=255)

                                               tgctgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAATCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCa                                                                                                      >  2:1818851/1‑150 (MQ=255)

                                               tgcggtttCTGCGGGCTCGTCTGCATTTTCTCGCGAACCTCCTGGCGCTGCTCAGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCa                                                                                                      <  2:1190994/144‑1 (MQ=38)

                                                gctgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCag                                                                                                     >  1:1421259/1‑150 (MQ=255)

                                                 ctgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCGTCTCTTAGCTGCTGGCTCTTAGGTTTAGGCTGCATGCTCTCGCGAACTGCCTGGAGCTGCTCagg                                                                                                    >  2:115510/1‑150 (MQ=21)

                                                  tgGTTCTGCGGGTTAGTCTGAATTTTCTCGCTAACCTCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCaggc                                                                                                   <  2:1749457/150‑1 (MQ=255)

                                                   ggTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCaggcg                                                                                                  >  1:643595/1‑150 (MQ=255)

                                                         gCGGGGTCGGCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATCCGCTCCCGCGCTGCCTCTCTTCGCTGCTGGATCTGCGGGTGAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCgg                                                                                            <  1:183985/150‑1 (MQ=18)

                                                         gCGGGCTAGTCTGCATTTTCTCCCTAACCGCCTGGCGCTGCGCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCgg                                                                                            <  2:435453/150‑1 (MQ=255)

                                                            ggTTCGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACt                                                                                         <  1:1810191/150‑1 (MQ=255)

                                                             gTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTAGATGCTGGATCTTAGGTTGAGGAGGCATTCTCTCGCGAACTGCCTCGCGCTGCGCAGGCGCTGGGGACTg                                                                                        >  2:1189698/1‑150 (MQ=25)

                                                              ctcggcggCATTTTCTCCCGAACCTCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTGCGGGTGAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGa                                                                                       <  1:90702/143‑1 (MQ=17)

                                                                 gTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGCTCTTCGGGTTAGTCTGCATTCTCTCGAGAACTGCCTGGCGCTTCTCAGGGGAGGCGGACTGATaa                                                                                    >  1:243227/1‑150 (MQ=255)

                                                                  tCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCTTTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAAc                                                                                   >  1:597568/1‑150 (MQ=255)

                                                                  tCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCGCAGGCGAGGCGGACTGATAAc                                                                                   >  1:63643/1‑150 (MQ=255)

                                                                    tGCATTTTCTCGCTAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGc                                                                                 <  1:1486795/150‑1 (MQ=255)

                                                                     gCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCt                                                                                >  1:1403588/1‑150 (MQ=255)

                                                                        ttttCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCCCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGGACTGCCTGGCGCTGCTCAGGAGTGGCGGACTGATAACGCTGaa                                                                             >  1:307693/1‑149 (MQ=255)

                                                                        ttttCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGAc                                                                             >  2:1856379/1‑150 (MQ=255)

                                                                            ctcGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgc                                                                         >  1:745928/1‑150 (MQ=255)

                                                                              cgcgAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcgg                                                                       >  1:1553350/1‑150 (MQ=255)

                                                                                cgAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcggcg                                                                     <  1:1162941/150‑1 (MQ=255)

                                                                                  aaCCGCCTCGCGCTGCTCAGGCGACGCGGACTGAATGCGCTCACGCGCTGCCTCGCTTCGCTGCTGGATCTTCGCGTTAGGCTGCATTCTCTCGAGAACTGCCTGGCTCTGCTCAGGCGAGGCGCACTGAATACGCTGCCGAGCGGCGTc                                                                   >  2:254185/1‑150 (MQ=255)

                                                                                   cccGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTcc                                                                  <  1:898982/149‑1 (MQ=255)

                                                                                      gCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCttt                                                               >  1:1053541/1‑150 (MQ=255)

                                                                                      gCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCttt                                                               >  1:1403770/1‑150 (MQ=255)

                                                                                        cTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTACGCTGCCGAGCGGCGTCCTTttg                                                             >  1:675644/1‑150 (MQ=255)

                                                                                         tGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGCGCGGCGTCCTTttgt                                                            >  1:1746772/1‑150 (MQ=255)

                                                                                          ggCGCTGCTCAGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTttgtt                                                           <  2:1447414/150‑1 (MQ=255)

                                                                                             gctgctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCt                                                        >  2:1779908/1‑150 (MQ=255)

                                                                                              ctgctCCGGCGAGGCGGACTGAATGCGCTGACGCGCTGCCTCTCTTCGCTGCTGGTTCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACAAGCGGCGTCCTTTTGTTGCTg                                                       <  1:16124/150‑1 (MQ=255)

                                                                                              ctgctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGACTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTg                                                       <  2:1710666/150‑1 (MQ=255)

                                                                                               ggctCTGGTGAGGCGGACTGAATTCGCGGACGGGCTGCCTGTCTTGGGTGCTGGTTCTTGGGGTTAGTCTGCCTTCTCTCGAGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTgg                                                      <  2:461881/149‑1 (MQ=17)

                                                                                                gctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGCCGAGCGGCGTCCGTTTGTTGCTggg                                                     >  1:482951/1‑150 (MQ=255)

                                                                                                 ctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGt                                                    >  2:468180/1‑150 (MQ=255)

                                                                                                  cctttctaggcgGACTGAATTCGATCACGGGCTGCCTCACTTTGCTGTTGGTTCTGTGGGTTATTCTGAATTTTCTCGATAAACTCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTc                                                   <  2:110982/143‑1 (MQ=12)

                                                                                                   cAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGtga                                                  >  1:1701887/1‑148 (MQ=255)

                                                                                                    aggcgaggcgGACTGAATCCGCTCCCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTCGTCTGCATTCTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAg                                                 <  2:1372523/150‑1 (MQ=255)

                                                                                                     ggcgaggcgGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTGTTGTTGCTGGGGCAGt                                                >  1:201394/1‑150 (MQ=255)

                                                                                                      gcgaggcgTACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGTTCTGCGGGTTTGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTg                                               <  1:1395008/150‑1 (MQ=21)

                                                                                                      gcgaggcgGACTTAATGCGCTCACGCTCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTg                                               <  1:1069275/150‑1 (MQ=255)

                                                                                                       ggaggcgGACTGAATCCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTgg                                              <  2:1819101/149‑1 (MQ=255)

                                                                                                          ggcgGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTg                                           <  1:536484/150‑1 (MQ=255)

                                                                                                           gctgACTGAATGCGCTCACGCGCGGGCTCTCTTCGCTGCGGTATCTTCGGGTTAGGCTGTATTCTCTCGCTAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTgg                                          <  1:948534/147‑1 (MQ=25)

                                                                                                            cgGACTGAATGCGCTCACGCGCGGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGc                                         <  2:1453126/150‑1 (MQ=255)

                                                                                                             ggACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCg                                        >  2:719036/1‑150 (MQ=255)

                                                                                                             ggACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCg                                        >  2:1524962/1‑150 (MQ=255)

                                                                                                                  gtatccgcTCACGCGCTGCCTCTCTTCTCTGCTGTTTCTTCGGGTTAGTCTGGATTGTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGc                                   <  2:84162/145‑1 (MQ=21)

                                                                                                                  gAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCTTTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGc                                   <  2:597568/150‑1 (MQ=255)

                                                                                                                  gAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGc                                   <  2:1446492/150‑1 (MQ=255)

                                                                                                                  gAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGc                                   >  1:1755070/1‑150 (MQ=255)

                                                                                                                  gAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCGTTAGTTGCTGGGTCAGTGGTTGGCGACGGc                                   >  1:1669937/1‑150 (MQ=255)

                                                                                                                   aaTGCGCTCACGCGCGGCCTCTCTCCGCTGCTGGAGCTTCCGGCTCGTCTGCATTCTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCt                                  <  1:671682/150‑1 (MQ=25)

                                                                                                                    actctcccACGCGCTGCCTCTCTTGTCTGATGTTGCTGCGGGGTATGCTGCATTTTTTCGTTAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTg                                 <  2:303810/143‑1 (MQ=18)

                                                                                                                         cTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTc                            >  1:1546241/1‑150 (MQ=255)

                                                                                                                            aCGCGCTGCCTCTCTTCGCTGCGGGTTCTTCGGGTTAGTCTGCATTCTCTCGCTAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTg                         <  2:1634533/150‑1 (MQ=255)

                                                                                                                             cgggcTGCCTTTCTTAGGTGCTGGATCTGCGGGTTAGTCTGCATTATATCGCGAACTTCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTgg                        <  2:1256649/147‑1 (MQ=37)

                                                                                                                             cgcgCTGCCTCTCTTCGCTGCTGTTTCTTCGGGTTAGTCTGCATTTTCTCGCTAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTgg                        <  2:681104/150‑1 (MQ=255)

                                                                                                                             cgcgCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTgg                        >  2:499789/1‑150 (MQ=255)

                                                                                                                             cgcgCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTgg                        >  1:321914/1‑150 (MQ=255)

                                                                                                                              gcgcTGCCTCTCTTCGCTGCTTGATCTGCGGGTTAGGCTGCAGTCTCTCGCGAACTTCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGa                       <  1:1211555/150‑1 (MQ=37)

                                                                                                                              gcgcTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGa                       <  1:1856379/150‑1 (MQ=255)

                                                                                                                                   gCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAgtcg                  >  2:297299/1‑150 (MQ=255)

                                                                                                                                     ctctctTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACGGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAgtcgtc                <  1:1466276/150‑1 (MQ=255)

                                                                                                                                       ctctTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGTGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGTTTGGCGACGGCTGAAGTCGTGGAAGTCGTCAt              >  2:89638/1‑150 (MQ=255)

                                                                                                                                          ttCGCTGCGGGATCTGCGGGCTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGc           <  1:20018/150‑1 (MQ=255)

                                                                                                                                           tCGCTGCTGGATCTTCGGGTTCGTCGGCATTCTCTCGCGAACTGCCTGGCGCTGATGAGTCGAGGCTGACTGAAAACGCTTACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCt          <  1:903952/150‑1 (MQ=255)

                                                                                                                                           tCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGATGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCt          >  2:1549936/1‑150 (MQ=255)

                                                                                                                                             gctgctGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGTGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTcc        >  2:881613/1‑150 (MQ=255)

                                                                                                                                             gctgctGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGTGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTcc        >  2:932198/1‑150 (MQ=255)

                                                                                                                                             gctgctGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTcc        >  2:892431/1‑150 (MQ=255)

                                                                                                                                                gctgGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACAGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATATCTCCCAt     >  2:1541626/1‑150 (MQ=255)

                                                                                                                                                 ctgGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATa    >  1:445083/1‑150 (MQ=255)

                                                                                                                                                  gggATCTGCGGGTTAGTCTTCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGTCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAg   <  1:580720/149‑1 (MQ=255)

GCTCTGGTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGT  >  NC_000913/129649‑129945

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