BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,008,642	C→T	A198A (GCC→GCT)	b2873 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,008,642	0	C	T	99.4%	637.2 / ‑5.5	178	A198A (GCC→GCT)	b2873	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (72/105); total (73/105)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 181 total aligned reads displayed.

ATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCAC  >  NC_000913/3008552‑3008730

aTTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATcgctgc                                                                                        >  2:3216497/1‑93 (MQ=255)

 ttACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCa                                                                                       <  2:1471690/93‑1 (MQ=255)

 ttACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCa                                                                                       <  2:144077/93‑1 (MQ=255)

   aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAgg                                                                                     <  2:1570628/93‑1 (MQ=255)

   aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAgg                                                                                     <  2:376294/93‑1 (MQ=255)

     gccgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGt                                                                                   <  2:2875255/93‑1 (MQ=255)

     gccgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGt                                                                                   <  2:1404148/93‑1 (MQ=255)

       cgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTa                                                                                 >  2:1658003/1‑93 (MQ=255)

       cgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTa                                                                                 <  1:38522/93‑1 (MQ=255)

       cgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGGTa                                                                                 >  1:3332302/1‑93 (MQ=255)

        gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTaa                                                                                <  1:4059151/93‑1 (MQ=255)

        gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTaa                                                                                >  2:658482/1‑93 (MQ=255)

         cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAc                                                                               <  1:1220765/93‑1 (MQ=255)

         cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAc                                                                               <  1:647623/93‑1 (MQ=255)

          ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAcc                                                                              <  1:2669824/93‑1 (MQ=255)

            gCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgc                                                                            <  1:4601725/93‑1 (MQ=255)

              aTGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgc                                                                          >  1:606790/1‑93 (MQ=255)

               tGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgcc                                                                         <  1:1465524/93‑1 (MQ=255)

                       gcgcgcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATc                                                                 <  2:2771654/93‑1 (MQ=255)

                         gcgcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCAt                                                               <  2:2125968/93‑1 (MQ=255)

                         gcgcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCAt                                                               <  1:44916/93‑1 (MQ=255)

                           gcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGc                                                             >  1:427340/1‑93 (MQ=255)

                           gcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGc                                                             >  1:3756406/1‑93 (MQ=255)

                           gcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGc                                                             >  1:2676826/1‑93 (MQ=255)

                           gcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGc                                                             >  2:2024247/1‑93 (MQ=255)

                            cTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                            <  1:2338260/93‑1 (MQ=255)

                            cTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                            <  1:1080353/93‑1 (MQ=255)

                              gACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCtt                                                          <  2:2460839/93‑1 (MQ=255)

                              gACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCtt                                                          <  2:2978905/93‑1 (MQ=255)

                               accaccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTg                                                         <  1:2742715/93‑1 (MQ=255)

                                ccaccGTGCGCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGa                                                        <  1:4158334/93‑1 (MQ=255)

                                 caccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAg                                                       <  2:4024563/93‑1 (MQ=255)

                                  accGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGt                                                      >  2:484571/1‑93 (MQ=255)

                                  accGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGt                                                      >  2:3529952/1‑93 (MQ=255)

                                  accGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGt                                                      >  1:793915/1‑93 (MQ=255)

                                  accGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGt                                                      <  1:2248803/93‑1 (MQ=255)

                                   ccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATTGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTc                                                     <  2:199772/93‑1 (MQ=255)

                                   ccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTc                                                     <  1:3410232/93‑1 (MQ=255)

                                   ccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTc                                                     >  1:4324294/1‑93 (MQ=255)

                                      tGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGcc                                                  >  1:243103/1‑93 (MQ=255)

                                      tGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGcc                                                  >  1:78944/1‑93 (MQ=255)

                                      tGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGcc                                                  >  1:201318/1‑93 (MQ=255)

                                        cACCCGGAAGATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCct                                                <  1:3022954/93‑1 (MQ=255)

                                        cACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCct                                                >  2:4115087/1‑93 (MQ=255)

                                         ccccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctc                                               <  2:4190876/92‑1 (MQ=255)

                                         aCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctc                                               >  2:2321313/1‑93 (MQ=255)

                                         aCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctc                                               <  1:1090869/93‑1 (MQ=255)

                                         aCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGACTTGAGTCGCCctc                                               >  2:1852380/1‑93 (MQ=255)

                                          cccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctct                                              <  1:3648572/93‑1 (MQ=255)

                                           ccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTg                                             <  1:414814/93‑1 (MQ=255)

                                              gAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGaa                                          <  2:3380151/93‑1 (MQ=255)

                                               aaaaTGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAAt                                         <  2:4160680/93‑1 (MQ=255)

                                                aaaTGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATg                                        >  1:2974144/1‑93 (MQ=255)

                                                 aaTGATGCGGCTACCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATgg                                       >  2:3150011/1‑92 (MQ=255)

                                                  atgatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCg                                      >  1:1335101/1‑93 (MQ=255)

                                                   tgatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGa                                     >  1:1558071/1‑93 (MQ=255)

                                                    gatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGaa                                    <  1:3788957/93‑1 (MQ=255)

                                                    gatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGaa                                    >  2:622569/1‑93 (MQ=255)

                                                    gatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGaa                                    <  2:3865656/93‑1 (MQ=255)

                                                     atgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAg                                   <  2:2226582/93‑1 (MQ=255)

                                                      tgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                  <  2:4545096/93‑1 (MQ=255)

                                                      tgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                  >  2:2658096/1‑93 (MQ=255)

                                                      tgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                  <  1:2484408/93‑1 (MQ=255)

                                                      tgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                  <  1:4514853/93‑1 (MQ=255)

                                                      tgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                  <  2:149324/93‑1 (MQ=255)

                                                       gCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCg                                 >  2:3506378/1‑93 (MQ=255)

                                                          gCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGaa                              <  1:1983142/93‑1 (MQ=255)

                                                           cTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAg                             >  1:2916731/1‑93 (MQ=255)

                                                           cTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAg                             >  1:2373454/1‑93 (MQ=255)

                                                            tATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGc                            <  1:301254/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        >  2:671104/1‑93 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTAAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:2594667/93‑1 (MQ=255)

                                                                 ccAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCg                       <  1:877644/93‑1 (MQ=255)

                                                                   aGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGcc                     <  1:3461725/93‑1 (MQ=255)

                                                                   aGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGcc                     >  2:1140079/1‑93 (MQ=255)

                                                                    gCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                    <  2:4630969/93‑1 (MQ=255)

                                                                       aGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 >  1:3208193/1‑93 (MQ=255)

                                                                       aGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 <  2:4388323/93‑1 (MQ=255)

                                                                       aGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 >  2:1034214/1‑93 (MQ=255)

                                                                          gggcggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGa              <  2:928973/93‑1 (MQ=255)

                                                                          gggcggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGa              <  1:2255045/93‑1 (MQ=255)

                                                                           ggcggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGAt             >  1:117965/1‑93 (MQ=255)

                                                                             cggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTa           <  1:3540168/93‑1 (MQ=255)

                                                                              ggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTaa          <  2:1156339/93‑1 (MQ=255)

                                                                              ggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTaa          <  2:1748033/93‑1 (MQ=255)

                                                                               gAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATTATTAAc         >  2:3067822/1‑93 (MQ=255)

                                                                               gAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAc         <  1:2497540/93‑1 (MQ=255)

                                                                               gAATTTATCGCTGCAGGGGTAACCGCGCCGCGCTATCATGCCATGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAc         >  1:3297197/1‑93 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        <  1:3955462/93‑1 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        >  1:3396193/1‑93 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        >  1:124960/1‑93 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        >  2:2071534/1‑93 (MQ=255)

                                                                                 aTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCt       <  1:1950848/93‑1 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTg      >  2:3207137/1‑93 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTg      <  2:1103894/93‑1 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTg      >  1:1661448/1‑93 (MQ=255)

                                                                                   ttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTgg     >  1:2145606/1‑93 (MQ=255)

                                                                                   ttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTgg     <  1:3094896/93‑1 (MQ=255)

                                                                                    tATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGc    >  1:1309632/1‑93 (MQ=255)

                                                                                    tATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGc    >  1:3322450/1‑93 (MQ=255)

ATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCAC  >  NC_000913/3008552‑3008730

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