BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,262,944	G→A	R242R (CGC→CGT)	tdcB ←	threonine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,262,944	0	G	A	99.1%	1184.7 / ‑3.6	339	R242R (CGC→CGT)	tdcB	threonine dehydratase
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (176/160); total (179/160)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 357 total aligned reads displayed.

GGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAACGTT  >  NC_000913/3262854‑3263036

 gCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCcggacgg                                                                                           <  2:7109487/93‑1 (MQ=255)

   aaTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGga                                                                                         <  2:6058767/93‑1 (MQ=255)

    aTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGgag                                                                                        >  2:4010880/1‑93 (MQ=255)

        tACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACAt                                                                                    >  1:5720068/1‑93 (MQ=255)

        tACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACAt                                                                                    >  1:62410/1‑93 (MQ=255)

         aCTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATc                                                                                   <  1:2965387/93‑1 (MQ=255)

             tttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACaa                                                                               <  2:3116005/93‑1 (MQ=255)

             tttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACaa                                                                               >  1:5243353/1‑93 (MQ=255)

             tttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACaa                                                                               >  1:5363121/1‑93 (MQ=255)

             tttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACaa                                                                               <  1:3149529/93‑1 (MQ=255)

               tCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTCAATTACCCGGACGGGAGACATCACAAcc                                                                             >  1:6675856/1‑93 (MQ=255)

               tCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAAcc                                                                             <  2:421655/93‑1 (MQ=255)

                cTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCa                                                                            >  1:7063851/1‑93 (MQ=255)

                 tGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCac                                                                           >  1:6540136/1‑92 (MQ=255)

                 tGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCAt                                                                           >  2:5774524/1‑93 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                          >  1:4335378/1‑93 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                          >  1:6178853/1‑93 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                          >  2:3899725/1‑93 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                          >  2:5928138/1‑93 (MQ=255)

                   aTTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATcc                                                                         <  2:7244514/93‑1 (MQ=255)

                   aTTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATcc                                                                         >  1:900035/1‑93 (MQ=255)

                    tttCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCg                                                                        >  2:622720/1‑93 (MQ=255)

                     ttCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTCCCCGGACGGGAGACATCACAACCATCCGc                                                                       <  2:3407174/93‑1 (MQ=255)

                     ttCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGc                                                                       >  2:295060/1‑93 (MQ=255)

                       cGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCa                                                                     <  1:7217075/93‑1 (MQ=255)

                        gTCTTCTCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  2:1682841/93‑1 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  1:4607416/93‑1 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  1:2174299/93‑1 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  2:6877159/93‑1 (MQ=255)

                        gTCTGCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  1:7362361/93‑1 (MQ=255)

                         tCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                    <  1:7541696/92‑1 (MQ=255)

                          cTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAggg                                                                  <  2:6394134/93‑1 (MQ=255)

                          cTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAggg                                                                  <  1:1583541/93‑1 (MQ=255)

                          cTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAggg                                                                  >  1:1155258/1‑93 (MQ=255)

                              gCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGcc                                                              <  1:2435061/93‑1 (MQ=255)

                               cTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCg                                                             <  2:7695892/93‑1 (MQ=255)

                                tGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCgg                                                            >  2:1445835/1‑93 (MQ=255)

                                tGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCgg                                                            >  1:6888376/1‑93 (MQ=255)

                                tGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCgg                                                            <  1:1771208/93‑1 (MQ=255)

                                  aCCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTa                                                          <  2:5256016/93‑1 (MQ=255)

                                   ccAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGCGACATCACAACCATCCGCCAGGGTGCCGGTAg                                                         >  1:3740073/1‑93 (MQ=255)

                                   ccAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAg                                                         <  1:4888614/93‑1 (MQ=255)

                                   ccAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAg                                                         <  1:2563699/93‑1 (MQ=255)

                                    cAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGt                                                        >  2:6514692/1‑93 (MQ=255)

                                     aGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGtt                                                       >  2:6594865/1‑93 (MQ=255)

                                     aGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGtt                                                       <  1:2271623/93‑1 (MQ=255)

                                      gCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTc                                                      <  2:6728521/93‑1 (MQ=255)

                                      gCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTc                                                      >  2:2512443/1‑93 (MQ=255)

                                       cACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCg                                                     <  2:5173731/93‑1 (MQ=255)

                                        aCGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCgg                                                    <  2:1356902/93‑1 (MQ=255)

                                        aCGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCgg                                                    >  2:4760261/1‑93 (MQ=255)

                                         cGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGt                                                   <  2:7516715/93‑1 (MQ=255)

                                            tGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCg                                                >  2:2706743/1‑93 (MQ=255)

                                              tCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTg                                              <  2:2661831/93‑1 (MQ=255)

                                              tCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTg                                              >  1:1233555/1‑93 (MQ=255)

                                               cATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTgg                                             <  1:3917467/93‑1 (MQ=255)

                                                aTCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGt                                            <  2:1918383/93‑1 (MQ=255)

                                                 tCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGtt                                           >  2:5293984/1‑93 (MQ=255)

                                                 tCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGtt                                           >  1:152331/1‑93 (MQ=255)

                                                  cGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTa                                          >  1:5605914/1‑93 (MQ=255)

                                                  cGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTa                                          <  2:7282303/93‑1 (MQ=255)

                                                   gACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGAAGGGAGACATCAAAACCATCCGCCAGGGTGCCGTAAGTTCGGTGCGTGGTTAt                                         >  2:6768144/1‑93 (MQ=255)

                                                       aaTTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTtct                                     <  1:2705588/93‑1 (MQ=255)

                                                        aTTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTtctc                                    <  2:3815411/93‑1 (MQ=255)

                                                         ttCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTcc                                   <  2:234160/93‑1 (MQ=255)

                                                           cACGAACTATTTCGTAAGTTAAATTACCCGGACGGGATACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCgg                                 <  2:4713469/93‑1 (MQ=255)

                                                           cACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCgg                                 <  2:6311535/93‑1 (MQ=255)

                                                           cACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCgg                                 <  1:2123345/93‑1 (MQ=255)

                                                            acgaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCTGGTAGTTCGGTGCGTGGTTATTTCTCCGGa                                <  2:1021711/93‑1 (MQ=255)

                                                            acgaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGa                                <  1:1055040/93‑1 (MQ=255)

                                                                 cgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGa                           >  2:5687930/1‑93 (MQ=255)

                                                                   aTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGaaa                         >  1:260242/1‑93 (MQ=255)

                                                                   aTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGaaa                         <  1:4207411/93‑1 (MQ=255)

                                                                     ttCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaaga                       <  2:2124268/93‑1 (MQ=255)

                                                                      tCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaa                      <  1:362596/93‑1 (MQ=255)

                                                                      tCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaa                      <  2:3687762/93‑1 (MQ=255)

                                                                       cGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaag                     >  2:791250/1‑93 (MQ=255)

                                                                        gTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgc                    >  2:3245193/1‑93 (MQ=255)

                                                                        gTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgc                    <  1:1307943/93‑1 (MQ=255)

                                                                        gTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGCGCGTGGTTATTTCTCCGGAGTGGAAAGAAgc                    <  1:6639183/93‑1 (MQ=255)

                                                                         tAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaagac                   >  2:5506948/1‑91 (MQ=255)

                                                                         tAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgcc                   >  1:5918585/1‑93 (MQ=255)

                                                                          aaGTTAAATTACCCGGACGGGAGACATCACACCCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgccg                  <  2:7138282/93‑1 (MQ=255)

                                                                           aGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgccgc                 >  2:2615376/1‑93 (MQ=255)

                                                                             ttAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCa               <  1:643410/93‑1 (MQ=255)

                                                                             ttAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCa               <  1:4526437/93‑1 (MQ=255)

                                                                             ttAAATTACCCGGACGGGAGACATCACAACCATCCGCCACGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCa               >  1:3303092/1‑93 (MQ=255)

                                                                              tAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCAt              >  1:1868640/1‑93 (MQ=255)

                                                                              tAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCAt              <  1:3405653/93‑1 (MQ=255)

                                                                                aaTTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGc            <  2:446765/93‑1 (MQ=255)

                                                                                 aTTACCCGGACGGGAGACATCACAACCATCCGCCAGGGGGCCGGGAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATgac           >  2:2772018/1‑91 (MQ=255)

                                                                                   tACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGt         >  2:3632662/1‑93 (MQ=255)

                                                                                   tACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGt         <  1:7680084/93‑1 (MQ=255)

                                                                                    aCCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGAGGAAAGAAGCCGCCATGCCGTg        >  1:4766516/1‑93 (MQ=255)

                                                                                      ccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGac      >  1:5807069/1‑92 (MQ=255)

                                                                                      ccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGaa      <  1:7423607/93‑1 (MQ=255)

                                                                                       cggacgggAGACATCACACCCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAAc     <  2:1591301/93‑1 (MQ=255)

                                                                                       cggacgggAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAAc     <  1:825166/93‑1 (MQ=255)

                                                                                          gcgggAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAACGtt  <  2:5409249/93‑1 (MQ=255)

                                                                                          gcgggAGACATCACAACCAGCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAACGtt  <  1:3202479/93‑1 (MQ=255)

GGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAACGTT  >  NC_000913/3262854‑3263036

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