BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,649	A→G	V201A (GTG→GCG)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,649	0	A	G	99.7%	1087.7 / ‑5.5	347	V201A (GTG→GCG)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (168/178); total (169/178)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 363 total aligned reads displayed.

TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738

tCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTTTCcgcg                                                                                         >  1:3201634/1‑93 (MQ=255)

  gTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCcgcgcg                                                                                       >  2:1946693/1‑93 (MQ=255)

  gTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCcgcgcg                                                                                       <  1:2199673/93‑1 (MQ=255)

   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      >  2:4945048/1‑93 (MQ=255)

        cATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaaa                                                                                 >  1:6005022/1‑93 (MQ=255)

         aTCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaaaa                                                                                >  2:4809344/1‑93 (MQ=255)

         aTCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGGGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTaaaa                                                                                >  2:1356508/1‑93 (MQ=255)

          tCTGCACAGTTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCACCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               <  2:5769366/93‑1 (MQ=255)

          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               >  1:6112948/1‑93 (MQ=255)

          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               >  1:2395925/1‑93 (MQ=255)

            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  1:6477862/93‑1 (MQ=255)

            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  1:3809551/93‑1 (MQ=255)

            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  2:5187662/93‑1 (MQ=255)

             gCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacgac                                                                            <  1:5716856/93‑1 (MQ=255)

               acaGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgc                                                                          <  2:5238513/93‑1 (MQ=255)

                caGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgcg                                                                         <  2:7269927/93‑1 (MQ=255)

                 aGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGa                                                                        <  1:4558505/93‑1 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  1:7121906/93‑1 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  2:1035867/93‑1 (MQ=255)

                   ctatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCATCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACt                                                                      <  2:6068020/93‑1 (MQ=255)

                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     <  1:7154181/93‑1 (MQ=255)

                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     <  2:6418063/93‑1 (MQ=255)

                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTg                                                                    >  1:269307/1‑93 (MQ=255)

                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTg                                                                    >  1:787987/1‑93 (MQ=255)

                      tctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGc                                                                   >  1:1301113/1‑93 (MQ=255)

                       ctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGcc                                                                  >  2:5378690/1‑93 (MQ=255)

                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGccc                                                                 <  2:5122035/93‑1 (MQ=255)

                          tttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCgg                                                               >  2:3416231/1‑93 (MQ=255)

                          tttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCgg                                                               <  2:6804006/93‑1 (MQ=255)

                          tttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACACGACTTGCCCgg                                                               >  2:391081/1‑93 (MQ=255)

                            tCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGCCTTGCCCGGCa                                                             >  2:1503215/1‑93 (MQ=255)

                            tCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCa                                                             <  2:1663895/93‑1 (MQ=255)

                             cACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCAt                                                            <  2:5712106/93‑1 (MQ=255)

                              ccGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATc                                                           <  1:3127906/92‑1 (MQ=255)

                               cGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCa                                                          <  1:7549713/93‑1 (MQ=255)

                               cGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCa                                                          <  1:6856830/93‑1 (MQ=255)

                                 ggAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                                        >  2:243931/1‑93 (MQ=255)

                                  gAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgt                                                       <  2:1844138/93‑1 (MQ=255)

                                  gAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgt                                                       >  1:4310147/1‑93 (MQ=255)

                                    aTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                                                     <  2:5902032/93‑1 (MQ=255)

                                    aTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                                                     <  1:1051585/93‑1 (MQ=255)

                                       aaCTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTaa                                                  <  1:6645359/93‑1 (MQ=255)

                                       aaCTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTaa                                                  >  2:7346873/1‑93 (MQ=255)

                                       aaCTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTaa                                                  <  1:1013192/93‑1 (MQ=255)

                                          ttATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                                               <  1:1652907/93‑1 (MQ=255)

                                           tATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct                                              >  1:2737245/1‑93 (MQ=255)

                                            aTCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtctc                                             <  2:3150640/93‑1 (MQ=255)

                                             tCACCCTGCCAAAGTGCACCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGTTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTctct                                            <  2:4513399/93‑1 (MQ=255)

                                             tCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTctct                                            <  1:4983548/93‑1 (MQ=255)

                                              cACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCtt                                           <  1:6974002/93‑1 (MQ=255)

                                               aCCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCttt                                          <  1:4216483/93‑1 (MQ=255)

                                                cccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGAAGTAATTCTCTTTa                                         >  1:7704510/1‑93 (MQ=255)

                                                 ccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTac                                        <  2:6359290/93‑1 (MQ=255)

                                                 ccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTac                                        <  1:6264243/93‑1 (MQ=255)

                                                 ccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaa                                        >  1:1845496/1‑92 (MQ=255)

                                                  ctcccAAAGTGCACCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                                       <  2:3992983/90‑1 (MQ=255)

                                                  cTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                                       <  2:1735790/93‑1 (MQ=255)

                                                   tGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                      <  2:6266042/93‑1 (MQ=255)

                                                   tGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                      <  1:5225810/93‑1 (MQ=255)

                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGCGTCCGCGCGGGAAAAACGACGCAACTTGCCGAGCAACAGTAAGCAATTCTCTTTACACTct                                   >  2:1296626/1‑92 (MQ=255)

                                                       aaaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACcctcct                                  <  2:2770902/93‑1 (MQ=255)

                                                       aaaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCt                                  <  2:6118273/93‑1 (MQ=255)

                                                       aaaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCt                                  >  2:702494/1‑93 (MQ=255)

                                                        aaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                 <  2:2152259/93‑1 (MQ=255)

                                                         aGTGCACCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTa                                <  2:78588/93‑1 (MQ=255)

                                                          gTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTAc                               <  1:4842657/93‑1 (MQ=255)

                                                            gCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCg                             <  2:4256611/93‑1 (MQ=255)

                                                             cAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGt                            <  2:7707913/93‑1 (MQ=255)

                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           >  1:1985240/1‑93 (MQ=255)

                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           <  1:489436/93‑1 (MQ=255)

                                                               aCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTg                          >  2:7656718/1‑93 (MQ=255)

                                                               aCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCGTACCGTTg                          >  1:4441033/1‑93 (MQ=255)

                                                                cccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTgg                         <  2:7463979/93‑1 (MQ=255)

                                                                 ccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggc                        >  1:3360590/1‑93 (MQ=255)

                                                                  cGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcg                       >  1:6213250/1‑93 (MQ=255)

                                                                   gTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcgg                      >  1:9906/1‑93 (MQ=255)

                                                                    tCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcggc                     >  2:2725458/1‑93 (MQ=255)

                                                                     ccTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCt                    >  2:325291/1‑93 (MQ=255)

                                                                     ccTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCt                    >  1:5667532/1‑93 (MQ=255)

                                                                      cTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCtt                   >  2:3555539/1‑93 (MQ=255)

                                                                      cTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCtt                   >  2:3625807/1‑93 (MQ=255)

                                                                       tttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCttt                  >  1:2580083/1‑93 (MQ=255)

                                                                       tttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCttt                  >  1:1128779/1‑93 (MQ=255)

                                                                           ttATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAgg              <  2:5572812/93‑1 (MQ=255)

                                                                           ttATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAgg              >  1:5838386/1‑93 (MQ=255)

                                                                              tCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCaa           <  2:2824323/93‑1 (MQ=255)

                                                                                cGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAgg         >  1:83660/1‑93 (MQ=255)

                                                                                 gCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGt        <  1:293283/93‑1 (MQ=255)

                                                                                 gCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGt        >  2:1119693/1‑93 (MQ=255)

                                                                                  cAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTc       <  1:1559398/93‑1 (MQ=255)

                                                                                   aGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTcc      >  2:6749067/1‑93 (MQ=255)

                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  1:5547176/1‑93 (MQ=255)

                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  2:3770275/1‑93 (MQ=255)

                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  1:4073219/1‑93 (MQ=255)

                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGGAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  1:2805174/1‑93 (MQ=255)

                                                                                      gtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTa   <  2:1190664/93‑1 (MQ=255)

                                                                                      gtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTa   <  1:4781692/93‑1 (MQ=255)

                                                                                      gtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTa   <  1:3673346/93‑1 (MQ=255)

                                                                                       tCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAg  >  1:3584503/1‑93 (MQ=255)

                                                                                       tCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAg  >  2:1590072/1‑93 (MQ=255)

TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738

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: `‑`