BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,517,484	C→G	17.1%	I32M (ATC→ATG)	ydcX →	DUF2566 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,517,484	0	C	G	17.1%	409.6 / 75.8	160	I32M (ATC→ATG)	ydcX	DUF2566 family protein
Reads supporting (aligned to +/- strand): ref base C (67/64); new base G (7/20); total (75/85)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 165 total aligned reads displayed.

GTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACG  >  NC_000913/1517385‑1517583

gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                     <  2:1776616/100‑1 (MQ=255)

 tACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAAGGCATcc                                                                                                    >  2:1925289/1‑100 (MQ=255)

 tACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATcc                                                                                                    >  1:569461/1‑100 (MQ=255)

      gTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTc                                                                                               <  1:4254469/100‑1 (MQ=255)

      gTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGATTc                                                                                               <  1:2416774/100‑1 (MQ=255)

          ctctATCTACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTTGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTat                                                                                           <  2:1288864/100‑2 (MQ=255)

          ctctATCAACCTATGCTTGTTTTTTATGTGGTTATTGCATCAATCGCATTTCTTATCACCTAGTTACTTTCTCACGATAAGCAACGTATTCGTTTCTTat                                                                                           <  1:1335329/100‑2 (MQ=255)

          ctctACCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTAACGATAAGAAACGGATGCGGGTCTTat                                                                                           <  1:3587415/100‑2 (MQ=255)

           tctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAg                                                                                          <  2:51423/100‑1 (MQ=255)

           tctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAg                                                                                          <  2:666123/100‑1 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:753379/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:3019401/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  1:616186/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  1:934475/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:2492542/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:2350160/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:70500/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:2178491/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  1:1660323/100‑4 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgc                                                                                         >  2:3668462/1‑100 (MQ=255)

            ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgc                                                                                         <  2:4159013/100‑1 (MQ=255)

             tATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcg                                                                                        <  2:3858711/100‑1 (MQ=255)

             tATCAACACATGCTTGTTCTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcg                                                                                        >  1:4404805/1‑100 (MQ=255)

              aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcgtg                                                                                       <  2:2750345/100‑6 (MQ=39)

              aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                         <  2:2194621/98‑4 (MQ=255)

              aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcgc                                                                                       >  1:4228834/1‑100 (MQ=255)

                naaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCt                                                                                     >  2:2364952/2‑100 (MQ=255)

                cAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCt                                                                                     >  2:149963/1‑100 (MQ=255)

                 aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcgtgaga                                                                                    <  2:3004934/100‑9 (MQ=38)

                 aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCt                                                                                     >  1:4370142/1‑99 (MQ=255)

                   cacaTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTc                                                                                  >  2:3709017/1‑100 (MQ=255)

                   cacaTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTc                                                                                  <  2:3409807/100‑1 (MQ=255)

                     caTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtg                                                                                >  2:1876828/1‑100 (MQ=255)

                      aTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtgg                                                                               <  2:3761063/100‑1 (MQ=255)

                         cTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtgg                                                                            >  2:3434747/1‑100 (MQ=255)

                         cTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtgg                                                                            <  1:1871604/100‑1 (MQ=255)

                           tGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgggg                                                                          >  1:2109813/1‑100 (MQ=255)

                           tGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgggg                                                                          <  2:525722/100‑1 (MQ=255)

                           tGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgggg                                                                          >  1:2017091/1‑100 (MQ=255)

                            gTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggggg                                                                         >  2:127115/1‑100 (MQ=255)

                             ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGc                                                                        >  1:3953231/1‑100 (MQ=255)

                             ttttttAAGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGc                                                                        <  1:274274/100‑1 (MQ=255)

                              tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                       >  1:1306993/1‑100 (MQ=255)

                              tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                       <  2:1711970/100‑1 (MQ=255)

                               ttttATGCGGTTATGGCAGCAATCGCATTTCTTCTCACCTGGGTTCTTTCTCAAGATAAGAAATGCCTACGGTTCTTAATCGCCTTTCTGGTGGAGGCaa                                                                      >  1:1381064/1‑100 (MQ=255)

                                 ttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACa                                                                    >  1:1765058/1‑100 (MQ=255)

                                   aTGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATg                                                                  >  2:1370785/1‑100 (MQ=255)

                                    tGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATgg                                                                 >  1:2291875/1‑100 (MQ=255)

                                     gCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGc                                                                <  1:296031/100‑1 (MQ=255)

                                      cGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGCAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGCGCAACATGAcc                                                               >  1:2798047/1‑100 (MQ=255)

                                      cGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATg                                                                  <  1:3573384/97‑1 (MQ=255)

                                      cGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGcc                                                               <  2:2109820/100‑1 (MQ=255)

                                      cGGTTATGGCAGCAATCGCATTTCTTATAACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGcc                                                               <  1:76577/100‑1 (MQ=255)

                                          tATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATg                                                           <  2:3241702/100‑1 (MQ=255)

                                          tATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATg                                                           <  2:3585625/100‑1 (MQ=255)

                                            tGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAg                                                         <  2:3243260/100‑1 (MQ=255)

                                            tGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAg                                                         <  2:4080915/100‑1 (MQ=255)

                                              gcagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGtt                                                       <  2:245442/100‑1 (MQ=255)

                                              gcagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGtt                                                       <  2:1219864/100‑1 (MQ=255)

                                                agcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTc                                                     >  1:638322/1‑100 (MQ=255)

                                                agcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTc                                                     >  1:3286685/1‑100 (MQ=255)

                                                   aaTCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCg                                                  <  1:2596628/100‑1 (MQ=255)

                                                       gCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTgg                                              <  2:1862094/100‑1 (MQ=255)

                                                         aTTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcg                                            >  1:4116633/1‑100 (MQ=255)

                                                          tttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcgc                                           <  2:4291451/100‑1 (MQ=255)

                                                               tATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgttg                                      <  2:2567942/100‑1 (MQ=255)

                                                                 tCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTgttgtt                                    <  2:2394386/100‑1 (MQ=255)

                                                                    ccTTGTTTCTTTCTCCCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGtttt                                  <  1:931778/99‑1 (MQ=255)

                                                                    ccTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTc                                 >  2:3789098/1‑100 (MQ=255)

                                                                    ccTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTc                                 >  1:3372173/1‑100 (MQ=255)

                                                                     cTGGTTTCTTTCTCACGATAAGAATCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCa                                <  2:4246164/100‑1 (MQ=255)

                                                                      tGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCAc                               <  2:2872455/100‑1 (MQ=255)

                                                                      tGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCAc                               >  1:3919713/1‑100 (MQ=255)

                                                                         tttctttcTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGt                            >  1:3827789/1‑100 (MQ=255)

                                                                            ctttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtttt                         >  2:4299761/1‑100 (MQ=255)

                                                                            ctttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtttt                         <  2:1806774/100‑1 (MQ=255)

                                                                             tttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTa                        <  2:1060203/100‑1 (MQ=255)

                                                                             tttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTa                        >  2:105229/1‑100 (MQ=255)

                                                                               tctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                      <  2:4363179/100‑1 (MQ=255)

                                                                               aaggcgcttAAGAAACGGATGCGTATCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTATTTTCACTGTTTTaaa                      >  1:2579035/9‑100 (MQ=38)

                                                                                 tcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATc                    >  2:914550/1‑100 (MQ=255)

                                                                                  cACGATAAGAAACGCATCCGTTTCTTCAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCt                   >  2:2677483/1‑100 (MQ=255)

                                                                                  cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCt                   >  1:4211937/1‑100 (MQ=255)

                                                                                   cgcttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCtt                  >  2:3490991/5‑100 (MQ=39)

                                                                                   aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCtt                  <  2:3896487/100‑1 (MQ=255)

                                                                                   aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCtt                  >  1:247672/1‑100 (MQ=255)

                                                                                     gATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTAc                >  2:3215634/1‑100 (MQ=255)

                                                                                     cttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCt                   >  1:3184011/3‑97 (MQ=255)

                                                                                     cttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTAc                >  1:1747313/3‑100 (MQ=255)

                                                                                       tAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTc              >  2:1233312/1‑100 (MQ=255)

                                                                                       tAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGCTGATTTCACTGTTTTAAATCTTACTc              >  1:346824/1‑100 (MQ=255)

                                                                                        aaGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCa             >  2:3538727/1‑100 (MQ=255)

                                                                                         aGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCAc            >  1:1493437/1‑100 (MQ=255)

                                                                                          gAAACGCATCGGTTTCTTATGTGCCTTTCTGGTGGGTGTAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACGCAc            <  1:3684594/99‑1 (MQ=255)

                                                                                          gAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACt           >  2:4013556/1‑100 (MQ=255)

                                                                                             aCGCTTCCGTTTCTTAAGCGCCTTTCTGGTGGGGGTAAGATGGCCGATGAGTTTCCCGCTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGc        <  1:2405785/100‑1 (MQ=255)

                                                                                             aCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGc        <  2:2558742/100‑1 (MQ=255)

                                                                                                  tCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGAc   >  2:2725030/1‑100 (MQ=255)

                                                                                                  tCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGAc   <  2:1168569/100‑1 (MQ=255)

                                                                                                   ccGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg  >  1:2798805/1‑100 (MQ=255)

GTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACG  >  NC_000913/1517385‑1517583

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: `‑`