BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,990,968	A→T	6.0%	D507V (GAC→GTC)	oplA →	5‑oxoprolinase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,990,968	0	A	T	6.0%	675.6 / 19.1	198	D507V (GAC→GTC)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand): ref base A (98/88); new base T (2/10); total (100/98)
Fisher's exact test for biased strand distribution p-value = 1.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

TGCGGCGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATG  >  NC_002947/3990819‑3991117

tGCGGCGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGa                                                                                                                                                       <  1:3349807/150‑1 (MQ=255)

    gcggGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGcc                                                                                                                                                   >  2:913267/1‑150 (MQ=255)

      gggCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCgg                                                                                                                                                 <  1:2833057/150‑1 (MQ=255)

          aTGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCAc                                                                                                                                             >  1:1068960/1‑150 (MQ=255)

          aTGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCAc                                                                                                                                             <  1:1299711/150‑1 (MQ=255)

           tGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGTTGCGCCAGGTCGACGTGCCCg                                                                                                                                            <  1:3018926/150‑1 (MQ=255)

           tGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACg                                                                                                                                            <  2:2365379/150‑1 (MQ=255)

            gCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGAGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGt                                                                                                                                           >  1:1782785/1‑150 (MQ=255)

             cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTc                                                                                                                                          <  2:4226373/150‑1 (MQ=255)

             cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTc                                                                                                                                          <  1:1766414/150‑1 (MQ=255)

             cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTc                                                                                                                                          <  2:2254193/150‑1 (MQ=255)

                  ggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTcca                                                                                                                                     <  2:3964157/150‑2 (MQ=255)

                   gtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTccag                                                                                                                                    <  2:1296523/150‑3 (MQ=255)

                       gcCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGgcgc                                                                                                                                >  2:806610/1‑150 (MQ=255)

                        ccTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCa                                                                                                                               >  2:3161815/1‑150 (MQ=255)

                         cTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                              <  1:545477/150‑2 (MQ=255)

                          tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                             <  1:1201070/150‑1 (MQ=255)

                          tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                             >  1:3439402/1‑150 (MQ=255)

                          tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                             <  2:2468368/150‑1 (MQ=255)

                          tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                             <  2:2506852/150‑1 (MQ=255)

                          tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCatc                                                                                                                             <  2:2497912/150‑3 (MQ=38)

                             ccGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCa                                                                                                                          <  2:1673992/150‑1 (MQ=255)

                               gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAAc                                                                                                                        >  1:4467217/1‑150 (MQ=255)

                               gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAAc                                                                                                                        >  1:1019034/1‑150 (MQ=255)

                                   aaaTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAgg                                                                                                                    <  1:2048632/150‑1 (MQ=255)

                                    aaTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGc                                                                                                                   <  2:1444328/150‑1 (MQ=255)

                                       tGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCt                                                                                                                <  1:651505/150‑1 (MQ=255)

                                        ggCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTc                                                                                                               <  2:4039707/150‑1 (MQ=255)

                                          cATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGc                                                                                                             >  2:4293317/1‑150 (MQ=255)

                                           aTGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGcc                                                                                                            >  2:2415892/1‑150 (MQ=255)

                                              cacaCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                         >  2:2289112/1‑150 (MQ=255)

                                              cacaCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                         >  1:2334300/1‑150 (MQ=255)

                                              cacaCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                         >  1:81394/1‑150 (MQ=255)

                                               acacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGa                                                                                                                                                       >  1:3079614/1‑103 (MQ=255)

                                               acacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                        <  1:1812894/150‑1 (MQ=255)

                                               acacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                        <  1:3709098/150‑1 (MQ=255)

                                                 acCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTa                                                                                                      >  1:639684/1‑150 (MQ=255)

                                                  ccGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAc                                                                                                     <  2:1101773/150‑1 (MQ=255)

                                                     tGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTc                                                                                                  >  2:1723730/1‑150 (MQ=255)

                                                        tGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGc                                                                                               >  2:3811605/1‑150 (MQ=255)

                                                              ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGcc                                                                                         >  1:3531957/1‑150 (MQ=255)

                                                                     ggCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCaa                                                                                  <  1:1701036/150‑1 (MQ=255)

                                                                      gCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAAc                                                                                 >  2:4480623/1‑150 (MQ=255)

                                                                      gCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAAc                                                                                 <  2:4161083/150‑1 (MQ=255)

                                                                        cTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCt                                                                               >  2:2690159/1‑150 (MQ=255)

                                                                         tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTg                                                                              >  1:3101829/1‑150 (MQ=255)

                                                                          cGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTgc                                                                             <  1:3007010/150‑1 (MQ=255)

                                                                          cGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTgc                                                                             <  2:4461241/150‑1 (MQ=255)

                                                                            ggCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTgcgc                                                                           <  1:4094645/150‑1 (MQ=255)

                                                                               ttGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                        >  1:1084534/1‑150 (MQ=255)

                                                                                 gtgtgCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGa                                                                      <  1:2528363/150‑1 (MQ=255)

                                                                                  tgtgCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGAc                                                                     <  2:2531852/150‑1 (MQ=255)

                                                                                  tgtgCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGAc                                                                     <  2:3263565/150‑1 (MQ=255)

                                                                                     gCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGc                                                                  >  2:470690/1‑150 (MQ=255)

                                                                                         tCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTg                                                              <  2:1510294/150‑1 (MQ=255)

                                                                                          cGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTgg                                                             <  2:2900922/150‑1 (MQ=255)

                                                                                           ggCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGc                                                            >  2:3920066/1‑150 (MQ=255)

                                                                                              cAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCCACCTGCGCGAAGACGGCTTTGGCGAc                                                         >  1:9925/1‑150 (MQ=255)

                                                                                               aGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGAcc                                                        >  1:196918/1‑150 (MQ=255)

                                                                                               aGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGAcc                                                        >  1:2504056/1‑150 (MQ=255)

                                                                                               aGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGAcc                                                        <  1:722413/150‑1 (MQ=255)

                                                                                                gATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCa                                                       <  1:2607815/150‑1 (MQ=255)

                                                                                                 aTCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                      <  1:3143020/150‑1 (MQ=255)

                                                                                                   ccTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGc                                                    >  2:702848/1‑150 (MQ=255)

                                                                                                    cTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGcc                                                   >  1:1423902/1‑150 (MQ=255)

                                                                                                       gccGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTCGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCcagc                                                >  2:2756166/2‑150 (MQ=255)

                                                                                                           aCGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGtt                                            <  1:1112290/150‑1 (MQ=255)

                                                                                                               gcgcTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCa                                        <  2:1068960/150‑1 (MQ=255)

                                                                                                                 gcTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACt                                      >  1:2090217/1‑150 (MQ=255)

                                                                                                                  cTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                     >  2:2205394/1‑150 (MQ=255)

                                                                                                                  cTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCGGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCGTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                     >  1:622831/1‑150 (MQ=255)

                                                                                                                   tACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTac                                    >  2:222860/1‑150 (MQ=255)

                                                                                                                   tACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTac                                    >  2:1580769/1‑150 (MQ=255)

                                                                                                                      gACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCCCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACAcc                                 >  1:3796976/1‑150 (MQ=255)

                                                                                                                      gACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACAcc                                 >  2:311876/1‑150 (MQ=255)

                                                                                                                         cAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                              <  1:61154/150‑1 (MQ=255)

                                                                                                                          aGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCGACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGGTGCCCCTACACCCTgg                             >  2:527414/1‑150 (MQ=255)

                                                                                                                          aGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTgg                             >  2:3506771/1‑150 (MQ=255)

                                                                                                                           gCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                            <  2:4351658/150‑1 (MQ=255)

                                                                                                                             tCACCAGCCTGTCGATGCGCCTGGACGCCGGGCCCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACGACGTCAGCCAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGgaga                          >  1:2259599/1‑150 (MQ=255)

                                                                                                                               aCCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATg                        >  2:623135/1‑150 (MQ=255)

                                                                                                                               aCCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATg                        >  1:1585059/1‑150 (MQ=255)

                                                                                                                                ccAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgc                       <  1:2557091/150‑1 (MQ=255)

                                                                                                                                ccAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgc                       <  1:547163/150‑1 (MQ=255)

                                                                                                                                ccAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgc                       <  1:278081/150‑1 (MQ=255)

                                                                                                                                 cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcg                      >  2:3186482/1‑150 (MQ=255)

                                                                                                                                 cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcg                      <  2:3065634/150‑1 (MQ=255)

                                                                                                                                  aGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcgc                     <  2:3754138/150‑1 (MQ=255)

                                                                                                                                  aGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcgc                     <  2:3357231/150‑1 (MQ=255)

                                                                                                                                   gCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCt                    >  1:4222620/1‑150 (MQ=255)

                                                                                                                                       gTCGATGCGCCTGGTCGCCGGGCACGTCGACCTGGCGCAACTCATCCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCt                >  2:1389952/1‑150 (MQ=255)

                                                                                                                                            tGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGc           >  1:316696/1‑150 (MQ=255)

                                                                                                                                             gcgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGtt                                            >  2:1479194/1‑116 (MQ=255)

                                                                                                                                             gcgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGtt                                            <  1:1479194/116‑1 (MQ=255)

                                                                                                                                             gcgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCa          >  2:4077914/1‑150 (MQ=255)

                                                                                                                                             gcgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCa          >  1:3568185/1‑150 (MQ=255)

                                                                                                                                             gcgcCAGGTCGACGTGCCCGGCGTCCAGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCa                                        >  1:2233129/1‑120 (MQ=38)

                                                                                                                                                ccTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGAt       >  2:3614353/1‑150 (MQ=255)

                                                                                                                                                   ggACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCa    >  1:2623387/1‑150 (MQ=255)

                                                                                                                                                    gACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAt   <  2:1939227/150‑1 (MQ=255)

TGCGGCGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATG  >  NC_002947/3990819‑3991117

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: `‑`