BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,093,283	G→A	16.2%	G116D (GGC→GAC)	yqgE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,283	0	G	A	16.2%	400.4 / 56.6	173	G116D (GGC→GAC)	yqgE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (74/70); new base A (19/9); total (94/79)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 177 total aligned reads displayed.

GTTTATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGCTAA  >  NC_000913/3093185‑3093381

   tATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCAcc                                                                                                >  2:2041828/1‑100 (MQ=255)

     ttttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGa                                                                                              <  1:204623/100‑1 (MQ=255)

      tttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAAAACCGTAATGCCCGCTTCCCGCGCTGGGCTAGAAACGCTCGCCAACGa                                                                                              >  2:3797118/1‑99 (MQ=255)

      ttagCATATTCCCCCATCCGATGTGACTTCCAGCATTCCCATTTCAGGCAACACGGTAATGACCACTTACCGCTATGTGCTGGAAACGCTCGGCAc                                                                                                 <  2:3178972/93‑1 (MQ=255)

       ttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATa                                                                                            >  2:1319181/1‑100 (MQ=255)

       ttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATa                                                                                            >  2:3064436/1‑100 (MQ=255)

       ttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATATCAGACAACACGGTAATGACCACTTCCCGCGAGGTGCTGGAAACGCGCGGCACCGATa                                                                                            >  2:3270645/1‑100 (MQ=255)

          catgccccGCCCTGAAATTTTGCTTCCAGCATTCGCACTTCCGACACCAAGGTCACAACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATaaa                                                                                          <  2:3272114/93‑1 (MQ=255)

          cATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaac                                                                                         >  1:1074926/1‑100 (MQ=255)

          cATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACAACTTCCCCCGATGTGCTGGAAACGCTCGGCACCGATaaa                                                                                          >  1:1486944/1‑99 (MQ=255)

          cATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGAACACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAaac                                                                                         >  1:1020354/1‑100 (MQ=255)

              cTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaacaac                                                                                      >  2:1289165/1‑99 (MQ=255)

              cTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAAcc                                                                                     >  2:2076858/1‑100 (MQ=255)

               tcctccCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCTCGATGTGCTGGAAACGCTCGGCACCGATAaacaac                                                                                      <  2:953016/98‑1 (MQ=255)

                ccGCCCTCCAATTTTGCTTCCAGCATTCGTATTTCAGATAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCg                                                                                    <  1:6207/99‑1 (MQ=255)

                 cGCCCTCCATTTTTGCCTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGGGATGTGCTGGAAACGCTCGACACCGATAaacaac                                                                                      <  2:1968713/96‑1 (MQ=255)

                 cGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTc                                                                                  <  2:1107256/100‑1 (MQ=255)

                 cGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTc                                                                                  <  1:3534308/100‑1 (MQ=255)

                  gCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGAGATGTGCTGGAAACGCTCGACACCGATAAACAACCg                                                                                    >  2:1097001/1‑97 (MQ=255)

                   cccTCCAATTTTGCTTCCAGCATTCGCATTTCCGACAACACGGTAATAACCACTTCCCGCGATGTGCTGGAAACGCTCGGAACCGATAAACAACCGTCTg                                                                                >  1:2441550/1‑100 (MQ=255)

                    ccTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGa                                                                               <  1:1951992/100‑1 (MQ=255)

                     cTCCAATTTTTCTTTTAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGAc                                                                              >  2:3356493/1‑100 (MQ=255)

                     cTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGAc                                                                              >  1:2199674/1‑100 (MQ=255)

                     cTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGAc                                                                              >  1:347249/1‑100 (MQ=255)

                     cTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGAc                                                                              >  1:1765980/1‑100 (MQ=255)

                      tCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACg                                                                             >  2:3646716/1‑100 (MQ=255)

                       ccAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGt                                                                            <  1:943316/100‑1 (MQ=255)

                         aaTTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTTATGACCACTTCCCGCGATGTGCTGGAAACGCACAGCACAGATAAACAAACATCTCAc                                                                              >  1:1877359/1‑96 (MQ=255)

                          aTTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTAtt                                                                         <  1:1588332/100‑1 (MQ=255)

                           ttttGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtg                                                                        <  1:1166180/100‑1 (MQ=255)

                              tGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtggtg                                                                     <  2:1445311/100‑1 (MQ=255)

                              tGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtggtg                                                                     <  2:2686310/100‑1 (MQ=255)

                               gCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtggtgg                                                                    <  2:454554/100‑1 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtggtgg                                                                    >  2:300790/1‑99 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATtggtgg                                                                    >  2:1406764/1‑99 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGc                                                                   <  1:3289710/100‑1 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGc                                                                   <  1:3356531/100‑1 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGc                                                                   >  1:1328273/1‑100 (MQ=255)

                                cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGc                                                                   >  1:3763220/1‑100 (MQ=255)

                                 ttCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGAAAACGCTCGGCACCGATAAACAACCGACTGACGTATTGATGGct                                                                  >  1:2329899/1‑100 (MQ=255)

                                  tCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGctc                                                                 >  2:2223575/1‑100 (MQ=255)

                                    cAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGctct                                                                >  1:2937947/1‑99 (MQ=255)

                                      gCATTCGCATTTCAGACAACACGGTAATGACCATTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTggg                                                             <  1:702853/100‑1 (MQ=255)

                                        aTTCTCATTTCAGACAACACGGTACTGACCACTTCCCGCGATGTTTTGGATACGCTCGGCACCGATAAACAACCTTCTGACGTATTGGTGGTTCTGGGtt                                                           <  2:1346374/100‑1 (MQ=255)

                                        aTTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGtt                                                           >  1:1457036/1‑100 (MQ=255)

                                         ttCGCATTTCAGACAACACGGTAATGAACACTTCCTGAGATGTGCTGTAAAGGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTa                                                          >  2:3869475/1‑100 (MQ=255)

                                           cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTAt                                                         >  1:1547678/1‑99 (MQ=255)

                                           cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                        <  1:953188/100‑1 (MQ=255)

                                           cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                        <  2:2722892/100‑1 (MQ=255)

                                           cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                        >  2:428766/1‑100 (MQ=255)

                                           cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                        <  1:570431/100‑1 (MQ=255)

                                                ttCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcctcc                                                   <  2:3055554/100‑1 (MQ=255)

                                                 tCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcctcct                                                  >  2:3479931/1‑100 (MQ=255)

                                                 tCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcctcct                                                  >  1:101059/1‑100 (MQ=255)

                                                  cAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGCAAACGCTCGACACCGATAAACAACCGTCTGAAGAATTAGTGGCTCTGGGTTATGCCTCCTg                                                 >  1:1352761/1‑100 (MQ=255)

                                                      caacaCGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGgag                                             <  2:2230862/100‑1 (MQ=255)

                                                       aacaCGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGgaga                                            >  1:2354029/1‑100 (MQ=255)

                                                       aacaCGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGgag                                             <  1:882597/99‑1 (MQ=255)

                                                        acacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTAGTGGCTCTGGGTTATGCCTCCTGGGAGaa                                           >  2:172362/1‑100 (MQ=255)

                                                         cacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaaa                                          <  1:1993254/100‑1 (MQ=255)

                                                         cacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaaa                                          >  1:2484694/1‑100 (MQ=255)

                                                          acGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCTGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAg                                         <  1:1164270/100‑1 (MQ=255)

                                                           cGGTAATGACAACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAgg                                        <  2:108668/100‑1 (MQ=255)

                                                                aTGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGt                                       >  2:1266632/1‑96 (MQ=255)

                                                                aTGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAAc                                   >  1:3686117/1‑100 (MQ=255)

                                                                 tGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAAc                                   >  2:1013956/1‑99 (MQ=255)

                                                                  gACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTg                                 >  1:1711984/1‑100 (MQ=255)

                                                                  gACCACTTCCCGCGACGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACt                                  <  2:2287203/99‑1 (MQ=255)

                                                                      aCTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAAc                             <  2:760781/100‑1 (MQ=255)

                                                                       cTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACa                            >  1:1418939/1‑100 (MQ=255)

                                                                        ttCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACaa                           >  1:2998568/1‑100 (MQ=255)

                                                                        ttCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACaa                           >  1:2995243/1‑100 (MQ=255)

                                                                        ttCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACaa                           >  1:584036/1‑100 (MQ=255)

                                                                          cccGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGa                         >  1:2265705/1‑100 (MQ=255)

                                                                           ccGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGACTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGaa                        >  2:1484210/1‑100 (MQ=255)

                                                                           ccGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGACTCTGGGTTATGCATCCTGGGAGAAAGGTCAACTGGAACAAGaa                        >  1:1273294/1‑100 (MQ=255)

                                                                            cgcgATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGaaa                       <  2:2783910/100‑1 (MQ=255)

                                                                            cgcgATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGaaa                       >  1:3138690/1‑100 (MQ=255)

                                                                            cgcgATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGaaa                       <  1:2381808/100‑1 (MQ=255)

                                                                             gcgATGTGCTGGAAACGCTCGGCATCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAAt                      <  1:1607897/100‑1 (MQ=255)

                                                                             gcgATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAAt                      <  1:1231883/100‑1 (MQ=255)

                                                                              cgATGTGCTGTAATCTTTCTGCACTGATAAACAACCGTCGGACATATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGa                         <  1:3828264/96‑1 (MQ=255)

                                                                                 tgtgCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtctc                  <  1:3887407/100‑1 (MQ=255)

                                                                                 atctcTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtctc                  <  1:67309/96‑1 (MQ=255)

                                                                                   tgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGa                >  2:1974046/1‑100 (MQ=255)

                                                                                    gCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGAt               >  1:1918065/1‑100 (MQ=255)

                                                                                    gCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGAt               <  2:3702768/100‑1 (MQ=255)

                                                                                    gCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGAt               >  2:2878381/1‑100 (MQ=255)

                                                                                    gCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTAATGGTGGCTCTGGGTTATGACTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGa                >  2:3746973/1‑99 (MQ=255)

                                                                                     cTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtctcagat              >  2:67306/1‑96 (MQ=255)

                                                                                      tGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATaa             >  2:1179619/1‑100 (MQ=255)

                                                                                       ggAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAc            >  2:2555333/1‑100 (MQ=255)

                                                                                          aaCGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAcgcg         <  2:1553123/100‑1 (MQ=255)

                                                                                            cGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTg       >  1:2498881/1‑100 (MQ=255)

                                                                                             gCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTgg      <  1:2893638/100‑1 (MQ=255)

                                                                                              cTCGGCAGCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGc     <  2:3098022/100‑1 (MQ=255)

                                                                                              cTCGGCACCGATAAACAACCTTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGc     <  1:795459/100‑1 (MQ=255)

                                                                                               tCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGc     <  1:1691293/99‑1 (MQ=255)

                                                                                                cGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCTTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGCTa   <  1:1008372/100‑1 (MQ=255)

                                                                                                 ggCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGCTaa  >  2:1057160/1‑100 (MQ=255)

GTTTATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGCTAA  >  NC_000913/3093185‑3093381

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