BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,063	A→G	D129G (GAC→GGC)	ygfR →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,063	0	A	G	99.7%	1057.5 / ‑5.2	334	D129G (GAC→GGC)	ygfR	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (161/171); total (162/172)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 345 total aligned reads displayed.

CTACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTA  >  NC_000913/3025973‑3026152

cgaCTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCt                                                                                         >  2:7028942/3‑93 (MQ=255)

cTACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCt                                                                                         >  2:2770684/1‑93 (MQ=255)

  aCTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCa                                                                                       <  2:3696683/93‑1 (MQ=255)

    tGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGt                                                                                     >  2:4213104/1‑93 (MQ=255)

     ggCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGGCCAGTGGCTCAGTa                                                                                    >  1:1813088/1‑93 (MQ=255)

      gCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTaa                                                                                   <  2:50939/93‑1 (MQ=255)

        aCCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGt                                                                                 <  2:776527/93‑1 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGAGCAGTGGCTCAGTAAGtt                                                                                >  2:7814599/1‑93 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGtt                                                                                >  1:6217078/1‑93 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGtt                                                                                >  2:6632071/1‑93 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGtt                                                                                <  1:3969715/93‑1 (MQ=255)

           aTCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTTACCAGTGGCTCAGTAAGTTCa                                                                              <  2:3781606/93‑1 (MQ=255)

           aTCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCa                                                                              <  1:809558/93‑1 (MQ=255)

            tCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCaa                                                                             >  1:2074912/1‑93 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           <  2:6097520/93‑1 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           <  1:5653342/93‑1 (MQ=255)

               gTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGCTCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAatat                                                                          <  2:3098591/93‑1 (MQ=255)

               gTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAatat                                                                          <  1:2787167/93‑1 (MQ=255)

                tGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATAtt                                                                         >  1:1661197/1‑93 (MQ=255)

                ggCCGGCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATAtt                                                                         <  1:6100702/92‑1 (MQ=255)

                 gCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtc                                                                        <  2:1762767/93‑1 (MQ=255)

                 gCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtc                                                                        <  1:193361/93‑1 (MQ=255)

                 gCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATAAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtc                                                                        >  2:3629025/1‑93 (MQ=255)

                   cGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtctc                                                                      <  1:3223620/93‑1 (MQ=255)

                   cGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtctc                                                                      <  1:6404511/93‑1 (MQ=255)

                    gTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTcc                                                                     <  1:146749/93‑1 (MQ=255)

                     tCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCg                                                                    <  2:2642313/93‑1 (MQ=255)

                     tCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCg                                                                    >  2:2065377/1‑93 (MQ=255)

                     tCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCg                                                                    >  1:5042514/1‑93 (MQ=255)

                      cGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCgg                                                                   <  2:4857052/93‑1 (MQ=255)

                       gCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGc                                                                  >  1:542511/1‑93 (MQ=255)

                         cGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCt                                                                <  2:4436950/93‑1 (MQ=255)

                          ggTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtg                                                               >  2:2211471/1‑93 (MQ=255)

                          ggTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtg                                                               >  1:250899/1‑93 (MQ=255)

                          ggTCAGGCGAATGTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtg                                                               <  1:4905360/93‑1 (MQ=255)

                           gTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtgg                                                              >  2:3162087/1‑93 (MQ=255)

                            tCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggt                                                             <  1:1569822/93‑1 (MQ=255)

                              aGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggtgg                                                           <  1:4814697/93‑1 (MQ=255)

                              aGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggtgg                                                           <  1:7364594/93‑1 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          <  2:3911392/93‑1 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          <  1:2933974/93‑1 (MQ=255)

                                  gAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGgcgc                                                       >  2:5210308/1‑93 (MQ=255)

                                    aTTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAg                                                     <  2:5788588/93‑1 (MQ=255)

                                    aTTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAg                                                     >  1:2159440/1‑93 (MQ=255)

                                     tttGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgc                                                    <  2:1529813/93‑1 (MQ=255)

                                      ttgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcg                                                   <  1:3368352/93‑1 (MQ=255)

                                      ttgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcg                                                   <  2:6383008/93‑1 (MQ=255)

                                       tgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcgc                                                  >  2:649764/1‑93 (MQ=255)

                                       tgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcgc                                                  >  2:2613804/1‑93 (MQ=255)

                                        gttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGcc                                                 >  1:1382620/1‑93 (MQ=255)

                                         ttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGccc                                                >  1:6721207/1‑93 (MQ=255)

                                          tgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCg                                               <  1:5921503/93‑1 (MQ=255)

                                            gATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCa                                             <  2:6441045/93‑1 (MQ=255)

                                            gATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCa                                             >  2:4008723/1‑93 (MQ=255)

                                              tAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGc                                           >  2:7317047/1‑93 (MQ=255)

                                              tAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGc                                           <  2:7648970/93‑1 (MQ=255)

                                              tAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGc                                           >  2:4763405/1‑93 (MQ=255)

                                                     aaCCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTtat                                    >  2:3332766/1‑93 (MQ=255)

                                                     aaCCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGGTtat                                    >  2:3518205/1‑93 (MQ=255)

                                                      aCCAGAGCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTtata                                   <  1:1251841/93‑1 (MQ=255)

                                                       ccAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTatat                                  <  2:845374/93‑1 (MQ=255)

                                                       ccAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTatat                                  >  1:1054123/1‑93 (MQ=255)

                                                       ccAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTatat                                  <  2:1091679/93‑1 (MQ=255)

                                                        cAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATc                                 >  2:6773421/1‑93 (MQ=255)

                                                         aGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCg                                <  2:4684072/93‑1 (MQ=255)

                                                          gATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGa                               >  2:7418257/1‑93 (MQ=255)

                                                           atcatcAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGaa                              <  2:2986159/93‑1 (MQ=255)

                                                           atcatcAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGaa                              <  1:6615350/93‑1 (MQ=255)

                                                            tcatcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAAt                             <  2:4636982/93‑1 (MQ=255)

                                                             catcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATc                            <  2:6311970/93‑1 (MQ=255)

                                                             catcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATc                            <  2:7184856/93‑1 (MQ=255)

                                                              atcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCa                           <  1:7685464/93‑1 (MQ=255)

                                                              atcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCa                           <  1:3113374/93‑1 (MQ=255)

                                                              atcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCa                           <  2:197044/93‑1 (MQ=255)

                                                               tcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAg                          >  2:1593588/1‑93 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         >  1:7127715/1‑93 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         <  1:7791224/93‑1 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         <  2:6855055/93‑1 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTAAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         >  1:7574850/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  2:2834508/1‑93 (MQ=255)

                                                                   cggcggCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGc                      <  1:6297044/93‑1 (MQ=255)

                                                                    ggcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCa                     <  2:1808457/93‑1 (MQ=255)

                                                                     gcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAg                    <  2:5484639/93‑1 (MQ=255)

                                                                       ggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGa                  <  1:4517363/93‑1 (MQ=255)

                                                                           aaGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccg              <  1:5599942/93‑1 (MQ=255)

                                                                            aGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgc             <  1:4729820/93‑1 (MQ=255)

                                                                            aGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgc             >  1:3750958/1‑93 (MQ=255)

                                                                              cccTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccggcg           >  1:173336/1‑93 (MQ=255)

                                                                                cTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgcc         >  2:6248818/1‑93 (MQ=255)

                                                                                cTGACCAGTGGCTCAGTAAGTTCAATATGCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgcc         <  1:4304739/93‑1 (MQ=255)

                                                                                 tGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgccg        <  1:694093/93‑1 (MQ=255)

                                                                                 tGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGCTTATATCGAATCAGCGGCAGGAAccgccgccg        >  2:1145275/1‑93 (MQ=255)

                                                                                  gACCAGTGGCGCAGTAAGTTCAATATGCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAGCCGCCGCCgg       <  2:1289863/93‑1 (MQ=255)

                                                                                   aCCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCggg      <  1:4112180/93‑1 (MQ=255)

                                                                                     cAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCggggg    <  1:4400849/93‑1 (MQ=255)

                                                                                      aGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCgggggg   >  2:5570221/1‑92 (MQ=255)

                                                                                      aGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCgggggg   >  1:1700554/1‑92 (MQ=255)

                                                                                       gTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTTGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTa  >  1:4050158/1‑93 (MQ=255)

                                                                                       gTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTa  >  1:2316575/1‑93 (MQ=255)

                                                                                       gTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTa  >  2:2407045/1‑93 (MQ=255)

CTACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTA  >  NC_000913/3025973‑3026152

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