BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,517,481	C→G	17.2%	R31R (CGC→CGG)	ydcX →	DUF2566 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,517,481	0	C	G	17.2%	417.0 / 77.7	164	R31R (CGC→CGG)	ydcX	DUF2566 family protein
Reads supporting (aligned to +/- strand): ref base C (70/65); new base G (8/20); total (78/86)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 164 total aligned reads displayed.

CGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCT  >  NC_000913/1517384‑1517578

 gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGAtg                                                                                                >  2:303813/1‑99 (MQ=255)

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

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

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

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

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

           ctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTaa                                                                                      >  1:3863770/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

             ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                    <  2:70500/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)

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

               aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcgc                                                                                  >  1:2493226/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)

                  aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCtt                                                                               >  1:525721/1‑100 (MQ=255)

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

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

                   acacATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCttt                                                                              <  1:3723135/100‑1 (MQ=255)

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

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

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

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

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

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

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

                           ttGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTATCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTtgg                                                                      <  2:2243785/100‑1 (MQ=255)

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

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

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

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

                             gTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggggg                                                                    <  2:2523756/100‑1 (MQ=255)

                             gTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggggg                                                                    <  2:2708822/100‑1 (MQ=255)

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

                             gTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCAACTGGTTTCTTTCTCACGATCAGAAACGCATCCGTTTCTTAAACGCCTTTCTGGTgggg                                                                     >  1:2836488/1‑99 (MQ=255)

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

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

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

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

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

                                ttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCaa                                                                 >  1:4240193/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)

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

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

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

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

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

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

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

                                                      tCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGt                                           >  2:34056/1‑100 (MQ=255)

                                                      tCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGt                                           >  2:3912206/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)

                                                                   aaCCTGGATTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTgttgtt                               >  1:4319947/2‑99 (MQ=255)

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

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

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

                                                                      cTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCa                           >  2:2010250/1‑100 (MQ=255)

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

                                                                         gTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTg                        >  2:1013575/1‑100 (MQ=255)

                                                                          tttctttcTCACGATAAGAAACGCATCCGTTACTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCAGGTGGCGCTGTTGTTTTCACt                         >  1:3315567/1‑98 (MQ=255)

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

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

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

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

                                                                               ttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaa                  >  2:1814478/1‑100 (MQ=255)

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

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

                                                                                  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             >  1:247672/1‑100 (MQ=255)

                                                                                    aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCtt             <  2:2859010/100‑1 (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)

                                                                                               cGCATCTGTTTCTTATGCTCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCCGTTGTTTTCACTGTTTTAAATCTTAc           <  1:582676/91‑1 (MQ=255)

CGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCT  >  NC_000913/1517384‑1517578

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