BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,678,204	G→A	P456L (CCG→CTG)	yhjJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,678,204	0	G	A	99.7%	1218.5 / ‑5.7	324	P456L (CCG→CTG)	yhjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (146/177); total (146/178)
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.84e-01

Only 100 of 328 total aligned reads displayed.

TCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTA  >  NC_000913/3678114‑3678293

 cATCCTTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCt                                                                                        <  1:3733458/93‑1 (MQ=255)

  aTCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTc                                                                                       >  2:814915/1‑93 (MQ=255)

     cGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGTCATGATTTGATCCCAGACCGCCTGCAAAGCTTTCATGTTAAATTCCAGCTCgac                                                                                    >  2:5018131/1‑91 (MQ=255)

     cGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGcc                                                                                    >  2:2805840/1‑93 (MQ=255)

     cGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGcc                                                                                    <  1:460642/93‑1 (MQ=255)

      gTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCt                                                                                   <  2:3346105/93‑1 (MQ=255)

      gTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCt                                                                                   <  1:6385247/93‑1 (MQ=255)

       tGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCtt                                                                                  >  2:4001862/1‑93 (MQ=255)

        tgCGACAGAGGTGGTCGCAGCGGCGTTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCttt                                                                                 <  1:1102677/92‑1 (MQ=255)

        ggCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCttt                                                                                 <  2:5377076/93‑1 (MQ=255)

        ggCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCttt                                                                                 <  1:6203105/93‑1 (MQ=255)

         gCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  2:2747575/93‑1 (MQ=255)

         gCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  1:2706930/93‑1 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGc                                                                             <  2:4201194/93‑1 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGc                                                                             >  2:3105766/1‑93 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGc                                                                             <  2:1546981/93‑1 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGc                                                                             >  1:4513289/1‑93 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGc                                                                             <  1:1037112/93‑1 (MQ=255)

             cAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCt                                                                            <  2:5872384/93‑1 (MQ=255)

              agagGTGGTCGCAGCGGCGGTAGATGGGGCCATGGTTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTg                                                                           >  2:216052/1‑93 (MQ=255)

              agagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCACCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGccg                                                                           <  1:6431362/93‑3 (MQ=255)

              agagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTg                                                                           >  2:204472/1‑93 (MQ=255)

              agagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTg                                                                           >  1:4975550/1‑93 (MQ=255)

               gagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgc                                                                          >  2:6201033/1‑93 (MQ=255)

                agGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgca                                                                         <  2:2106142/93‑1 (MQ=255)

                agGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgca                                                                         >  2:4724351/1‑93 (MQ=255)

                  gtggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgcagc                                                                       <  2:5607092/93‑1 (MQ=255)

                  gtggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgcagc                                                                       >  2:2857709/1‑93 (MQ=255)

                     gtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGt                                                                    <  2:6098368/93‑1 (MQ=255)

                     gtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGt                                                                    <  2:38652/93‑1 (MQ=255)

                      gcGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTa                                                                   <  2:1033543/92‑1 (MQ=255)

                       cGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTAt                                                                  <  1:4685840/93‑1 (MQ=255)

                        gCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTAtt                                                                 >  1:4398465/1‑93 (MQ=255)

                        gCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTAtt                                                                 <  1:1775474/93‑1 (MQ=255)

                        gCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTAtt                                                                 <  1:15879/93‑1 (MQ=255)

                          aGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTaa                                                               <  2:4358215/93‑1 (MQ=255)

                           gcggcgGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAAc                                                              <  2:4090983/93‑1 (MQ=255)

                           gcggcgGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAAc                                                              >  2:1708611/1‑93 (MQ=255)

                             ggcggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGc                                                            >  2:449514/1‑93 (MQ=255)

                              gcggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGcc                                                           >  2:2501996/1‑93 (MQ=255)

                               cggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCa                                                          <  2:6368435/93‑1 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        <  1:5935626/93‑1 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        <  2:1390115/93‑1 (MQ=255)

                                  tAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCatat                                                       >  1:5191855/1‑93 (MQ=255)

                                   aGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATc                                                      >  1:3467880/1‑93 (MQ=255)

                                   aGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATc                                                      >  2:1159232/1‑93 (MQ=255)

                                    gATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCa                                                     >  1:2297372/1‑93 (MQ=255)

                                     aTGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCAt                                                    >  2:1275057/1‑93 (MQ=255)

                                        gggCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCg                                                 <  1:4798219/93‑1 (MQ=255)

                                          gCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGAc                                               >  2:4842936/1‑93 (MQ=255)

                                            cATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACaa                                             <  1:891607/93‑1 (MQ=255)

                                            cATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACaa                                             <  1:6136063/93‑1 (MQ=255)

                                                aTTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGt                                         <  1:2969780/93‑1 (MQ=255)

                                                aTTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGt                                         <  2:1788012/93‑1 (MQ=255)

                                                  ttGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTg                                       <  2:1178164/93‑1 (MQ=255)

                                                   tGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGa                                      <  2:5305282/93‑1 (MQ=255)

                                                   tGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGa                                      >  1:3134415/1‑93 (MQ=255)

                                                    gATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     >  1:100187/1‑93 (MQ=255)

                                                    gATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     <  1:2850352/93‑1 (MQ=255)

                                                    gATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCCGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     <  2:3573873/93‑1 (MQ=255)

                                                     aTCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACg                                    >  2:6538122/1‑93 (MQ=255)

                                                           gACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATc                              >  1:6077413/1‑93 (MQ=255)

                                                           gACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATc                              >  2:6715367/1‑93 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  1:5643298/93‑1 (MQ=255)

                                                               gCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGa                          >  2:1964511/1‑93 (MQ=255)

                                                               gCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGa                          <  2:4792989/93‑1 (MQ=255)

                                                                ccTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAt                         >  1:1806148/1‑93 (MQ=255)

                                                                 cTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACGACTGTTGACGCAGATCCTGAtt                        >  2:3928109/1‑93 (MQ=255)

                                                                 cTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAtt                        <  1:215752/93‑1 (MQ=255)

                                                                 cTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAtt                        <  1:356375/93‑1 (MQ=255)

                                                                  tGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAttt                       <  2:2918900/93‑1 (MQ=255)

                                                                   gCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTa                      >  2:3154318/1‑93 (MQ=255)

                                                                    cAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTaa                     <  1:2653434/93‑1 (MQ=255)

                                                                     aaCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAAc                    <  1:1407455/93‑1 (MQ=255)

                                                                     aaCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAAc                    <  1:1740102/93‑1 (MQ=255)

                                                                     aaCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAAc                    >  1:5760342/1‑93 (MQ=255)

                                                                     aaCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGGTGACGCAGATCCTGATTTAAc                    >  1:3277179/1‑93 (MQ=255)

                                                                      aCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACa                   <  1:3877441/93‑1 (MQ=255)

                                                                      aCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACa                   >  1:4252413/1‑93 (MQ=255)

                                                                         cTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtct                >  2:5133543/1‑93 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               >  1:2544620/1‑93 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               <  1:2402903/93‑1 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               >  1:6487615/1‑93 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               >  1:4071846/1‑93 (MQ=255)

                                                                            tCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTa             <  2:3775025/93‑1 (MQ=255)

                                                                             cATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTAc            <  1:5440193/93‑1 (MQ=255)

                                                                              aTGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTAcc           >  2:59969/1‑93 (MQ=255)

                                                                               tGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCg          >  1:1954449/1‑93 (MQ=255)

                                                                               tGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCg          >  1:6603375/1‑93 (MQ=255)

                                                                                gTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGt         >  2:3213393/1‑93 (MQ=255)

                                                                                gTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGt         >  1:5138798/1‑93 (MQ=255)

                                                                                 ttAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTc        <  2:1831674/93‑1 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCa       <  2:232764/93‑1 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCa       <  2:2198667/93‑1 (MQ=255)

                                                                                   aaaTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAg      <  1:4359144/93‑1 (MQ=255)

                                                                                    aaTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAgg     >  1:4190422/1‑93 (MQ=255)

                                                                                     aTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTCACATCTCTACCGTCAGGc    >  2:5882608/1‑93 (MQ=255)

                                                                                     aTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGc    <  1:6290129/93‑1 (MQ=255)

                                                                                      ttCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCt   >  1:2507057/1‑93 (MQ=255)

                                                                                       tCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTa  <  2:3405100/93‑1 (MQ=255)

TCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTA  >  NC_000913/3678114‑3678293

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