BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,009,624	A→C	K287N (AAA→AAC)	ygeY →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,009,624	0	A	C	97.5%	540.1 / 2.1	158	K287N (AAA→AAC)	ygeY	putative peptidase
Reads supporting (aligned to +/- strand): ref base A (1/3); new base C (71/83); total (72/86)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 158 total aligned reads displayed.

CAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAA  >  NC_000913/3009527‑3009720

cAGACAGCTGCGCAGTCTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGc                                                                                                <  2:1589770/100‑1 (MQ=255)

  gACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTa                                                                                              >  1:4069133/1‑100 (MQ=255)

      gCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGCAATCCGCGCCCTGCCTGCAGTACAGAACGCTAAcgc                                                                                          >  1:919995/1‑100 (MQ=255)

       cTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAAcgcg                                                                                         <  2:670649/100‑1 (MQ=255)

         gcgcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgt                                                                                       <  1:1105814/100‑1 (MQ=255)

         gcgcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgt                                                                                       >  1:2733790/1‑100 (MQ=255)

         gcgcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgt                                                                                       >  1:141435/1‑100 (MQ=255)

           gcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGCAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgttg                                                                                     <  1:3877656/100‑1 (MQ=255)

           gcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgttg                                                                                     <  2:1322653/100‑1 (MQ=255)

           gcAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCGGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgttg                                                                                     <  2:924593/100‑1 (MQ=255)

             aGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGgttgtt                                                                                   <  1:4071460/100‑1 (MQ=255)

              gTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGttt                                                                                  >  1:2577192/1‑100 (MQ=255)

              gTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGttt                                                                                  >  1:3626390/1‑100 (MQ=255)

                ttCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCt                                                                                <  1:3951254/100‑1 (MQ=255)

                  cAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTAt                                                                              <  1:2077158/100‑1 (MQ=255)

                   aaTCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATg                                                                             >  1:2626000/1‑100 (MQ=255)

                    aTCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGt                                                                            <  2:139096/100‑1 (MQ=255)

                     tCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTa                                                                           >  2:2299879/1‑100 (MQ=255)

                         ccGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAAc                                                                       <  1:3997921/100‑1 (MQ=255)

                           gtcgtcTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACtc                                                                     >  2:1045491/1‑99 (MQ=255)

                           gtcgtcTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTa                                                                     >  2:2334854/1‑100 (MQ=255)

                           gtcgtcTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTa                                                                     >  2:3734155/1‑100 (MQ=255)

                            tcgtcTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTAc                                                                    <  2:1109429/100‑1 (MQ=255)

                             cgtcTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACg                                                                   <  1:1093004/100‑1 (MQ=255)

                                 tGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGAccg                                                               <  2:2291765/100‑1 (MQ=255)

                                   aCCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTAcgtccgtc                                                             <  2:3523563/100‑1 (MQ=255)

                                   aCCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGAccgtc                                                             <  1:2388901/100‑1 (MQ=255)

                                    ccTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGAccgtcc                                                            <  2:85777/100‑1 (MQ=255)

                                      tGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGAccgtccgt                                                          <  1:1492638/100‑1 (MQ=255)

                                       ggggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACcgtccgtc                                                         >  2:243298/1‑100 (MQ=255)

                                        gggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCgtccgtcc                                                        <  2:2218964/100‑1 (MQ=255)

                                         ggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCt                                                       >  2:2635682/1‑100 (MQ=255)

                                         ggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCt                                                       >  2:1682327/1‑100 (MQ=255)

                                         ggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCt                                                       >  2:1752454/1‑100 (MQ=255)

                                         ggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCt                                                       >  2:3249301/1‑100 (MQ=255)

                                          gCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTg                                                      <  1:2127496/100‑1 (MQ=255)

                                                ccTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTgg                                                <  2:3745173/100‑1 (MQ=255)

                                                 cTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACTACCGTCCGTCCTGGACTGGc                                               >  1:3301829/1‑100 (MQ=255)

                                                 cTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGc                                               >  1:117589/1‑100 (MQ=255)

                                                  taggAAGGCGCGCTGGACGAAAACCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTAGGTACAACTACGACCGTCCGTCCTGGACTGGcc                                              <  2:3191616/98‑1 (MQ=255)

                                                  tGGGAAGGCGCGCTGGACGAAATCCGCGCGCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGcc                                              <  2:1935912/100‑1 (MQ=255)

                                                  tGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGcc                                              <  1:1882559/100‑1 (MQ=255)

                                                   gggAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCt                                             >  2:299036/1‑100 (MQ=255)

                                                   gggAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCt                                             <  2:3228309/100‑1 (MQ=255)

                                                   gggAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCt                                             <  2:599545/100‑1 (MQ=255)

                                                    ggAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTg                                            >  1:2211939/1‑100 (MQ=255)

                                                    ggAAGGCGCGCTGAAAGAAAACAGCGACAAGAATGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTg                                            <  2:2809052/100‑1 (MQ=255)

                                                      aaGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGt                                          <  1:1012014/100‑1 (MQ=255)

                                                       aGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGtt                                         <  1:3775253/100‑1 (MQ=255)

                                                        ggCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGttt                                        >  1:1493467/1‑100 (MQ=255)

                                                          cgcgcTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTAc                                      <  2:1347535/100‑1 (MQ=255)

                                                            cgcTGGACGAAATCCGCGCACTGACAGCAATACAGAACGCTAACGCGGATGTTCCCATGTACAACTCGGACCGTCCGTCCTGGACTGGCCTGGTTTAccc                                    >  2:883179/1‑100 (MQ=255)

                                                               tGGACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAc                                 <  2:1828581/100‑1 (MQ=255)

                                                                ggACGAACTCCGCGCCCTGCCTGCAGTACAGAACGCTACCGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAcc                                <  2:2612256/100‑1 (MQ=255)

                                                                ggACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAcc                                >  2:233428/1‑100 (MQ=255)

                                                                ggACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAcc                                >  1:1307075/1‑100 (MQ=255)

                                                                ggACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAcc                                >  2:191660/1‑100 (MQ=255)

                                                                ggACGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAAcc                                <  1:2119280/100‑1 (MQ=255)

                                                                  aCGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGa                              >  1:1265316/1‑100 (MQ=255)

                                                                  aCGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGa                              >  1:2617552/1‑100 (MQ=255)

                                                                  aCGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGa                              >  1:541403/1‑100 (MQ=255)

                                                                   cGAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGaa                             >  2:2671292/1‑100 (MQ=255)

                                                                    gAAATCCGCGCCCTGCGTGCGGTACAAAACGCGAACGCGGTTGGTTCTATGTACAACTACGACCCTCCGTCCGGGAGTGGCCTGGTTTACCCAACCGAAt                            >  2:1441711/1‑100 (MQ=255)

                                                                    gAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAAt                            <  2:2692272/100‑1 (MQ=255)

                                                                    gAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAAt                            <  1:811179/100‑1 (MQ=255)

                                                                    gAAATCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAAt                            >  2:1198790/1‑100 (MQ=255)

                                                                     aaaTCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATg                           <  2:2976637/100‑1 (MQ=255)

                                                                      aaTCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGc                          <  1:3693560/100‑1 (MQ=255)

                                                                       aTCCGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCt                         >  2:3158326/1‑100 (MQ=255)

                                                                         ccGCGCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTAc                       >  1:507249/1‑100 (MQ=255)

                                                                          cgcgCCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACt                      <  2:766055/100‑1 (MQ=255)

                                                                          cgcgCCCTGCCCGCAGTACAGGACGGGAACGGGGGTGTTTGTAGGTACAATTACGACAGTCCTTCCTGGAATGGGCTGGTTTACCCAACCGAATGCTACt                      >  2:409504/1‑100 (MQ=255)

                                                                           gcgcCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACtt                     <  1:3146938/100‑1 (MQ=255)

                                                                             gcCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTcc                   <  1:655400/100‑1 (MQ=255)

                                                                             gcCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTcc                   <  1:1440094/100‑1 (MQ=255)

                                                                             gcCCTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGCCCTGGTTTACCCAACCGACTGCTACTccc                   >  2:1971984/1‑100 (MQ=255)

                                                                              cccTGCCTGCAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTccc                  >  2:3044516/1‑100 (MQ=255)

                                                                                 tgcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAc               <  1:103403/100‑1 (MQ=255)

                                                                                 tgcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAc               >  1:2906693/1‑100 (MQ=255)

                                                                                 tgcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAc               >  1:668992/1‑100 (MQ=255)

                                                                                  gcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAcc              <  1:301338/100‑1 (MQ=255)

                                                                                  gcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAcc              >  1:3634715/1‑100 (MQ=255)

                                                                                  gcctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGAcc              <  1:247269/100‑1 (MQ=255)

                                                                                   cctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCt             <  1:3361707/100‑1 (MQ=255)

                                                                                    ctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTg            <  1:154380/100‑1 (MQ=255)

                                                                                    ctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTg            <  2:323472/100‑1 (MQ=255)

                                                                                    ctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTg            <  2:1283172/100‑1 (MQ=255)

                                                                                    ctgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTg            <  1:649092/100‑1 (MQ=255)

                                                                                     tgcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTgg           >  2:2124339/1‑100 (MQ=255)

                                                                                      gcAGTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGa          >  1:3506093/1‑100 (MQ=255)

                                                                                         gTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       <  2:1126419/100‑1 (MQ=255)

                                                                                         gTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       <  1:1606088/100‑1 (MQ=255)

                                                                                         gTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       >  2:2739043/1‑100 (MQ=255)

                                                                                         gTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       >  1:2055195/1‑100 (MQ=255)

                                                                                         gTACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       >  1:3809465/1‑100 (MQ=255)

                                                                                         gTACAGAACACTAACGCGGTTGTTGCTATGTACAACTACGACCGTCCGTCCTAGACTGGCCTTGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAg       >  2:739393/1‑100 (MQ=255)

                                                                                          tACAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGt      >  1:2646470/1‑100 (MQ=255)

                                                                                            cAGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTgg    <  2:3548557/100‑1 (MQ=255)

                                                                                             aGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTCCTTCCCGACCTGGAAAGTGga   >  2:1169570/1‑100 (MQ=255)

                                                                                             aGAACGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGga   <  2:164403/100‑1 (MQ=255)

CAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAA  >  NC_000913/3009527‑3009720

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