BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,410,416	A→G	17.8%	Y95C (TAC→TGC)	yjfI →	DUF2170 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,416	0	A	G	17.8%	360.5 / 69.8	168	Y95C (TAC→TGC)	yjfI	DUF2170 family protein
Reads supporting (aligned to +/- strand): ref base A (66/72); new base G (16/14); total (82/86)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 172 total aligned reads displayed.

TCAGCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTAC  >  NC_000913/4410317‑4410515

tCAGCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTa                                                                                                     <  2:3272085/100‑1 (MQ=255)

  aGCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTa                                                                                                     >  2:1375961/1‑98 (MQ=255)

  aGCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACa                                                                                                   >  1:3180569/1‑100 (MQ=255)

   gCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACAt                                                                                                  >  2:2345902/1‑100 (MQ=255)

     aaTCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTg                                                                                                <  2:3115116/100‑1 (MQ=255)

        ccAGATGGATTTAATGCCTTCTTATTAAGAAATCAGAAGATGATGCCGTGGTTAGCGGTATGGATTTCTAGCGTACAACAGGAAGAGTATTACAGTGttt                                                                                             <  2:160703/100‑1 (MQ=255)

        ccAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGttt                                                                                             <  1:227097/100‑1 (MQ=255)

        ccAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGttt                                                                                             <  2:682908/100‑1 (MQ=255)

         cAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTTTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGttt                                                                                             >  2:790023/1‑99 (MQ=255)

         cAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGtttt                                                                                            >  2:1189726/1‑100 (MQ=255)

             tGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCggg                                                                                        <  2:2736956/100‑1 (MQ=255)

             tGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCggg                                                                                        <  1:3557211/100‑1 (MQ=255)

              gAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCgggg                                                                                       >  2:1919361/1‑100 (MQ=255)

              gAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCgggg                                                                                       <  1:858533/100‑1 (MQ=255)

                aTTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCg                                                                                     >  1:616205/1‑100 (MQ=255)

                   tAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTa                                                                                  <  2:2119239/100‑1 (MQ=255)

                      tACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCg                                                                               >  1:3034719/1‑100 (MQ=255)

                      tACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCg                                                                               <  1:2732429/100‑1 (MQ=255)

                       aCCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGc                                                                              <  1:338935/100‑1 (MQ=255)

                        ccTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCt                                                                             >  2:3613856/1‑100 (MQ=255)

                         cttcttATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCtt                                                                            >  1:3501946/1‑100 (MQ=255)

                         cttcttATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCtt                                                                            >  2:1511507/1‑100 (MQ=255)

                            cttATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTaaa                                                                         <  2:1755961/100‑1 (MQ=255)

                               attaAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAAtct                                                                      >  1:2711313/1‑100 (MQ=255)

                                ttaAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAACGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAAtct                                                                      >  2:1510935/1‑99 (MQ=255)

                                  aaGAAATCAGCAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACCGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCCTAACtcttc                                                                    >  2:3353329/1‑99 (MQ=255)

                                  aaGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGGATTGTTTTCGGGGCGTTATCGCTTAAAtcttct                                                                   >  2:3798762/1‑100 (MQ=255)

                                   aGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTtctc                                                                  >  1:323747/1‑100 (MQ=255)

                                   aGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAAtcttct                                                                   >  2:2652154/1‑99 (MQ=255)

                                    gAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTctct                                                                 >  2:2642977/1‑100 (MQ=255)

                                     aaaTCAGAAGAGGATGCCGTTGTCATCGGTAGGGATCGCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAAtcttcaaga                                                                <  2:2135843/100‑5 (MQ=255)

                                      aaTCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTTCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTg                                                               <  2:2838328/100‑1 (MQ=255)

                                        tCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGaa                                                             >  1:458207/1‑100 (MQ=255)

                                         cAGAAGATGATGACGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAg                                                            <  2:127824/100‑1 (MQ=255)

                                          agaagaTGATGCCGTTGTCATCGGTAGGTATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGa                                                           <  1:1318850/100‑1 (MQ=255)

                                          agaagaTGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGa                                                           <  2:1018200/100‑1 (MQ=255)

                                          agaagaTGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGa                                                           <  2:1383517/100‑1 (MQ=255)

                                              gatgatGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATc                                                       <  2:3476861/100‑1 (MQ=255)

                                                 gatgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTg                                                    >  2:251359/1‑100 (MQ=255)

                                                 gatgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTg                                                    >  2:3280213/1‑100 (MQ=255)

                                                   tgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCt                                                  >  2:2487864/1‑100 (MQ=255)

                                                   tgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCt                                                  >  2:628518/1‑100 (MQ=255)

                                                   tgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCt                                                  >  2:79942/1‑100 (MQ=255)

                                                    gCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCtt                                                 >  2:3166717/1‑100 (MQ=255)

                                                     ccGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                                                >  2:1952237/1‑100 (MQ=255)

                                                     ccGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                                                <  1:2369721/100‑1 (MQ=255)

                                                     ccGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                                                <  1:3475683/100‑1 (MQ=255)

                                                      cGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                                                <  2:1937950/99‑1 (MQ=255)

                                                       gTTGTCATCGGTAGGGATCTCCAGCGTACAACTGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTgag                                              <  1:364067/100‑1 (MQ=255)

                                                       gTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTgag                                              <  1:903271/100‑1 (MQ=255)

                                                        ttGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTgaga                                             <  2:2278323/100‑1 (MQ=255)

                                                         tGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTgaga                                             <  2:2967582/99‑1 (MQ=255)

                                                         tGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGAt                                            <  2:3412443/100‑1 (MQ=255)

                                                            cATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAAc                                         >  2:1001378/1‑100 (MQ=255)

                                                            cATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAAc                                         <  2:3760480/100‑1 (MQ=255)

                                                             aTCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAAcc                                        <  1:1337841/100‑1 (MQ=255)

                                                                ggTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTCCTCTTGAAGATATCCTGCTTGAGATAACCTc                                      >  2:3035543/1‑99 (MQ=255)

                                                                 gTAGGGATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGc                                    <  1:2164370/100‑1 (MQ=255)

                                                                 gTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGc                                    >  1:3416833/1‑100 (MQ=255)

                                                                 gTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGc                                    <  2:403309/100‑1 (MQ=255)

                                                                 gTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGc                                    <  1:891807/100‑1 (MQ=255)

                                                                  tAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCt                                   >  1:773941/1‑100 (MQ=255)

                                                                    gggATCTCCAGCGTACAACAGGAATAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTgg                                 <  2:53610/100‑1 (MQ=255)

                                                                    gggATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTgg                                 <  1:3632392/100‑1 (MQ=255)

                                                                    gggATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTgg                                 <  1:967230/100‑1 (MQ=255)

                                                                    gggATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTgg                                 <  1:548333/100‑1 (MQ=255)

                                                                     ggATCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGt                                >  1:1641806/1‑100 (MQ=255)

                                                                     ggATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGt                                <  2:525290/100‑1 (MQ=255)

                                                                     ggATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGt                                <  2:2702593/100‑1 (MQ=255)

                                                                       aTCTCCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAg                              <  1:1490784/100‑1 (MQ=255)

                                                                       aTCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAg                              <  1:1164195/100‑1 (MQ=255)

                                                                         ctcCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGAt                            >  1:2658366/1‑100 (MQ=255)

                                                                         ctcCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGAt                            >  1:1587807/1‑100 (MQ=255)

                                                                          tcCAGCGTACAACAGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATa                           >  1:1980303/1‑100 (MQ=255)

                                                                          tcCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATa                           >  1:2458125/1‑100 (MQ=255)

                                                                           ccAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATaa                          >  1:598402/1‑100 (MQ=255)

                                                                            cAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAAc                         <  1:642478/100‑1 (MQ=255)

                                                                             agggTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCGTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAAc                         <  2:263803/96‑1 (MQ=255)

                                                                             aGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACg                        >  2:3551373/1‑100 (MQ=255)

                                                                             aGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACg                        >  1:1352654/1‑100 (MQ=255)

                                                                                 tACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCAtt                    >  1:1688585/1‑100 (MQ=255)

                                                                                 tACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCAtt                    >  1:1688476/1‑100 (MQ=255)

                                                                                  acaacaGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTg                   <  2:901182/100‑1 (MQ=255)

                                                                                   caacaGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTgg                  >  1:311010/1‑100 (MQ=255)

                                                                                    aacaGGAAGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGa                 <  1:1130558/100‑1 (MQ=255)

                                                                                    aacaGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGa                 >  2:1589179/1‑100 (MQ=255)

                                                                                     acaGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGAt                <  1:543025/100‑1 (MQ=255)

                                                                                        ggAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTg             >  1:2921603/1‑100 (MQ=255)

                                                                                         gAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTGCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTgg            <  2:2855493/100‑1 (MQ=255)

                                                                                          aaGAGTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGc           <  2:2383332/100‑1 (MQ=255)

                                                                                          aaGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGc           <  2:3049092/100‑1 (MQ=255)

                                                                                          aaGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGc           <  1:1236108/100‑1 (MQ=255)

                                                                                          aaGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGc           <  2:261012/100‑1 (MQ=255)

                                                                                            gagTATTGCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTg         <  1:1600159/100‑1 (MQ=255)

                                                                                            gagTATTACATTGTTTTCGGGGTGTTATCGCTTAAATCGTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTg         >  2:2346237/1‑100 (MQ=255)

                                                                                            gagTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGGAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTg         <  2:830738/100‑1 (MQ=255)

                                                                                              gTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGaa       >  1:3021114/1‑100 (MQ=255)

                                                                                              gTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGaa       <  1:696775/100‑1 (MQ=255)

                                                                                                   aCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTAc  >  1:2820627/1‑100 (MQ=255)

                                                                                                   aCATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAACATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTAc  >  2:2086560/1‑100 (MQ=255)

TCAGCAATCCAGATGAATTTAATACCTTCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTAC  >  NC_000913/4410317‑4410515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`