BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	161,581	C→T	E151K (GAG→AAG)	yadP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	161,581	0	C	T	99.6%	889.9 / NA	239	E151K (GAG→AAG)	yadP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (106/132); minor base A (1/0); total (107/132)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

ACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACG  >  NC_000913/161491‑161670

aCAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCg                                                                                         >  2:4512454/1‑93 (MQ=255)

 cAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGt                                                                                        >  2:605943/1‑93 (MQ=255)

 cAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGt                                                                                        >  1:4917995/1‑93 (MQ=255)

 cAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGt                                                                                        >  2:5564564/1‑93 (MQ=255)

 cAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGt                                                                                        >  2:3602535/1‑93 (MQ=255)

   aTCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCa                                                                                      <  2:69093/93‑1 (MQ=255)

   aTCCTTGTTATTGCGTTAGCGCCCAGCGGTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCa                                                                                      >  1:3886321/1‑93 (MQ=255)

    tCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAc                                                                                     <  1:3415652/93‑1 (MQ=255)

    tCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAc                                                                                     <  1:2152025/93‑1 (MQ=255)

     ccTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAcc                                                                                    <  2:1286401/93‑1 (MQ=255)

       ttGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGc                                                                                  >  1:1806010/1‑93 (MQ=255)

         gTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCAt                                                                                >  2:3997133/1‑93 (MQ=255)

         gTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGCGAACTTCGTCACCGCAt                                                                                >  2:4647162/1‑93 (MQ=255)

          ttATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATa                                                                               >  1:1752753/1‑93 (MQ=255)

          ttATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATa                                                                               >  2:5162242/1‑93 (MQ=255)

           tATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATAc                                                                              <  1:1600039/93‑1 (MQ=255)

           tATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATAc                                                                              <  1:1002302/93‑1 (MQ=255)

             ttGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGa                                                                            <  1:4368707/93‑1 (MQ=255)

                 gTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAg                                                                        <  1:2391881/93‑1 (MQ=255)

                 gTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAg                                                                        >  1:5750749/1‑93 (MQ=255)

                  ttAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGt                                                                       >  2:5755673/1‑93 (MQ=255)

                    aGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTa                                                                     <  2:149043/93‑1 (MQ=255)

                       gcCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTaaaa                                                                  <  1:1104359/93‑1 (MQ=255)

                        cccAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAAc                                                                 <  2:829435/93‑1 (MQ=255)

                         ccAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAAcc                                                                >  2:4115204/1‑93 (MQ=255)

                         ccAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAAcc                                                                <  1:861915/93‑1 (MQ=255)

                          cAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCt                                                               <  2:4382051/93‑1 (MQ=255)

                          cAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCt                                                               <  1:1382987/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              <  1:1318994/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              >  2:1028749/1‑93 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              <  2:666943/93‑1 (MQ=255)

                            gCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTgg                                                             <  1:4591427/93‑1 (MQ=255)

                            gCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTgg                                                             <  2:3210070/93‑1 (MQ=255)

                             cGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggc                                                            >  1:912469/1‑93 (MQ=255)

                               ttttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcgg                                                          >  2:1657341/1‑93 (MQ=255)

                               ttttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcgg                                                          >  1:5846537/1‑93 (MQ=255)

                               ttttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcgg                                                          >  2:3921857/1‑93 (MQ=255)

                                tttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcggc                                                         <  2:5144856/93‑1 (MQ=255)

                                tttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcggc                                                         >  2:3957954/1‑93 (MQ=255)

                                 ttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCCCATACGACCAGTTAAAACCTAgcggcg                                                        >  2:2334268/1‑93 (MQ=255)

                                  tAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGcggcgg                                                       >  2:4894158/1‑93 (MQ=255)

                                    gcggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGa                                                     <  2:1402422/93‑1 (MQ=255)

                                    gcggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGa                                                     <  2:3824206/93‑1 (MQ=255)

                                    gcggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGa                                                     <  1:5039320/93‑1 (MQ=255)

                                     cggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGAt                                                    >  2:5080585/1‑93 (MQ=255)

                                      ggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGAtt                                                   <  1:4644954/93‑1 (MQ=255)

                                          tgtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCa                                               >  2:6022336/1‑93 (MQ=255)

                                          tgtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCa                                               >  2:4670171/1‑93 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAc                                              <  1:5856856/93‑1 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAc                                              <  2:2316563/93‑1 (MQ=255)

                                            tAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAcc                                             >  1:5979444/1‑93 (MQ=255)

                                                gtgtGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCt                                         <  1:5536042/93‑1 (MQ=255)

                                                 tgtgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTc                                        <  2:3145339/93‑1 (MQ=255)

                                                   tgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGc                                      <  1:6222681/93‑1 (MQ=255)

                                                   tgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGc                                      <  1:617081/93‑1 (MQ=255)

                                                   tgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGc                                      <  2:2989217/93‑1 (MQ=255)

                                                     cGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCGGGCGGCGGGATTGTCACCGCCTCGCTg                                    <  1:3527364/93‑1 (MQ=255)

                                                      gTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTgg                                   <  2:5404161/93‑1 (MQ=255)

                                                      gTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTgg                                   >  2:777642/1‑93 (MQ=255)

                                                      gTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTgg                                   <  1:3205905/93‑1 (MQ=255)

                                                       tCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGc                                  >  1:6151769/1‑93 (MQ=255)

                                                       tCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGc                                  >  1:6046067/1‑93 (MQ=255)

                                                       tCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGc                                  >  1:5279808/1‑93 (MQ=255)

                                                         cACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGt                                <  2:4353908/93‑1 (MQ=255)

                                                          aCGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTc                               <  1:1512415/93‑1 (MQ=255)

                                                           cGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATGGTCACGGCCTCGCTGGCGTcg                              >  2:5365603/1‑93 (MQ=255)

                                                             ggCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTcgcg                            >  1:550134/1‑93 (MQ=255)

                                                             ggCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTcgcg                            >  2:4199484/1‑93 (MQ=255)

                                                              gCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCgcgc                           >  2:4899731/1‑93 (MQ=255)

                                                               cAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCa                          >  1:2767644/1‑93 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        <  1:2722520/93‑1 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        <  2:1615704/93‑1 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        <  1:1897296/93‑1 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        <  2:383752/93‑1 (MQ=255)

                                                                  aCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATa                       >  2:4469076/1‑93 (MQ=255)

                                                                    gaggagGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAg                     <  1:3051815/93‑1 (MQ=255)

                                                                    gaggagGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAg                     <  1:985031/93‑1 (MQ=255)

                                                                     aggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAgg                    >  1:1466034/1‑93 (MQ=255)

                                                                     aggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAgg                    <  1:6277138/93‑1 (MQ=255)

                                                                      ggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGt                   <  1:4118594/93‑1 (MQ=255)

                                                                      ggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGt                   <  1:1217581/93‑1 (MQ=255)

                                                                         ggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAat                >  1:11272/1‑93 (MQ=255)

                                                                          gCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAata               <  1:4055207/93‑1 (MQ=255)

                                                                          gCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAata               <  2:272930/93‑1 (MQ=255)

                                                                           cGTAAAGGGTGATCTTCGTCACCGCATACGACCAGTTAAACCCCGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAatat              <  1:5417389/93‑1 (MQ=255)

                                                                           cGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAatat              <  2:2544199/93‑1 (MQ=255)

                                                                           cGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAatat              >  1:2670834/1‑93 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGCAATAtgt            >  1:5236922/1‑93 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  2:3481399/93‑1 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  2:1877201/93‑1 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  2:3224461/93‑1 (MQ=255)

                                                                               aaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTgg          <  2:1959120/93‑1 (MQ=255)

                                                                                aGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGa         <  2:2089751/93‑1 (MQ=255)

                                                                                 gggTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGAt        >  1:1455070/1‑93 (MQ=255)

                                                                                  ggTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATg       <  1:4545436/93‑1 (MQ=255)

                                                                                  ggTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATg       >  1:2620047/1‑93 (MQ=255)

                                                                                    tGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaa     >  2:2194072/1‑93 (MQ=255)

                                                                                    tGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGACACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaa     <  1:1198335/93‑1 (MQ=255)

                                                                                     gAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaaa    <  1:4699958/93‑1 (MQ=255)

                                                                                       aCTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACg  >  2:3062299/1‑93 (MQ=255)

ACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACG  >  NC_000913/161491‑161670

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