BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,387,305	G→A	P124P (CCG→CCA)	yfbM →	DUF1877 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,387,305	0	G	A	96.6%	678.4 / 5.9	205	P124P (CCG→CCA)	yfbM	DUF1877 family protein
Reads supporting (aligned to +/- strand): ref base G (4/3); new base A (100/97); total (105/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 208 total aligned reads displayed.

AGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCGA  >  NC_000913/2387208‑2387401

aGTTAGCCGCGACAATCGACCGACTGGAATCACTTGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCagg                                                                                                <  2:1705003/100‑1 (MQ=255)

aGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCagg                                                                                                <  2:2844042/100‑1 (MQ=255)

 gTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGa                                                                                               >  1:391239/1‑100 (MQ=255)

   tAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGt                                                                                             >  1:4966081/1‑100 (MQ=255)

    aGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTg                                                                                            >  2:2471534/1‑100 (MQ=255)

    aGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTg                                                                                            >  2:3761736/1‑100 (MQ=255)

      ccGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGAc                                                                                          <  2:5184575/100‑1 (MQ=255)

        gcgACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCt                                                                                        <  2:3590759/100‑1 (MQ=255)

         cgACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCtt                                                                                       <  1:3539183/100‑1 (MQ=255)

             aaTCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCa                                                                                   <  2:4561346/100‑1 (MQ=255)

             aaTCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCa                                                                                   <  2:2568480/100‑1 (MQ=255)

             aaTCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCa                                                                                   <  2:3556289/100‑1 (MQ=255)

               tCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTGAGCGACTCAATGAAATGGAGATTTAACCAGGAGTGACCTTTTCaga                                                                                 <  2:3249999/100‑1 (MQ=255)

               tCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCaga                                                                                 <  2:1099887/100‑1 (MQ=255)

                cgaccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCagaa                                                                                >  2:865389/1‑100 (MQ=255)

                 gaccgacTGGAATCACTGGATCGAAATGAAATGCGTAAACAGTTGTCCATTAAGCGACTCAATGAATTGTATATTTATCCAGGAGTGACCTTTTCagaag                                                                               <  2:1907394/100‑1 (MQ=255)

                  accgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCagaaga                                                                              <  1:4759010/100‑1 (MQ=255)

                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                             >  1:3242560/1‑100 (MQ=255)

                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                             >  1:2192608/1‑100 (MQ=255)

                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                             <  1:2179059/100‑1 (MQ=255)

                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                             <  2:487747/100‑1 (MQ=255)

                    cgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGt                                                                            <  1:2724595/100‑1 (MQ=255)

                      acTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTa                                                                          >  1:5494574/1‑100 (MQ=255)

                      acTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTa                                                                          <  2:637785/100‑1 (MQ=255)

                       cTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTag                                                                         >  2:545691/1‑100 (MQ=255)

                        tGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTaga                                                                        <  2:2973276/100‑1 (MQ=255)

                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                       >  1:3764076/1‑100 (MQ=255)

                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                       <  1:2661980/100‑1 (MQ=255)

                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                       <  1:4459172/100‑1 (MQ=255)

                          gAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgg                                                                      >  2:49036/1‑100 (MQ=255)

                          gAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgg                                                                      <  1:5159599/100‑1 (MQ=255)

                          gAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgg                                                                      >  1:5262848/1‑100 (MQ=255)

                          gAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTGTCAGATGAGTTAGAgg                                                                      <  1:4872538/100‑1 (MQ=255)

                           aaTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAggg                                                                     >  2:2391799/1‑100 (MQ=255)

                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                    >  1:921737/1‑100 (MQ=255)

                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                    <  1:4606508/100‑1 (MQ=255)

                              cACTGGATCGAAATGAACTGCGTAAACAGTCTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCa                                                                  <  2:1276657/100‑1 (MQ=255)

                              aaCTGGGTCGAAATAAACGGCGTAATCAGTTTTCCATTAAGCGACTGAATGAAATTGGGATGTATCCAGGAGTGGCCTTTTCAGAAGAGTTAGAGAGGCa                                                                  <  2:4654035/99‑1 (MQ=255)

                               aCTGGATCGAAATGAACTGCGTAAACAGTTCTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCaa                                                                 >  2:5111337/1‑100 (MQ=255)

                                cTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAc                                                                >  1:5042829/1‑100 (MQ=255)

                                cTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAc                                                                >  2:1454171/1‑100 (MQ=255)

                                  ggATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAActc                                                              >  1:3658506/1‑100 (MQ=255)

                                  ggATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAActc                                                              >  1:5080247/1‑100 (MQ=255)

                                   gATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAActct                                                             >  1:3784177/1‑100 (MQ=255)

                                    aTCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCtt                                                            >  2:3499934/1‑100 (MQ=255)

                                       gAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGGCCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGc                                                         >  2:618012/1‑100 (MQ=255)

                                       gAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGc                                                         <  2:3877283/100‑1 (MQ=255)

                                           tGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCg                                                     <  1:3822907/100‑1 (MQ=255)

                                              aCTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATc                                                  >  1:5415077/1‑100 (MQ=255)

                                               cTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCa                                                 <  1:166130/100‑1 (MQ=255)

                                               cTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCa                                                 >  1:4231592/1‑100 (MQ=255)

                                               cTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCa                                                 >  1:21400/1‑100 (MQ=255)

                                                 gCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATg                                               >  2:1569774/1‑100 (MQ=255)

                                                 gCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGCTCATg                                               >  2:3711892/1‑100 (MQ=255)

                                                   gTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCt                                             <  1:3532491/100‑1 (MQ=255)

                                                   gTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCt                                             <  1:663017/100‑1 (MQ=255)

                                                      aaCAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGa                                          >  1:3132805/1‑100 (MQ=255)

                                                       aCAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGat                                         >  2:2669844/1‑100 (MQ=255)

                                                       aCAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGat                                         <  2:479596/100‑1 (MQ=255)

                                                         aGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGatat                                       <  2:3448534/100‑1 (MQ=255)

                                                          gTTTTACATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAATAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATg                                      <  1:5284953/100‑1 (MQ=255)

                                                              tCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGaaa                                  <  2:3661098/100‑1 (MQ=255)

                                                               ccATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGaaaa                                 >  2:1884005/1‑100 (MQ=255)

                                                               ccATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGaaaa                                 <  1:755786/100‑1 (MQ=255)

                                                                cATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGTGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                                <  2:5368256/100‑1 (MQ=255)

                                                                cATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                                >  2:2620993/1‑100 (MQ=255)

                                                                cATAAAGCGACTCAATGAAATGGAGATGTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                                <  2:5339418/100‑1 (MQ=255)

                                                                  ttAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCt                              <  2:3385420/100‑1 (MQ=255)

                                                                    aaGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGa                            <  2:2247489/100‑1 (MQ=255)

                                                                     aGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGAt                           >  1:1589101/1‑100 (MQ=255)

                                                                     aGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGAt                           >  1:3839533/1‑100 (MQ=255)

                                                                      gCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGAtt                          >  1:2080790/1‑100 (MQ=255)

                                                                      gCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGAtt                          <  2:2163605/100‑1 (MQ=255)

                                                                       cGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTa                         >  1:3974510/1‑100 (MQ=255)

                                                                        gACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAg                        >  1:1037529/1‑100 (MQ=255)

                                                                         aCTCAATGAAATGGAGATTTATCCAGGAGTGGCCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgc                       <  1:625453/100‑1 (MQ=255)

                                                                         aCTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgc                       >  2:640514/1‑100 (MQ=255)

                                                                          cTCACTGAAATGGAGATTTATCCAGGAGTGACCTTTTCGGACGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGACATGGAATAGCTGATTAggg                      >  2:1103442/1‑98 (MQ=255)

                                                                          cTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgcg                      >  1:1726568/1‑100 (MQ=255)

                                                                          cTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgcg                      >  2:1150051/1‑100 (MQ=255)

                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                    >  1:531837/1‑100 (MQ=255)

                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                    <  1:1999288/100‑1 (MQ=255)

                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                    >  1:2332627/1‑100 (MQ=255)

                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                    >  1:4330891/1‑100 (MQ=255)

                                                                             aaTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGt                   >  1:818104/1‑100 (MQ=255)

                                                                             aaTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGt                   >  1:982097/1‑100 (MQ=255)

                                                                              aTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTa                  <  2:4297107/100‑1 (MQ=255)

                                                                               tGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTAt                 <  1:1415670/100‑1 (MQ=255)

                                                                                  aaTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGa              >  1:2032880/1‑100 (MQ=255)

                                                                                   aTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGAc             <  1:659792/100‑1 (MQ=255)

                                                                                   aTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGAc             >  1:2506512/1‑100 (MQ=255)

                                                                                    tGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACg            <  1:5346214/100‑1 (MQ=255)

                                                                                    tGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACg            >  1:127121/1‑100 (MQ=255)

                                                                                     ggAGATTTATCCAGGAGTGACCTTTTCACAAGAGTTCGAGGGGCAACTCTTTGCCCCGATCATGCTTGCTATGGCAAAGCTGATTAGCGCGTATCGACGt           <  2:2763192/100‑1 (MQ=255)

                                                                                       agaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAt         <  2:4005369/100‑1 (MQ=255)

                                                                                        gaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTAGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtg        >  1:4933303/1‑100 (MQ=255)

                                                                                           ttATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCAATCATGCTTGATATGGAAAAGATGATTAGCGCGTATCGACGTAtgctg     >  2:4981855/1‑100 (MQ=255)

                                                                                            tATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtgctgc    >  1:3169185/1‑100 (MQ=255)

                                                                                            tATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtgctgc    >  2:1463036/1‑100 (MQ=255)

                                                                                             aTCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCg   >  2:1792186/1‑100 (MQ=255)

AGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCGA  >  NC_000913/2387208‑2387401

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