BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,790,787	C→T	5.6%	S153N (AGC→AAC)	ydiU ←	UPF0061 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,790,787	0	C	T	5.6%	332.4 / 5.0	105	S153N (AGC→AAC)	ydiU	UPF0061 family protein
Reads supporting (aligned to +/- strand): ref base C (43/56); new base T (3/3); total (46/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCC  >  NC_000913/1790638‑1790934

gACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGc                                                                                                                                                     <  1:1488790/150‑1 (MQ=255)

  cGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTATAAATGTTCGAAATCACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTCGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTg                                                                                                                                                   >  1:1530457/1‑150 (MQ=255)

  cGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCTGATACACTGGCGAATCGCTg                                                                                                                                                   >  1:76568/1‑150 (MQ=255)

  cGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTg                                                                                                                                                   <  2:1252347/150‑1 (MQ=255)

  cGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGGAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATAAACTGGCGAATCGCTg                                                                                                                                                   >  1:39272/1‑150 (MQ=255)

   taaTAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTgg                                                                                                                                                  <  1:830130/149‑1 (MQ=255)

   gAATAGCAAATTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAATCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTgg                                                                                                                                                  <  1:137875/150‑1 (MQ=255)

   gAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTgg                                                                                                                                                  <  1:1153050/150‑1 (MQ=255)

        gCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCGGGTGACg                                                                                                                                             <  2:857212/150‑1 (MQ=255)

        gCAAAGTCAGCTAACTGACGAACCTATTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGAGCAGATGAATAGGTGCCACACGCATCAGCATCGCTCATGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACg                                                                                                                                             >  1:226848/1‑150 (MQ=255)

           aaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGCCCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACCCTGGCGAATCGCTGGTGACGATa                                                                                                                                          >  1:581993/1‑150 (MQ=255)

            aGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATAc                                                                                                                                         <  1:997913/150‑1 (MQ=255)

                 gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAc                                                                                                                                    >  1:1690642/1‑150 (MQ=255)

                  cTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcg                                                                                                                                   >  2:359897/1‑150 (MQ=255)

                      cTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgcg                                                                                                                               >  2:1003520/1‑150 (MQ=255)

                       tGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAACGCGCAGATGACTTGGTGCCACACCCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACACGCgg                                                                                                                              >  2:763873/1‑150 (MQ=255)

                       tGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCgg                                                                                                                              >  1:137246/1‑150 (MQ=255)

                        gACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAACGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCggg                                                                                                                             <  1:436078/150‑1 (MQ=255)

                        gACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCggg                                                                                                                             >  1:388672/1‑150 (MQ=255)

                         aCGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAACGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgt                                                                                                                            <  1:1471644/150‑1 (MQ=255)

                          cGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCCCACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGCCGATACTTAACGCGCGGgtc                                                                                                                           >  2:485312/1‑150 (MQ=255)

                            aaCCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGCCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGCATCGCTGGTGACGATACTTAACGCGCGGgtcgt                                                                                                                         >  2:308570/1‑150 (MQ=255)

                            aaCCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgt                                                                                                                         <  1:1724786/150‑1 (MQ=255)

                               cTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg                                                                                                                      <  2:770776/150‑1 (MQ=255)

                               cTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg                                                                                                                      <  2:438287/150‑1 (MQ=255)

                                 tttCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAATGCGCTGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGACGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGaa                                                                                                                    <  2:80889/150‑1 (MQ=255)

                                     cGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc                                                                                                                <  2:1469552/150‑1 (MQ=255)

                                       aCTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCa                                                                                                              >  2:1181168/1‑150 (MQ=255)

                                        cTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAg                                                                                                             <  1:1563320/150‑1 (MQ=255)

                                        cTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATCCACAg                                                                                                             >  2:812329/1‑150 (MQ=255)

                                          ctcggcggTAGTAAAAATGTTCGAACTGACCAAAGCGCAGTTGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTTGTGACGTTACTTAACGCGCGGGTCGTCGGAATGCCCAGAt                                                                                                           <  1:1295582/148‑1 (MQ=255)

                                          cgcgGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACCCGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGAt                                                                                                           >  1:1140810/1‑150 (MQ=255)

                                             ggcggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCCTCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAAt                                                                                                        >  2:1556196/1‑150 (MQ=255)

                                              gcggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATg                                                                                                       <  2:1086507/150‑1 (MQ=255)

                                                    gtaAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGc                                                                                                 >  2:574627/1‑150 (MQ=255)

                                                     taAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACAGGCATCAGCAGCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGcc                                                                                                <  2:1533521/150‑1 (MQ=255)

                                                        aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTTGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCa                                                                                             <  2:1600503/150‑1 (MQ=255)

                                                        aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCa                                                                                             <  2:651089/150‑1 (MQ=255)

                                                           tGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTg                                                                                          <  2:564733/150‑1 (MQ=255)

                                                              tCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCg                                                                                       <  1:37334/150‑1 (MQ=255)

                                                              tCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCg                                                                                       <  1:557840/150‑1 (MQ=255)

                                                               cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCCTCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCga                                                                                      <  2:507414/150‑1 (MQ=255)

                                                               cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCga                                                                                      <  2:137246/150‑1 (MQ=255)

                                                                gAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGTTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCgag                                                                                     >  2:733822/1‑150 (MQ=255)

                                                                gAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGTTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCgag                                                                                     >  2:734983/1‑150 (MQ=255)

                                                                gAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGTTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCgag                                                                                     >  2:740249/1‑150 (MQ=255)

                                                                gAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCgag                                                                                     >  2:974977/1‑150 (MQ=255)

                                                                gAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCgag                                                                                     <  1:330255/150‑1 (MQ=255)

                                                                      aCCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGACTGCCCAGATAATGCATCGCCTCACTGGCGAGACTTtc                                                                               >  1:1528864/1‑150 (MQ=255)

                                                                      aCCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAAGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTtc                                                                               >  2:1287643/1‑150 (MQ=255)

                                                                        cAAAGCGCCGATGACTTGGTGCCAAACGCATCATCATCGCGCCTGGCTCCGCAGTTTCCGGATACACTCGCGAATCGCGGGTGACGATACTTCACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTtctc                                                                             <  2:116842/150‑1 (MQ=255)

                                                                        cAAAGCGCAGATGACTTTGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGATTCCCGATACACTGGCGAATCGCTGGTGACGATACTTCCCGCGCGGGTCGTCGGAATGCCCAGA‑AATGCATCGCCTCACTGGCGAGACTTTCTCg                                                                            <  2:85025/150‑1 (MQ=255)

                                                                         acaGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGAGACGATACATAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCg                                                                            <  1:264235/148‑1 (MQ=255)

                                                                         aaaGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGGGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCg                                                                            <  1:1805768/150‑1 (MQ=255)

                                                                               cAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCCCTGGCGAGACTGTCTCGTAtcgt                                                                      >  1:804448/1‑150 (MQ=255)

                                                                                  aTGACTTGGAGCCACACGCATCAGCATCGCGCCTGGCTCAGCCGTTTCCCTATACACTGGCGAATCTCTTGTGACGATCCTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGa                                                                   <  1:1504392/150‑1 (MQ=255)

                                                                                     aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATTCACTGGCGAATCGCTGGTGACGAAAGTTAACGCGCGGCTCGTCGGCCTGCCCAGATACTGCCTCGCCGCCCTGGCGAGCCTTACTCGAATCGTCGCACg                                                                >  1:100089/1‑150 (MQ=255)

                                                                                      cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTACCGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGt                                                               >  1:1619582/1‑150 (MQ=255)

                                                                                      cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGATTGCCCAGATAACGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGt                                                               >  1:145994/1‑150 (MQ=255)

                                                                                      cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGt                                                               >  1:1240710/1‑150 (MQ=255)

                                                                                        tGGTGACACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTTGTGACGATACTTAAGGCGCGGGTCGTCGGAATGACCAGATAATGCATCGCCTCTCTGGCGTGACTTTCTCGTTTCGTGGAACGTat                                                             >  1:889214/1‑149 (MQ=255)

                                                                                          gTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTaaaa                                                           >  2:1064/1‑150 (MQ=255)

                                                                                             ccACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCg                                                        <  2:131310/150‑1 (MQ=255)

                                                                                              cacaCGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCGAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGc                                                       <  1:1335185/150‑1 (MQ=255)

                                                                                              cacaCGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGc                                                       >  1:1665505/1‑150 (MQ=255)

                                                                                                cacGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGGCGTCGGAATGCCCAGATAATGCAGCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGccc                                                     <  1:1139850/150‑1 (MQ=255)

                                                                                                 acGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCg                                                    <  2:478278/150‑1 (MQ=255)

                                                                                                   gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTc                                                  <  2:1351063/150‑1 (MQ=255)

                                                                                                   gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTc                                                  <  2:1530457/150‑1 (MQ=255)

                                                                                                   gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTc                                                  <  1:729562/150‑1 (MQ=255)

                                                                                                     aTCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGACGTCGGAATGCCCAGATAATGCATCGCCTCCCTGGCGAGACTTTCTCGTATCGTCGAACGTAACACCGCCCGTcaa                                                >  1:347405/1‑148 (MQ=255)

                                                                                                       cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATc                                              <  2:169927/150‑1 (MQ=255)

                                                                                                       cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCACACTGGCGAGACTGACTCGTATCGTCGAACGTAAAACCGCCCGTCCATc                                              >  2:112475/1‑150 (MQ=255)

                                                                                                        aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGCACGTAAAACCGCCCGTCCATcc                                             >  1:1386499/1‑149 (MQ=255)

                                                                                                          cATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTCGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCAcc                                           <  2:1556565/150‑1 (MQ=255)

                                                                                                           aTCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACGGTCTCGTATCGGCGACCGTAAAACCGCCCGTCCATCAccc                                          >  1:1750889/1‑150 (MQ=255)

                                                                                                             cgcgCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCAt                                        >  2:1678641/1‑150 (MQ=255)

                                                                                                                   gggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGTACGTAAAACCGCCCGTCCATCACCCATTCGCGa                                  <  1:1427516/149‑1 (MQ=255)

                                                                                                                    ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGaa                                 <  2:850894/150‑1 (MQ=255)

                                                                                                                     gCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAAt                                >  1:1415667/1‑150 (MQ=255)

                                                                                                                      cTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATa                               >  2:808485/1‑150 (MQ=255)

                                                                                                                         cgccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTCTCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAgg                            >  2:1544216/1‑150 (MQ=255)

                                                                                                                           ccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGCCTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCg                          >  1:894868/1‑150 (MQ=255)

                                                                                                                                tCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAgg                     <  1:1181168/150‑1 (MQ=255)

                                                                                                                                  ccGATACACTGGCGACTCGCGGGGGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTAGCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGcc                   <  1:893877/150‑1 (MQ=255)

                                                                                                                                  ccGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGcc                   >  1:1758746/1‑150 (MQ=255)

                                                                                                                                  ccGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGcc                   <  2:1690642/150‑1 (MQ=255)

                                                                                                                                   cGATACACTGGAGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGATTAAGGCGTCAGGCCa                  <  2:1223808/150‑1 (MQ=255)

                                                                                                                                     aTACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGc                <  1:1809820/150‑1 (MQ=255)

                                                                                                                                          cTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCtt           >  2:584482/1‑150 (MQ=255)

                                                                                                                                            ggCGAATCGCTGGTGACGATACTTAACGCGCGGTGCGTCGTAATGCCCTGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTc         <  2:226848/150‑1 (MQ=255)

                                                                                                                                            ggCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGCTCGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCttta         >  2:1584733/1‑149 (MQ=255)

                                                                                                                                            ggCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTTAAACCGCCCGGCCATGACCCCTTCGGGAAGAAGGCGGCAGGCCCGCACCtttg         >  2:906998/1‑149 (MQ=255)

                                                                                                                                             gCGAATCGTTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCATCAGGCCAGCACCTTTCa        <  1:733822/150‑1 (MQ=255)

                                                                                                                                             gCGAATCGTTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCATCAGGCCAGCACCTTTCa        <  1:740249/150‑1 (MQ=255)

                                                                                                                                             gCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCa        <  1:1669514/150‑1 (MQ=255)

                                                                                                                                                 aTCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATg    <  1:686255/150‑1 (MQ=255)

                                                                                                                                                   cGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGcc  <  2:763854/150‑1 (MQ=255)

                                                                                                                                                   cGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGcc  <  1:308570/150‑1 (MQ=255)

                                                                                                                                                   cGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGcc  <  2:1531412/150‑1 (MQ=255)

GACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCC  >  NC_000913/1790638‑1790934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`