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%	795.3 / NA	239	K138R (AAA→AGA)	b1757	putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (125/114); total (125/114)

Only 100 of 244 total aligned reads displayed.

TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992

tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         <  1:4820988/93‑1 (MQ=255)

tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         <  1:5542444/93‑1 (MQ=255)

tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         <  2:2167053/93‑1 (MQ=255)

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

  gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       <  1:1076578/93‑1 (MQ=255)

  gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       >  1:4631279/1‑93 (MQ=255)

  gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       <  2:1710097/93‑1 (MQ=255)

   aCGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGc                                                                                      <  2:2367704/93‑1 (MQ=255)

   aCGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGc                                                                                      <  1:5451800/93‑1 (MQ=255)

   aCGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGc                                                                                      >  1:3737800/1‑93 (MQ=255)

    cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGcc                                                                                     >  1:4199710/1‑93 (MQ=255)

    cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGcc                                                                                     >  2:1271327/1‑93 (MQ=255)

       cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  <  2:3255565/93‑1 (MQ=255)

       cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  <  1:63578/93‑1 (MQ=255)

         tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                <  1:5708212/93‑1 (MQ=255)

         tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAACGAACCTTCCCGTCTGCAACGACTGCCGCAtt                                                                                >  1:2419348/1‑93 (MQ=255)

           aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              <  1:1625938/93‑1 (MQ=255)

              aCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGag                                                                           >  2:581653/1‑93 (MQ=255)

              aCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGag                                                                           >  1:782010/1‑93 (MQ=255)

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

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

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

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

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

                      gCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGt                                                                   <  1:1687039/93‑1 (MQ=255)

                        aaaaaaCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGttt                                                                 <  1:5294707/93‑1 (MQ=255)

                        aaaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGttt                                                                 <  1:4512754/93‑1 (MQ=255)

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

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

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

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

                           aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  1:1840629/93‑1 (MQ=255)

                            aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  1:3438628/93‑1 (MQ=255)

                            aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  1:2546442/93‑1 (MQ=255)

                            aGCAGGTTTAACACATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  2:4323787/93‑1 (MQ=255)

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

                                  tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       >  2:4169521/1‑93 (MQ=255)

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

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

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

                                    tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  2:860421/93‑1 (MQ=255)

                                     aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    >  2:5647084/1‑93 (MQ=255)

                                     aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    >  2:1339330/1‑93 (MQ=255)

                                     aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    >  2:5067888/1‑93 (MQ=255)

                                        gCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCa                                                 >  2:3958884/1‑93 (MQ=255)

                                        gCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCa                                                 >  2:6243698/1‑93 (MQ=255)

                                        gCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCa                                                 >  2:191814/1‑93 (MQ=255)

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

                                          atatCTCCATCCTGAGTGACACGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  2:2509964/1‑93 (MQ=255)

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

                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             >  2:165656/1‑93 (MQ=255)

                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             >  2:1057409/1‑93 (MQ=255)

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

                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            <  1:1337407/93‑1 (MQ=255)

                                              ctcCATCCTGAGTGACGCGCTAGGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  1:694845/93‑1 (MQ=255)

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

                                              ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  2:4528560/93‑1 (MQ=255)

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

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

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

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

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

                                                       gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  2:6551633/1‑93 (MQ=255)

                                                       gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  2:4831469/1‑93 (MQ=255)

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

                                                        aGTCACGCGCTAAGCGAACGTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGt                                 <  1:406930/93‑1 (MQ=255)

                                                         gTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTa                                <  2:3342813/93‑1 (MQ=255)

                                                         gTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTa                                >  1:4379679/1‑93 (MQ=255)

                                                         gTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTa                                <  1:3180474/93‑1 (MQ=255)

                                                           gACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              <  2:5965340/93‑1 (MQ=255)

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

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

                                                               cgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAgg                          <  2:1243405/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                   aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  1:4992678/93‑1 (MQ=255)

                                                                   aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      >  2:366086/1‑93 (MQ=255)

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

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

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

                                                                        aaCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGa                 >  1:1349425/1‑93 (MQ=255)

                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             >  2:2330613/1‑93 (MQ=255)

                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  2:4966942/93‑1 (MQ=255)

                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             >  2:318348/1‑93 (MQ=255)

                                                                             tCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCt            >  2:5254834/1‑93 (MQ=255)

                                                                              cccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAAGGCTGCACGACCTGCAACAAGGTGAAGAGGTTACGGCGAAACCTg           <  2:753964/93‑1 (MQ=255)

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

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

                                                                                   cTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGt      >  1:1307385/1‑93 (MQ=255)

                                                                                    tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     <  2:5905648/93‑1 (MQ=255)

                                                                                    tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     <  1:619681/93‑1 (MQ=255)

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

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

                                                                                      cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   <  2:1771564/93‑1 (MQ=255)

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

                                                                                       aaagaCTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  >  1:3559331/1‑93 (MQ=255)

TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992

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