BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,787,955	A→G	W340R (TGG→CGG)	tdh ←	threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,787,955	0	A	G	85.9%	810.1 / 124.9	340	W340R (TGG→CGG)	tdh	threonine dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (23/25); new base G (140/152); total (163/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 343 total aligned reads displayed.

AACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCG  >  NC_000913/3787863‑3788047

aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCa                                                                                              >  2:4855894/1‑93 (MQ=255)

  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGc                                                                                            <  2:4511599/93‑1 (MQ=255)

  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGc                                                                                            >  2:3316545/1‑93 (MQ=255)

  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGc                                                                                            <  1:6076925/93‑1 (MQ=255)

  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGc                                                                                            <  1:2034011/93‑1 (MQ=255)

     tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCa                                                                                         >  1:3464913/1‑93 (MQ=255)

         ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                                     >  2:4355995/1‑93 (MQ=255)

           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATaa                                                                                   <  1:7512091/93‑1 (MQ=255)

           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATaa                                                                                   >  1:5918430/1‑93 (MQ=255)

            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAAc                                                                                  >  1:4191173/1‑93 (MQ=255)

            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                  >  2:6536500/1‑93 (MQ=255)

            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                  <  1:5206383/93‑1 (MQ=255)

                gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTc                                                                              >  1:2652980/1‑93 (MQ=255)

                gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTc                                                                              >  2:2210780/1‑93 (MQ=255)

                gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTc                                                                              >  1:6672455/1‑93 (MQ=255)

                 aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTcc                                                                             <  2:4883999/93‑1 (MQ=255)

                 aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTcc                                                                             >  2:3077870/1‑93 (MQ=255)

                  cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTccc                                                                            <  2:1337892/93‑1 (MQ=255)

                  cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTccc                                                                            <  1:1816926/93‑1 (MQ=255)

                  cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTccc                                                                            <  1:7351217/93‑1 (MQ=255)

                     gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGa                                                                         >  1:6133590/1‑93 (MQ=255)

                       ttttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACt                                                                       >  2:2598033/1‑93 (MQ=255)

                       ttttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACt                                                                       >  2:6620947/1‑93 (MQ=255)

                       ttttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACt                                                                       >  2:1570264/1‑93 (MQ=255)

                       ttttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                                       <  1:6802834/93‑1 (MQ=255)

                         ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTgg                                                                     <  1:3629853/93‑1 (MQ=255)

                           aCATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGcc                                                                   <  1:4402892/93‑1 (MQ=255)

                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGccc                                                                  <  1:7437465/93‑1 (MQ=255)

                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGccc                                                                  <  2:5758383/93‑1 (MQ=255)

                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGccc                                                                  >  2:6265878/1‑93 (MQ=255)

                             atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCg                                                                 >  1:701442/1‑93 (MQ=255)

                             atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                                                 >  2:754305/1‑93 (MQ=255)

                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGa                                                                >  1:3940920/1‑93 (MQ=255)

                                tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                                              >  2:4090994/1‑93 (MQ=255)

                                 ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                             >  2:2855958/1‑93 (MQ=255)

                                 ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                             >  1:2141425/1‑93 (MQ=255)

                                  tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                                            <  2:496962/93‑1 (MQ=255)

                                     cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                         <  2:5281073/93‑1 (MQ=255)

                                       gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                                       <  2:5107606/93‑1 (MQ=255)

                                          tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                                    >  2:5434465/1‑93 (MQ=255)

                                            tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                  <  2:4752805/93‑1 (MQ=255)

                                            tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                  <  1:1071766/93‑1 (MQ=255)

                                            tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                  >  2:7414530/1‑93 (MQ=255)

                                             ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                 <  2:120230/93‑1 (MQ=255)

                                             ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                 <  2:768399/93‑1 (MQ=255)

                                              cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                                                >  2:6914419/1‑93 (MQ=255)

                                               tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                                               <  2:2122305/93‑1 (MQ=255)

                                                 aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                                             >  2:1504703/1‑93 (MQ=255)

                                                 aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                                             >  1:2207209/1‑93 (MQ=255)

                                                 aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                                             >  1:2897531/1‑93 (MQ=255)

                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                                            <  2:172201/93‑1 (MQ=255)

                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                                            <  1:808737/93‑1 (MQ=255)

                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                                            >  1:1723429/1‑93 (MQ=255)

                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                                           >  1:1637258/1‑93 (MQ=255)

                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                                           <  1:2994500/93‑1 (MQ=255)

                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                          >  2:498034/1‑93 (MQ=255)

                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                          >  2:3578964/1‑93 (MQ=255)

                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                          >  2:2933526/1‑93 (MQ=255)

                                                     aaaaGGGCTGGAATACCAGCCCTTGTTCGTTTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                                         <  1:3893439/93‑1 (MQ=255)

                                                     aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                                         <  1:5360065/93‑1 (MQ=255)

                                                      aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                                        >  2:3172497/1‑93 (MQ=255)

                                                      aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                                        >  2:4625185/1‑93 (MQ=255)

                                                          ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                    >  1:3493842/1‑93 (MQ=255)

                                                           gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                   >  1:3827956/1‑93 (MQ=255)

                                                            cTGGAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                                  <  2:2728210/93‑1 (MQ=255)

                                                              ggAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                >  2:5375041/1‑93 (MQ=255)

                                                              ggAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                <  2:5883995/93‑1 (MQ=255)

                                                               gAATACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                               <  2:6338160/93‑1 (MQ=255)

                                                                 aTACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                             <  1:5667135/93‑1 (MQ=255)

                                                                 aTACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                             <  2:5599110/93‑1 (MQ=255)

                                                                  tACCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt                            <  1:3103937/93‑1 (MQ=255)

                                                                   aCCAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa                           >  2:4145888/1‑93 (MQ=255)

                                                                     cAGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                         <  2:6141384/93‑1 (MQ=255)

                                                                      aGCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag                        >  1:40222/1‑93 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       >  1:6364569/1‑93 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       <  1:1279144/93‑1 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       <  2:2343980/93‑1 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       >  1:3521797/1‑93 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       <  1:1911522/93‑1 (MQ=255)

                                                                       gCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                       >  1:4523207/1‑93 (MQ=255)

                                                                        cccTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaa                      <  2:3854716/93‑1 (MQ=255)

                                                                         ccTTGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa                     >  2:7229094/1‑93 (MQ=255)

                                                                           ttGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                   >  2:5963955/1‑93 (MQ=255)

                                                                           ttGTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                   >  2:6081519/1‑93 (MQ=255)

                                                                             gTTCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGAt                 <  1:6549139/93‑1 (MQ=255)

                                                                              ttCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATg                <  2:3535613/93‑1 (MQ=255)

                                                                              ttCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATg                <  2:5819099/93‑1 (MQ=255)

                                                                               tCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg               >  1:6116559/1‑93 (MQ=255)

                                                                               tCGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg               >  2:1909769/1‑93 (MQ=255)

                                                                               tCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg               <  1:4985363/93‑1 (MQ=255)

                                                                                cGTGTTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATggg              >  1:1134379/1‑93 (MQ=255)

                                                                                 gtgtTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGt             >  1:6480661/1‑93 (MQ=255)

                                                                                 gtgtTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGt             >  2:1460483/1‑93 (MQ=255)

                                                                                  tgtTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTa            >  2:7117544/1‑93 (MQ=255)

                                                                                   gtTAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTaa           <  1:7136016/93‑1 (MQ=255)

                                                                                    ttAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAAt          <  2:4033549/93‑1 (MQ=255)

                                                                                     tAATCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATg         >  1:3312667/1‑93 (MQ=255)

                                                                                       aTCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAt       >  2:3774327/1‑93 (MQ=255)

                                                                                        tCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATc      >  1:5079629/1‑93 (MQ=255)

                                                                                        tCCCGGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATc      >  1:952032/1‑93 (MQ=255)

AACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCG  >  NC_000913/3787863‑3788047

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