BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,837,903	A→G	K138R (AAA→AGA)	b1757 →	putative thiosulfate sulfur transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,837,903	0	A	G	100.0%	790.6 / NA	241	K138R (AAA→AGA)	b1757	putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (121/120); total (121/120)

Only 100 of 246 total aligned reads displayed.

TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992

tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         >  2:5326462/1‑93 (MQ=255)

 cGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACt                                                                                        <  1:9096992/93‑1 (MQ=255)

 cGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACt                                                                                        <  1:7273845/93‑1 (MQ=255)

 cGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACt                                                                                        <  1:6121929/93‑1 (MQ=255)

 cGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCCGCAAAGACt                                                                                        <  1:6664901/93‑1 (MQ=255)

    cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGcc                                                                                     <  2:9786200/93‑1 (MQ=255)

     gCCGTCAAAACGCGACTGCAAAAAGCTGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCg                                                                                    <  2:2484339/93‑1 (MQ=255)

      ccGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGc                                                                                   <  2:7250511/93‑1 (MQ=255)

         tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  2:5377276/1‑93 (MQ=255)

               cgcgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagc                                                                          <  1:4995975/93‑1 (MQ=255)

                gcgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTTCAAAGACTGCCGCATTTTGagca                                                                         >  2:5600807/1‑93 (MQ=255)

                gcgACTGCAAAAAGCAGGGTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACATTCCCGTCTGCCAAGACTGCCGCATTCTGagcg                                                                         >  2:4288209/1‑92 (MQ=255)

                  gACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcagc                                                                       >  2:3168764/1‑93 (MQ=255)

                  gACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcagc                                                                       <  1:3821925/93‑1 (MQ=255)

                  gACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcagc                                                                       >  2:9433401/1‑93 (MQ=255)

                   aCTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCt                                                                      <  2:1179482/93‑1 (MQ=255)

                   aCTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCt                                                                      >  1:1872741/1‑93 (MQ=255)

                    cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     <  2:3556131/93‑1 (MQ=255)

                    cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     >  1:2409434/1‑93 (MQ=255)

                    cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     >  2:8363016/1‑93 (MQ=255)

                       cAAAAAGCAGGTTTAACGCATTTCGCCATCCTGAGTGGCGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  2:8255899/93‑1 (MQ=255)

                        aaaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGttt                                                                 >  2:6686707/1‑93 (MQ=255)

                         aaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTa                                                                <  1:9697303/93‑1 (MQ=255)

                         aaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTa                                                                >  2:3991539/1‑93 (MQ=255)

                          aaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTAt                                                               >  1:5703404/1‑93 (MQ=255)

                          aaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTAt                                                               >  2:5445619/1‑93 (MQ=255)

                              cAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGc                                                           >  2:6395643/1‑93 (MQ=255)

                              cAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGc                                                           <  1:3714577/93‑1 (MQ=255)

                                ggTTTCCCGCTTCTCTCCCTCCTGAGTGCCGCGCTCCGCGGACCTTCCCGTCTGCCACGCCTGCCGCGTTTTGCGCAGCTGGTTTATCCGccc                                                         <  1:7702125/93‑3 (MQ=255)

                                ggTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCaa                                                         <  2:6458415/93‑1 (MQ=255)

                                 gTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAt                                                        <  1:7132821/93‑1 (MQ=255)

                                 gTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAt                                                        >  1:2889303/1‑93 (MQ=255)

                                  tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       <  1:7471302/93‑1 (MQ=255)

                                  tttAACGCATATCTCCATCCTGAGTGACGCCCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       <  1:8524232/93‑1 (MQ=255)

                                   ttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATgg                                                      <  2:9327301/93‑1 (MQ=255)

                                    tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     >  2:3431430/1‑93 (MQ=255)

                                      aCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTg                                                   >  1:2470137/1‑93 (MQ=255)

                                       cGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGc                                                  <  2:1644258/93‑1 (MQ=255)

                                         cATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAc                                                >  1:1264692/1‑93 (MQ=255)

                                         cATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAc                                                <  1:6786092/93‑1 (MQ=255)

                                         cATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAc                                                <  2:3070168/93‑1 (MQ=255)

                                          atatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  1:1009006/1‑93 (MQ=255)

                                           tatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGa                                              >  2:2213920/1‑93 (MQ=255)

                                           tatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGa                                              <  1:4067300/93‑1 (MQ=255)

                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             <  2:8921563/93‑1 (MQ=255)

                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             <  1:8885355/93‑1 (MQ=255)

                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             <  2:5732625/93‑1 (MQ=255)

                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:4848089/1‑93 (MQ=255)

                                              ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           >  2:8692806/1‑93 (MQ=255)

                                               tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          >  1:1241901/1‑93 (MQ=255)

                                                ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         >  1:1456490/1‑93 (MQ=255)

                                                ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         >  1:5028820/1‑93 (MQ=255)

                                                  ctCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       <  1:1631729/92‑1 (MQ=255)

                                                  aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:8785333/1‑93 (MQ=255)

                                                  aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  2:4501989/1‑93 (MQ=255)

                                                   tCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaac                                      <  2:4716003/93‑1 (MQ=255)

                                                     cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    <  1:4874412/93‑1 (MQ=255)

                                                     cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  1:6217799/1‑93 (MQ=255)

                                                       gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACACgg                                  >  1:3482319/1‑93 (MQ=255)

                                                        aGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGt                                 <  2:3314178/93‑1 (MQ=255)

                                                          tGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaa                               >  2:8766289/1‑93 (MQ=255)

                                                          tGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaa                               >  1:6320028/1‑93 (MQ=255)

                                                           gACGCGCTAAGCGACCCTTCCCGTCTGCAAAGACTGCCGCATTGTGAGCAGCTGGTTTAGCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              <  2:5135066/93‑1 (MQ=255)

                                                           gACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              >  2:2643334/1‑93 (MQ=255)

                                                            aCGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAag                             <  2:6572507/93‑1 (MQ=255)

                                                            aCGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAag                             <  1:5808731/93‑1 (MQ=255)

                                                             cgcgCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAaga                            <  2:997196/93‑1 (MQ=255)

                                                              gcgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAagag                           >  1:9272976/1‑93 (MQ=255)

                                                              gcgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAagag                           <  1:5205168/93‑1 (MQ=255)

                                                               cgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAgg                          >  2:1840760/1‑93 (MQ=255)

                                                                gcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGt                         <  2:7423597/93‑1 (MQ=255)

                                                                 cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        >  2:9634665/1‑93 (MQ=255)

                                                                 cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        <  1:5377277/93‑1 (MQ=255)

                                                                  tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       <  1:8476980/93‑1 (MQ=255)

                                                                  tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  1:8391729/1‑93 (MQ=255)

                                                                    aGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACg                     <  1:3000248/93‑1 (MQ=255)

                                                                     gCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGGTACgg                    >  1:4191345/1‑93 (MQ=255)

                                                                     gCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGGAAAGAGGTTACgg                    >  1:7682708/1‑93 (MQ=255)

                                                                      cGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGc                   <  1:802860/93‑1 (MQ=255)

                                                                           cTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAc              <  1:1955566/93‑1 (MQ=255)

                                                                           cTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAc              >  1:4404357/1‑93 (MQ=255)

                                                                           cTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAc              <  2:7809663/93‑1 (MQ=255)

                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGATTACGGCGAAAcc             <  1:8006723/93‑1 (MQ=255)

                                                                               ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          >  2:2397722/1‑93 (MQ=255)

                                                                               ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  2:3715852/93‑1 (MQ=255)

                                                                               ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          >  2:815188/1‑93 (MQ=255)

                                                                                cGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGcc         >  2:1072593/1‑93 (MQ=255)

                                                                                  tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       <  2:6998789/93‑1 (MQ=255)

                                                                                  tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       >  1:5413995/1‑93 (MQ=255)

                                                                                  tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       <  1:7623715/93‑1 (MQ=255)

                                                                                   cTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCCGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGt      <  2:4224554/93‑1 (MQ=255)

                                                                                   cTGCAAAGACTGCCGCATTTTGAGCAGCTGGCTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACCGCGAAAACTGccagg      >  2:280816/1‑90 (MQ=255)

                                                                                    tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     >  1:7501726/1‑93 (MQ=255)

                                                                                     gCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGa    >  1:3943885/1‑93 (MQ=255)

                                                                                     gCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGa    >  2:1525616/1‑93 (MQ=255)

                                                                                     gCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGa    <  1:6446391/93‑1 (MQ=255)

                                                                                     gCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGa    >  2:9967993/1‑93 (MQ=255)

                                                                                      cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGCAGTTACGGCGAACCCTGCCGGAGAc   >  2:5224514/1‑93 (MQ=255)

                                                                                      cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   >  2:9059187/1‑93 (MQ=255)

                                                                                      cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGATGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   <  1:8441912/93‑1 (MQ=255)

TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992

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