BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,059,421	(A)_7→8	coding (426/1644 nt)	yeeO ←	FMN/FAD exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,059,414	1	.	A	99.7%	1286.4 / ‑3.2	296	coding (433/1644 nt)	yeeO	FMN/FAD exporter
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (146/149); total (147/149)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 343 total aligned reads displayed.

ATGATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑AAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAA  >  NC_000913/2059314‑2059509

aTGATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATagcagc                                                                                                  >  2:3769057/1‑101 (MQ=255)

 tGATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGCCAGTAGTACCAAGATCGCTAgcagc‑a                                                                                                <  1:3768499/101‑1 (MQ=255)

 tGATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAgcagc‑a                                                                                                <  1:3768575/101‑1 (MQ=255)

  gATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aa                                                                                               >  1:3769060/1‑101 (MQ=255)

   aTCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaa                                                                                              >  2:3769062/1‑101 (MQ=255)

     cACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaa                                                                                            <  2:3768332/101‑1 (MQ=255)

     cACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaa                                                                                             <  1:3768542/100‑1 (MQ=255)

      aCCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaa                                                                                           >  1:3769069/1‑101 (MQ=255)

       ccAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaa                                                                                          >  2:3769070/1‑101 (MQ=255)

         aaTGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaaag                                                                                        <  2:3768674/101‑3 (MQ=255)

          aTGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaaagc                                                                                       <  2:3768414/101‑4 (MQ=255)

           tGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaaagcc                                                                                      >  2:3769073/1‑97 (MQ=255)

            gACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaaagcca                                                                                     >  1:3769077/1‑96 (MQ=255)

            gACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaaaaagcca                                                                                     <  2:3768467/101‑6 (MQ=255)

               tGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATaa                                                                                   >  2:3769080/1‑101 (MQ=255)

               tGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATaa                                                                                   <  1:3768729/101‑1 (MQ=255)

                 ccGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATg                                                                                 <  1:3768403/101‑1 (MQ=255)

                 ccGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATg                                                                                 >  2:3769085/1‑101 (MQ=255)

                   gcgTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGAc                                                                               <  2:3768621/101‑1 (MQ=255)

                   gcgTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGAc                                                                               <  2:3768421/101‑1 (MQ=255)

                    cgTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGAcc                                                                              <  2:3768755/101‑1 (MQ=255)

                     gTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCa                                                                             >  1:3769087/1‑101 (MQ=255)

                      tCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCat                                                                            <  2:3768589/101‑1 (MQ=255)

                      tCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCat                                                                            <  1:3768599/101‑1 (MQ=255)

                      tCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑aaaaa                                                                                            <  2:3768198/84‑1 (MQ=255)

                         ccACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATAtt                                                                         >  2:3769094/1‑101 (MQ=255)

                          cACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTg                                                                        >  1:3769097/1‑101 (MQ=255)

                           aCCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCata                                                                           >  2:3769103/1‑97 (MQ=255)

                           aCCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGa                                                                       >  2:3769099/1‑101 (MQ=255)

                            cccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGaa                                                                      <  1:3768646/101‑1 (MQ=255)

                            cccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGa                                                                       <  2:3768554/100‑1 (MQ=255)

                             ccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAg                                                                     >  1:3769105/1‑101 (MQ=255)

                             ccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAg                                                                     <  2:3768773/101‑1 (MQ=255)

                             ccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAg                                                                     <  1:3768534/101‑1 (MQ=255)

                             ccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAg                                                                     <  1:3768638/101‑1 (MQ=255)

                             ccTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAg                                                                     <  1:3768583/101‑1 (MQ=255)

                                 gctcgTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGt                                                                 >  1:3769110/1‑101 (MQ=255)

                                    cgtcgtcgATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCg                                                              >  1:3769115/1‑101 (MQ=255)

                                     gtcgtcgATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGc                                                             >  2:3769116/1‑101 (MQ=255)

                                      tcgtcgATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGcc                                                            <  1:3768349/101‑1 (MQ=255)

                                      tcgtcgATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGcc                                                            <  1:3768468/101‑1 (MQ=255)

                                        gtcgATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCaa                                                          <  1:3768524/101‑1 (MQ=255)

                                          cgATCCCGCTTACCGAGACTAAATGCCACAACGCCAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATc                                                        >  1:3769122/1‑101 (MQ=255)

                                           gATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATcc                                                       >  2:3769130/1‑101 (MQ=255)

                                           gATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATcc                                                       >  2:3769124/1‑101 (MQ=255)

                                           gATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATcc                                                       >  1:3769129/1‑101 (MQ=255)

                                             tCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCa                                                     >  2:3769133/1‑101 (MQ=255)

                                             tCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCa                                                     <  2:3768962/101‑1 (MQ=255)

                                             tCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCa                                                     <  2:3768711/101‑1 (MQ=255)

                                              cccGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCAc                                                    >  2:3769137/1‑101 (MQ=255)

                                                cGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGc                                                  <  1:3768688/101‑1 (MQ=255)

                                                cGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGc                                                  >  2:3769140/1‑101 (MQ=255)

                                                 gCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGcc                                                 <  1:3768782/101‑1 (MQ=255)

                                                  cTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCg                                                <  1:3768933/101‑1 (MQ=255)

                                                  cTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCg                                                <  1:3768517/101‑1 (MQ=255)

                                                   ttACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCgg                                               <  1:3768280/101‑1 (MQ=255)

                                                   ttACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCgg                                               >  2:3769144/1‑101 (MQ=255)

                                                   ttACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCgg                                               <  1:3768757/101‑1 (MQ=255)

                                                      ccGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCa                                            >  1:3769151/1‑101 (MQ=255)

                                                        gagaCTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATc                                          >  2:3769153/1‑101 (MQ=255)

                                                          gaCTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGc                                        >  2:3769156/1‑101 (MQ=255)

                                                          gaCTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGc                                        >  1:3769158/1‑101 (MQ=255)

                                                           aCTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCa                                       <  1:3768829/101‑1 (MQ=255)

                                                            cTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCAt                                      <  1:3768795/101‑1 (MQ=255)

                                                            cTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCAt                                      >  2:3769162/1‑101 (MQ=255)

                                                            cTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCAt                                      <  2:3768504/101‑1 (MQ=255)

                                                              aaaTGCCACAACGACAGTAGTACCAAGATCGCTAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTg                                                                  >  2:3769166/1‑71 (MQ=255)

                                                                aTGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCttt                                  <  1:3768316/101‑1 (MQ=255)

                                                                 tGCCACAACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCtttt                                 >  1:3769171/1‑101 (MQ=255)

                                                                     acaACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCcca                             <  2:3768794/101‑1 (MQ=255)

                                                                     acaACGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCcca                             >  2:3769174/1‑101 (MQ=255)

                                                                       aaCGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagc                           <  2:3768692/101‑1 (MQ=255)

                                                                       aaCGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagc                           >  2:3769179/1‑101 (MQ=255)

                                                                       aaCGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagc                           >  1:3769178/1‑101 (MQ=255)

                                                                        aCGACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagcc                          >  2:3769183/1‑101 (MQ=255)

                                                                          gACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagccag                        >  1:3769186/1‑101 (MQ=255)

                                                                          gACAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCccagccag                        <  2:3768910/101‑1 (MQ=255)

                                                                           aCAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCcagccagc                       >  1:3769193/1‑101 (MQ=255)

                                                                           aCAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCcagccagc                       >  2:3769192/1‑101 (MQ=255)

                                                                           aCAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCcagccagc                       <  1:3768559/101‑1 (MQ=255)

                                                                            cAGTAGTACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCt                      <  1:3768669/101‑1 (MQ=255)

                                                                             agtagtACCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTg                     >  1:3769197/1‑101 (MQ=255)

                                                                              gtagtaCCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTg                     >  2:3769199/1‑100 (MQ=255)

                                                                              gtagtaCCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGa                    <  2:3768838/101‑1 (MQ=255)

                                                                                 gtaCCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCa                 >  2:3769203/1‑101 (MQ=255)

                                                                                  taCCAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAg                >  2:3769207/1‑101 (MQ=255)

                                                                                    ccAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGaa              <  1:3769019/101‑1 (MQ=255)

                                                                                     cAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGaaa             >  1:3769212/1‑101 (MQ=255)

                                                                                     cAAGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGaa              <  1:3768611/100‑1 (MQ=255)

                                                                                      aaGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGaaaa            >  1:3769216/1‑101 (MQ=255)

                                                                                      aaGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGaaaa            >  2:3769218/1‑101 (MQ=255)

                                                                                       aGATCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAg           <  1:3768604/101‑1 (MQ=255)

                                                                                         aTCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTg         <  1:3768846/101‑1 (MQ=255)

                                                                                          tCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGc        <  1:3768675/101‑1 (MQ=255)

                                                                                          tCGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGc        >  1:3769222/1‑101 (MQ=255)

                                                                                           cGATAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCAAGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCt       >  2:3769225/1‑101 (MQ=255)

                                                                                              tAGCAGCAAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAg    <  2:3768636/101‑1 (MQ=255)

                                                                                               agcagcaAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGa   <  2:3768658/101‑1 (MQ=255)

                                                                                               agcagcaAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGa   <  2:3768603/101‑1 (MQ=255)

                                                                                                gcagcaAAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGaa  <  1:3768183/101‑1 (MQ=255)

ATGATCACCAATGACTGCCGCGTCGCCACCCTCGCTCGTCGTCGATCCCGCTTACCGAGACTAAATGCCACAACGACAGTAGTACCAAGATCGATAGCAGC‑AAAAAAAGCCATAATGACCATATTGAAGCTGTCCGCCAATCCCACGCCGGCCATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAA  >  NC_000913/2059314‑2059509

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