BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,658,003	C→G	15.2%	intergenic (+436/‑363)	hdeD → / → gadE	acid‑resistance membrane protein/gad regulon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,658,003	0	C	G	15.2%	368.1 / 53.4	132	intergenic (+436/‑363)	hdeD/gadE	acid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand): ref base C (62/50); new base G (12/8); total (74/58)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 138 total aligned reads displayed.

CATAAATATCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGAT  >  NC_000913/3657904‑3658102

cATAAATATCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATcc                                                                                                     >  2:21261/1‑100 (MQ=255)

cATAAATATCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATcc                                                                                                     >  2:497395/1‑100 (MQ=255)

 aTAAATATCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCa                                                                                                    >  2:1491219/1‑100 (MQ=255)

     atatCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTg                                                                                                <  1:508284/100‑1 (MQ=255)

     atatCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTg                                                                                                <  2:2351664/100‑1 (MQ=255)

        tCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCa                                                                                             >  2:599673/1‑100 (MQ=255)

         ccGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAg                                                                                            <  2:731646/100‑1 (MQ=255)

         ccGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAg                                                                                            <  2:1022918/100‑1 (MQ=255)

          cGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGc                                                                                           <  2:1492814/100‑1 (MQ=255)

          cGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGc                                                                                           <  2:1113698/100‑1 (MQ=255)

          cGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAg                                                                                            >  1:646155/1‑99 (MQ=255)

           gtgtTTCCAGACGCTATATGAAAACCTGAAGACATGAATGCTTTAGTTGCTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGc                                                                                           <  1:2040296/99‑1 (MQ=255)

           gtgtCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAg                                                                                            <  2:1337987/98‑1 (MQ=255)

            tgtCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCtt                                                                                         <  1:2069412/100‑1 (MQ=255)

              tctcCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGCtttc                                                                                       <  2:281066/100‑3 (MQ=39)

               ctcCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAg                                                                                      <  1:1301452/100‑1 (MQ=255)

                tcCAGACGGTATATAAAAACCTGAAGACATGAATGCGGTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAg                                                                                      <  2:2501148/99‑1 (MQ=255)

                tcCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGt                                                                                     >  2:1850032/1‑100 (MQ=255)

                tcCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGt                                                                                     >  2:1290819/1‑100 (MQ=255)

                 ccAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTc                                                                                    >  2:2327518/1‑100 (MQ=255)

                  cAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGCtttctttt                                                                                   <  1:724350/100‑7 (MQ=38)

                  cAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCg                                                                                   >  2:1500611/1‑100 (MQ=255)

                   aGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGCtttctttta                                                                                  <  2:1233948/100‑8 (MQ=37)

                    gACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGaa                                                                                 <  2:1651492/100‑1 (MQ=255)

                    gACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGaa                                                                                 >  2:2594908/1‑100 (MQ=255)

                        cTATATAAAAACCTGAAGACATGAATACGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACaa                                                                             >  1:943715/1‑100 (MQ=255)

                         tatataAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAg                                                                            >  1:635494/1‑100 (MQ=255)

                         tatataAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAg                                                                            >  2:1978365/1‑100 (MQ=255)

                         tatataAAAACCTGAAGACATGAATACGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATTGGATTGCCAGCtttcttttaac                                                                                <  1:1079578/96‑10 (MQ=25)

                           tataAAAACCTGAAGATATGAATACGTTATTTACTCAGGTAATTTAAATGCGTTAAAAGAAAGCTGGAAATCCAATTGCCAGCTTAAGTCGAAACAAGga                                                                          <  1:2251035/100‑1 (MQ=255)

                           tataAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGga                                                                          <  2:117700/100‑1 (MQ=255)

                            ataAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGgag                                                                         >  1:2372856/1‑100 (MQ=255)

                            ataAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGgag                                                                         >  1:1937495/1‑100 (MQ=255)

                             taAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGgaga                                                                        >  2:1328082/1‑100 (MQ=255)

                                 aaCCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCg                                                                    <  1:1024183/100‑1 (MQ=255)

                                   ccTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGat                                                                  <  2:2354476/100‑1 (MQ=255)

                                    cTGAAGACATGAATGCGTTATTTACTCAGGTAATTTTAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGata                                                                 >  1:1079969/1‑100 (MQ=255)

                                     tGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGatat                                                                <  2:626060/100‑1 (MQ=255)

                                     tGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGatat                                                                <  2:34767/100‑1 (MQ=255)

                                       aaGACATGAATGCGTTCTTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATAttt                                                              <  1:1513427/100‑1 (MQ=255)

                                       aaGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCAAGCTTAAGTCGAAACAAGGAGACTCGATAttt                                                              <  1:404824/100‑1 (MQ=255)

                                        aGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTa                                                             <  2:8039/100‑1 (MQ=255)

                                           cATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAAt                                                          >  1:139452/1‑100 (MQ=255)

                                             tGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCg                                                        >  2:929255/1‑100 (MQ=255)

                                              gAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCgg                                                       <  1:660158/100‑1 (MQ=255)

                                                 tGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGAtt                                                    >  1:1334969/1‑100 (MQ=255)

                                                    gTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACa                                                 >  2:2118925/1‑100 (MQ=255)

                                                       aTTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACAttt                                              <  2:1015747/100‑1 (MQ=255)

                                                        tttACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACAtttt                                             >  2:190560/1‑100 (MQ=255)

                                                        tttACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACAtttt                                             <  1:997672/100‑1 (MQ=255)

                                                          tACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTaa                                           >  1:194478/1‑100 (MQ=255)

                                                          tACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTaa                                           >  2:43599/1‑100 (MQ=255)

                                                          tACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTaa                                           >  1:649634/1‑100 (MQ=255)

                                                           aCTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAAc                                          <  2:2216585/100‑1 (MQ=255)

                                                           aCTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAAc                                          <  2:2216240/100‑1 (MQ=255)

                                                              cAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACttt                                       >  1:2480864/1‑100 (MQ=255)

                                                              cAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACttt                                       >  2:766146/1‑100 (MQ=255)

                                                               aGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACtt                                        >  1:2487388/1‑98 (MQ=255)

                                                                ngTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAg                                     >  2:811646/2‑100 (MQ=255)

                                                                ggTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAg                                     <  1:1000384/100‑1 (MQ=255)

                                                                ggTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAg                                     >  1:889816/1‑100 (MQ=255)

                                                                 gTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGt                                    >  2:599972/1‑100 (MQ=255)

                                                                  tAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTa                                   >  1:1786252/1‑100 (MQ=255)

                                                                     tttCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTACCGTTAGTAata                                >  1:1937633/1‑100 (MQ=255)

                                                                     tttCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTaag                                 >  1:1568063/1‑98 (MQ=255)

                                                                     tttCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAata                                >  2:696516/1‑100 (MQ=255)

                                                                       tCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATAtt                              >  1:2224905/1‑100 (MQ=255)

                                                                        cAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATAttc                             <  2:2370532/100‑1 (MQ=255)

                                                                           tGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATAttcttc                          <  2:139453/100‑1 (MQ=255)

                                                                             ctttAATAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGCAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCag                        <  1:1250768/98‑1 (MQ=255)

                                                                              gTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCaga                       >  1:374026/1‑100 (MQ=255)

                                                                              gTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTAAAATCGTATTACATTTTAACTTTAGTAATATTCTTCaga                       >  1:1817538/1‑100 (MQ=255)

                                                                               ttAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCagag                      <  1:1460017/100‑1 (MQ=255)

                                                                               ttAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCagag                      >  1:1760717/1‑100 (MQ=255)

                                                                                 aaaaGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGAt                    >  1:1049882/1‑100 (MQ=255)

                                                                                  gacttaagCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATc                   >  1:2252268/6‑100 (MQ=38)

                                                                                  aaagaaagCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAgaga                     >  1:1551487/1‑98 (MQ=255)

                                                                                  aaagaaagCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATc                   <  2:1871388/100‑1 (MQ=255)

                                                                                    agaaagCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATcac                 >  2:242340/1‑100 (MQ=255)

                                                                                    agaaagCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATcac                 >  2:2472958/1‑100 (MQ=255)

                                                                                     gaaagCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACGTTAGTAATATTCTTCAGAGATcaca                <  1:342768/100‑1 (MQ=255)

                                                                                        agCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAAc             >  2:1591081/1‑100 (MQ=255)

                                                                                        agCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAAc             <  2:1937505/100‑1 (MQ=255)

                                                                                         gCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCCGAGATCACAAACt            >  1:1271746/1‑100 (MQ=255)

                                                                                         gCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACt            >  1:1528520/1‑100 (MQ=255)

                                                                                         gCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACt            >  2:1455794/1‑100 (MQ=255)

                                                                                         gCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACt            >  1:1648726/1‑100 (MQ=255)

                                                                                         gCTGGCAATTGGATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACt            >  2:530796/1‑100 (MQ=255)

                                                                                         gCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACt            <  2:2053441/100‑1 (MQ=255)

                                                                                           tGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTgg          <  2:1191624/100‑1 (MQ=255)

                                                                                           tGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTgg          <  1:235967/100‑1 (MQ=255)

                                                                                              cAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTa       >  2:1146122/1‑100 (MQ=255)

                                                                                              cAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTa       >  2:2379471/1‑100 (MQ=255)

                                                                                                aTCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTAtt     <  1:798931/100‑1 (MQ=255)

                                                                                                aTCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTAtt     >  2:1956843/1‑100 (MQ=255)

                                                                                                 tCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTg    >  2:1605211/1‑100 (MQ=255)

                                                                                                  ccAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGa   >  2:1862420/1‑100 (MQ=255)

                                                                                                   cAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGAt  <  1:2393321/100‑1 (MQ=255)

                                                                                                   cAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGAt  >  2:918345/1‑100 (MQ=255)

                                                                                                   cAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGAt  >  1:556244/1‑100 (MQ=255)

CATAAATATCCGTGTCTCCAGACGCTATATAAAAACCTGAAGACATGAATGCGTTATTTACTCAGGTAATTTCAATGCGTTAAAAGAAAGCTGGCAATCCAATTGCCAGCTTAAGTCGAAACAAGGAGACTCGATATTTAAATCGGATTACATTTTAACTTTAGTAATATTCTTCAGAGATCACAAACTGGTTATTGAT  >  NC_000913/3657904‑3658102

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