BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,013,946	G→A	Q26Q (CAG→CAA)	fliH →	FliH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,013,946	0	G	A	97.6%	435.4 / ‑2.5	126	Q26Q (CAG→CAA)	fliH	FliH
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (58/65); total (60/66)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 127 total aligned reads displayed.

CGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAT  >  NC_000913/2013858‑2014035

 gagaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCaga                                                                                      >  2:2912725/1‑93 (MQ=255)

  agaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCagag                                                                                     <  1:2565273/93‑1 (MQ=255)

  agaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCagag                                                                                     <  1:1782114/93‑1 (MQ=255)

  agaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCagag                                                                                     <  1:886755/93‑1 (MQ=255)

   gaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGt                                                                                    >  1:4297165/1‑93 (MQ=255)

   gaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGt                                                                                    >  2:1968753/1‑93 (MQ=255)

    aTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGtt                                                                                   <  1:1804468/93‑1 (MQ=255)

    aTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGtt                                                                                   >  1:3127870/1‑93 (MQ=255)

     tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  <  1:1056219/93‑1 (MQ=255)

     tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  <  1:2678034/93‑1 (MQ=255)

     tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  >  1:3789162/1‑93 (MQ=255)

     tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  <  1:2203628/93‑1 (MQ=255)

     tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  >  2:3043042/1‑93 (MQ=255)

      ggTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtg                                                                                 <  1:752475/93‑1 (MQ=255)

      ggTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtg                                                                                 <  2:2546988/93‑1 (MQ=255)

       gTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgt                                                                                <  1:1995381/93‑1 (MQ=255)

        tAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgtg                                                                               >  1:1446339/1‑93 (MQ=255)

        tAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgtg                                                                               >  2:354006/1‑93 (MQ=255)

           ttGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGccc                                                                            >  1:28130/1‑93 (MQ=255)

                aGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGt                                                                       <  1:515274/93‑1 (MQ=255)

                aGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGt                                                                       <  1:4465279/93‑1 (MQ=255)

                   ggcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGa                                                                    >  2:1684490/1‑93 (MQ=255)

                    gcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAg                                                                   <  1:3574851/93‑1 (MQ=255)

                     cgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGc                                                                  >  1:4620225/1‑93 (MQ=255)

                     cgAGGATACCTATGTCTGATAAGCTGCCGTGGAAAAGCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGc                                                                  <  2:1566716/93‑1 (MQ=255)

                      gAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGcc                                                                 >  1:2917286/1‑93 (MQ=255)

                        ggATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCgg                                                               >  2:1974845/1‑93 (MQ=255)

                          aTACATATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaa                                                             <  1:204755/93‑1 (MQ=255)

                            aCCTATGTCTGATAATCTGCCGTGGAAACCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           <  2:4052855/93‑1 (MQ=255)

                            aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           <  1:1051065/93‑1 (MQ=255)

                            aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           >  2:2000146/1‑93 (MQ=255)

                            aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           <  1:1341042/93‑1 (MQ=255)

                            aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           >  2:3061793/1‑93 (MQ=255)

                             ccTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaagaa                                                          <  2:3194853/93‑1 (MQ=255)

                              cTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGaaa                                                         >  1:1844942/1‑93 (MQ=255)

                              cTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGaaa                                                         >  1:3768326/1‑93 (MQ=255)

                               tATGTCTGATAATCTGCCTTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        <  2:1167383/93‑1 (MQ=255)

                               tATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        <  1:626541/93‑1 (MQ=255)

                               tATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        <  1:1113037/93‑1 (MQ=255)

                               tATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        >  2:4074320/1‑93 (MQ=255)

                                aTGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAcc                                                       >  1:797335/1‑93 (MQ=255)

                                 tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtg                                                                                 >  2:3828515/1‑66 (MQ=255)

                                 tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACca                                                      >  1:510183/1‑93 (MQ=255)

                                 tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACca                                                      >  2:4194447/1‑93 (MQ=255)

                                   tCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatc                                                    <  2:3141965/93‑1 (MQ=255)

                                    cTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGGTTGTGCCCATAGTCGAGCCGGAAGAAACcatca                                                   >  1:1625321/1‑93 (MQ=255)

                                     tGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatcat                                                  <  2:1909/93‑1 (MQ=255)

                                       aTAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTg                                                <  1:371190/93‑1 (MQ=255)

                                        tgaTCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGa                                               >  2:4486842/3‑93 (MQ=255)

                                        tAATCTGCCGTGGAAAACCTGGACGCCGGAAGATCTCGCGCCACCACAAGCAGAGCTTGTGCCCATAGTCGAGCCGGCAGAAACCATCATTGa                                               >  2:1835824/1‑93 (MQ=255)

                                         aaTCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGaa                                              <  2:1127128/93‑1 (MQ=255)

                                            cTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAagag                                           <  2:1713557/93‑1 (MQ=255)

                                             tGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAgg                                          <  2:2811869/93‑1 (MQ=255)

                                              gCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGc                                         >  2:4189148/1‑93 (MQ=255)

                                              gCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGc                                         <  2:317850/93‑1 (MQ=255)

                                              gCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGc                                         >  1:2835822/1‑93 (MQ=255)

                                               ccGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCt                                        <  1:3368107/93‑1 (MQ=255)

                                                   ggAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAc                                    <  1:2306193/93‑1 (MQ=255)

                                                   ggAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAc                                    <  2:1476055/93‑1 (MQ=255)

                                                   ggAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAc                                    <  1:2097834/93‑1 (MQ=255)

                                                       aaCCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAACCCATCATTGAAGAGGCTGAACCCAg                                <  1:3917466/93‑1 (MQ=255)

                                                       aaCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAg                                >  1:3258247/1‑93 (MQ=255)

                                                        aCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGc                               >  2:1879367/1‑93 (MQ=255)

                                                         ccTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGcc                              <  1:4390168/93‑1 (MQ=255)

                                                         ccTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGcc                              <  2:2175249/93‑1 (MQ=255)

                                                          cTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCt                             <  2:3126584/93‑1 (MQ=255)

                                                          cTGGACGCAGGACGATCTCGCGCCACCACAAGCAGAGGTTGTGCCCATAGTCGAGCCGGAAGAAACGATAATTGAAGAGGCTGAGTCGAGCCt                             >  2:1831454/1‑93 (MQ=255)

                                                           tGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCtt                            >  2:925850/1‑93 (MQ=255)

                                                             gACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGa                          <  1:236844/93‑1 (MQ=255)

                                                               cGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagc                        >  2:2558884/1‑93 (MQ=255)

                                                                gCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagca                       <  2:2000373/93‑1 (MQ=255)

                                                                  cGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagcagc                     <  1:902652/93‑1 (MQ=255)

                                                                    gACGGTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGCGGCTGAACCCAGGCATGCGCAGCaa                   >  2:1754805/1‑93 (MQ=255)

                                                                    gACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCaa                   >  2:2942421/1‑93 (MQ=255)

                                                                    gACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCaa                   >  2:672550/1‑93 (MQ=255)

                                                                    gACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCaa                   >  1:2038458/1‑93 (MQ=255)

                                                                    gACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCaa                   >  2:4377927/1‑93 (MQ=255)

                                                                      cGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACt                 <  1:1837897/93‑1 (MQ=255)

                                                                      cGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACt                 <  1:2708459/93‑1 (MQ=255)

                                                                       gATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTg                <  1:4300879/93‑1 (MQ=255)

                                                                        aTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTgg               >  2:2843601/1‑93 (MQ=255)

                                                                        aTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTgg               <  1:1750222/93‑1 (MQ=255)

                                                                        aTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTgg               >  2:218340/1‑93 (MQ=255)

                                                                        aTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTgg               >  1:2656559/1‑93 (MQ=255)

                                                                        aTCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTgg               >  1:2114526/1‑93 (MQ=255)

                                                                          ctcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACTATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgcg             <  1:2513678/93‑1 (MQ=255)

                                                                          ctcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgcg             >  1:3287742/1‑93 (MQ=255)

                                                                           tcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgcgc            >  2:4669657/1‑93 (MQ=255)

                                                                           tcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgcgc            >  2:975207/1‑93 (MQ=255)

                                                                           tcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgcgc            <  1:3908963/93‑1 (MQ=255)

                                                                            cgcgCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCa           >  1:4243459/1‑93 (MQ=255)

                                                                              cgcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAAc         <  2:573505/93‑1 (MQ=255)

                                                                               gcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACt        >  2:1028562/1‑93 (MQ=255)

                                                                                    cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   <  2:3863961/93‑1 (MQ=255)

                                                                                    cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   <  1:637259/93‑1 (MQ=255)

                                                                                    cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   <  1:3919082/93‑1 (MQ=255)

                                                                                    cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   >  1:3362479/1‑93 (MQ=255)

                                                                                    cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   >  1:2129898/1‑93 (MQ=255)

                                                                                     acaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAt  >  2:2888650/1‑93 (MQ=255)

                                                                                     acaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAt  <  1:2205707/93‑1 (MQ=255)

CGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAT  >  NC_000913/2013858‑2014035

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