BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	433,889	G→A	6.6%	V438V (GTC→GTT)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	433,889	0	G	A	6.6%	698.4 / 26.1	227	V438V (GTC→GTT)	glcB	malate synthase G
Reads supporting (aligned to +/- strand): ref base G (119/93); new base A (10/5); total (129/98)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 227 total aligned reads displayed.

TCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTACACGCTGC  >  NC_002947/433744‑434037

   aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                                                               >  1:4549856/1‑150 (MQ=255)

   aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATAGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                                                               <  1:2265192/150‑1 (MQ=255)

       ctcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACg                                                                                                                                           >  1:1409251/1‑150 (MQ=255)

        tcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACgg                                                                                                                                          >  1:1664916/1‑150 (MQ=255)

          gTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCg                                                                                                                                        >  1:4352148/1‑150 (MQ=255)

             cttcATGTCGCCCTTGCGGATCATCGGGCCAGCTTCCATCGAGGTATGGATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                                                                     <  1:5436363/150‑1 (MQ=255)

                cATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTc                                                                                                                                  >  1:2964071/1‑150 (MQ=255)

                   ggCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTc                                                                                                                               <  1:1811072/150‑1 (MQ=255)

                   ggCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTc                                                                                                                               >  2:847357/1‑150 (MQ=255)

                        cTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatga                                                                                                                          <  1:706354/150‑1 (MQ=255)

                         tttGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatgat                                                                                                                         >  2:4944688/1‑150 (MQ=255)

                            gcgcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGcc                                                                                                                      <  1:4944659/150‑1 (MQ=255)

                                accacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGAcc                                                                                                                  <  2:1388522/150‑1 (MQ=255)

                                 ccacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCt                                                                                                                 <  2:1822831/150‑1 (MQ=255)

                                    cGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGa                                                                                                              <  2:4087711/150‑1 (MQ=255)

                                       cgcCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACgg                                                                                                           <  1:2028319/150‑1 (MQ=255)

                                         ccGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtg                                                                                                         <  1:5251070/150‑1 (MQ=255)

                                          cGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtgt                                                                                                        >  2:4939896/1‑150 (MQ=255)

                                          cGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtgt                                                                                                        <  1:1083657/150‑1 (MQ=255)

                                           ggCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGtt                                                                                                       <  2:944505/150‑1 (MQ=255)

                                            gCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTg                                                                                                      >  1:3860075/1‑150 (MQ=255)

                                              ttCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTgcg                                                                                                    >  2:4399202/1‑150 (MQ=255)

                                              ttCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTgcg                                                                                                    >  2:1949076/1‑150 (MQ=255)

                                                 cATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCgg                                                                                                 >  1:2804632/1‑150 (MQ=255)

                                                  aTCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGc                                                                                                >  2:2776986/1‑150 (MQ=255)

                                                     gAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATg                                                                                             >  2:4205425/1‑150 (MQ=255)

                                                      aGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGc                                                                                            <  2:1379755/150‑1 (MQ=255)

                                                       ggTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                                                           <  2:4352148/150‑1 (MQ=255)

                                                        gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                                                           <  1:533572/149‑1 (MQ=255)

                                                        gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                                                           >  2:533572/1‑149 (MQ=255)

                                                            ggATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCaagca                                                                                      <  1:5369012/150‑1 (MQ=255)

                                                            ggATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcagca                                                                                      <  2:515348/150‑1 (MQ=255)

                                                             gATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcagcag                                                                                     >  1:3923005/1‑150 (MQ=255)

                                                               tttCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAggc                                                                                   >  2:1120062/1‑149 (MQ=255)

                                                               tttCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGt                                                                                   >  2:5433183/1‑150 (MQ=255)

                                                                ttCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                                                                     >  1:3288647/1‑99 (MQ=255)

                                                                ttCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                                                                     <  2:3288647/99‑1 (MQ=255)

                                                                 tCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCt                                                                                 >  1:1189535/1‑150 (MQ=255)

                                                                      gCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGAc                                                                            <  1:481001/150‑1 (MQ=255)

                                                                      gCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggt                                                                            >  2:2346006/1‑148 (MQ=255)

                                                                       ccAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggtg                                                                           >  1:5025322/1‑147 (MQ=255)

                                                                          gTGCGGTCAAGGAAGCCGGGGTTGCTGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTgg                                                                        <  1:2581345/150‑1 (MQ=255)

                                                                            gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                                                                      <  1:3478759/150‑1 (MQ=255)

                                                                              ggTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGa                                                                    <  1:1234482/150‑1 (MQ=255)

                                                                                tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                                                                  >  1:4128625/1‑150 (MQ=255)

                                                                                tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                                                                  >  2:5343243/1‑150 (MQ=255)

                                                                                tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                                                                  >  1:4388050/1‑150 (MQ=255)

                                                                                  aaGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAt                                                                >  1:4066354/1‑150 (MQ=255)

                                                                                    ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGtt                                                                                                       >  1:1001494/1‑109 (MQ=255)

                                                                                    ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                                                              >  1:2716422/1‑150 (MQ=255)

                                                                                    ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                                                              >  2:364522/1‑150 (MQ=255)

                                                                                    ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                                                              >  2:796833/1‑150 (MQ=255)

                                                                                     gAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtctc                                                             >  1:2410380/1‑150 (MQ=255)

                                                                                        gCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGt                                                          >  2:3706390/1‑150 (MQ=255)

                                                                                          cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACg                                                                                                            <  1:1768360/98‑1 (MQ=255)

                                                                                          cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACg                                                                                                            >  2:1768360/1‑98 (MQ=255)

                                                                                          cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTgg                                                        <  2:4762382/150‑1 (MQ=255)

                                                                                          cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTgg                                                        <  2:2621375/150‑1 (MQ=255)

                                                                                           ggTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGc                                                       <  1:4004692/150‑1 (MQ=255)

                                                                                            gtgtTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCg                                                      >  2:5007949/1‑150 (MQ=255)

                                                                                                 gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc                                                 >  2:2810356/1‑150 (MQ=255)

                                                                                                  atgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGcc                                                <  1:4606711/150‑1 (MQ=255)

                                                                                                   tgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCa                                               >  2:1523162/1‑150 (MQ=255)

                                                                                                    gaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCAc                                              >  2:559144/1‑150 (MQ=255)

                                                                                                     aaCACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACt                                             <  1:1466788/150‑1 (MQ=255)

                                                                                                      acacCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGTTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACtt                                            >  2:1689516/1‑150 (MQ=255)

                                                                                                       caccacGCTCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTc                                           >  1:4300007/1‑150 (MQ=255)

                                                                                                         ccacGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCt                                         <  1:254494/150‑1 (MQ=255)

                                                                                                          cacGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTc                                        >  2:2839141/1‑150 (MQ=255)

                                                                                                          cacGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTc                                        >  1:3963284/1‑150 (MQ=255)

                                                                                                          cacGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTc                                        >  2:4012005/1‑150 (MQ=255)

                                                                                                            cgcgCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAg                                      >  1:594088/1‑150 (MQ=255)

                                                                                                               gcTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGc                                   >  1:1294477/1‑150 (MQ=255)

                                                                                                                cTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGcc                                  <  1:249233/150‑1 (MQ=255)

                                                                                                                 tCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCg                                 >  2:3192999/1‑150 (MQ=255)

                                                                                                                 tCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCg                                 >  1:1254872/1‑150 (MQ=255)

                                                                                                                 tCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCg                                 >  2:3953632/1‑150 (MQ=255)

                                                                                                                  cggcggCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGt                                <  2:2770906/150‑1 (MQ=255)

                                                                                                                  cggcggCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGt                                <  2:2292340/150‑1 (MQ=255)

                                                                                                                   ggcggcTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTg                               >  2:158781/1‑150 (MQ=255)

                                                                                                                    gcggcTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGc                              <  1:4182858/150‑1 (MQ=255)

                                                                                                                          gCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTc                        <  2:4549856/150‑1 (MQ=255)

                                                                                                                           ccTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCg                       >  1:3816650/1‑150 (MQ=255)

                                                                                                                             ttGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtg                                                                                                         <  2:1916609/66‑1 (MQ=255)

                                                                                                                             ttGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGc                     <  1:929530/150‑1 (MQ=255)

                                                                                                                             ttGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACATTGACGgtg                                                                                                         >  1:1916609/1‑66 (MQ=255)

                                                                                                                               gATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCtt                   >  2:2484785/1‑150 (MQ=255)

                                                                                                                               gATGCAGGACTTGAGGTTGACCGTGGTACGGCGGTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCtt                   <  2:3610532/150‑1 (MQ=255)

                                                                                                                                  gCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGAc                >  2:3130116/1‑150 (MQ=255)

                                                                                                                                   cAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACg               >  1:3142613/1‑150 (MQ=255)

                                                                                                                                    aGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGa              >  2:5245056/1‑150 (MQ=255)

                                                                                                                                       aCTTGAGGTTGACTGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGt           <  1:2902945/150‑1 (MQ=255)

                                                                                                                                        cTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTa          <  2:4825029/150‑1 (MQ=255)

                                                                                                                                         ttGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTac         >  1:4964566/1‑150 (MQ=255)

                                                                                                                                         ttGAGGTTGACCGTGGTAAGGCGTTCATCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTAGGCTTGACGATGTac         >  1:3661867/1‑150 (MQ=255)

                                                                                                                                            aGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCCTCTTCGGCGTGACGATGTACAcg      >  2:2961374/1‑150 (MQ=255)

                                                                                                                                             ggTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTACAcgc     <  2:1146725/150‑1 (MQ=255)

                                                                                                                                              gTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc                                                 >  2:4683760/1‑105 (MQ=255)

                                                                                                                                              gTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCACGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTc                                           <  1:1539901/111‑1 (MQ=255)

                                                                                                                                                tGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTACAcgctgc  >  1:3254092/1‑150 (MQ=255)

TCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTCAGGGCCGTGCATCTTCGGCTTGACGATGTACACGCTGC  >  NC_002947/433744‑434037

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