BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,654,977	T→G	5.7%	H311P (CAC→CCC)	PP_4113 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,654,977	0	T	G	5.7%	783.6 / 24.6	228	H311P (CAC→CCC)	PP_4113	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (97/117); new base G (5/8); total (103/125)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 236 total aligned reads displayed.

TTGACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4654830‑4655126

ttGACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGc                                                                                                                                                     >  1:4426055/1‑150 (MQ=255)

  gACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCtt                                                                                                                                                   >  1:4588807/1‑150 (MQ=255)

     aTCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                >  1:3801586/1‑150 (MQ=255)

      tCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGccc                                                                                                                                               <  1:789904/150‑1 (MQ=255)

        gcgcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAg                                                                                                                                             >  1:2020521/1‑150 (MQ=255)

          gcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGcc                                                                                                                                           >  1:1966457/1‑150 (MQ=255)

           cgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGGGCAAGCACAGccc                                                                                                                                          <  2:2542923/150‑1 (MQ=255)

              gcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTgg                                                                                                                                       <  2:269587/150‑1 (MQ=255)

               ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGc                                                                                                                                      <  1:1271356/150‑1 (MQ=255)

                    ccGGTGAATGCCCCAGCTCCAGGTCAGCACCCGCACGCCATCTTCGACCAGGTTCCCCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCaaaac                                                                                                                                 >  1:2249585/1‑149 (MQ=255)

                     cGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGGGCAAGCACAGCCCTGGcttcgag                                                                                                                                <  2:3924669/150‑7 (MQ=255)

                       gTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                              <  1:1010745/150‑1 (MQ=255)

                       gTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                              <  1:2805788/150‑1 (MQ=255)

                        tGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                             <  1:4976115/150‑1 (MQ=255)

                        tGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                             <  1:4658071/150‑1 (MQ=255)

                         gAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCt                                                                                                                            >  1:2986875/1‑150 (MQ=255)

                            tGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCAc                                                                                                                         >  1:5580454/1‑150 (MQ=255)

                                ccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTg                                                                                                                     >  1:1781667/1‑150 (MQ=255)

                                ccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTg                                                                                                                     >  2:3733440/1‑150 (MQ=255)

                                     tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAg                                                                                                                >  1:1912204/1‑150 (MQ=255)

                                     tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAg                                                                                                                >  2:5298208/1‑150 (MQ=255)

                                      ccAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                               <  1:4129090/150‑1 (MQ=255)

                                      ccAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                               <  2:1589369/150‑1 (MQ=255)

                                       cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCg                                                                                                              >  1:5230353/1‑150 (MQ=255)

                                        aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCga                                                                                                             <  1:907623/150‑1 (MQ=255)

                                         gTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCgag                                                                                                            >  1:113862/1‑150 (MQ=255)

                                            cAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAAt                                                                                                         <  2:1545025/150‑1 (MQ=255)

                                              gCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACAGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGATGAGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCATGCGCAGCGAGAATGc                                                                                                       >  2:5257730/1‑150 (MQ=255)

                                              gCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGGTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGc                                                                                                       >  2:3209043/1‑150 (MQ=255)

                                                aCCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCt                                                                                                     <  1:3690169/150‑1 (MQ=255)

                                                aCCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCt                                                                                                     <  1:2641886/150‑1 (MQ=255)

                                                   cGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCg                                                                                                  >  1:917812/1‑150 (MQ=255)

                                                    gCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGa                                                                                                 <  1:5477646/150‑1 (MQ=255)

                                                    gCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGa                                                                                                 >  1:5140550/1‑150 (MQ=255)

                                                       cGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                              <  2:2573078/150‑1 (MQ=255)

                                                        gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATg                                                                                             <  1:995625/150‑1 (MQ=255)

                                                        gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATg                                                                                             <  2:2378327/150‑1 (MQ=255)

                                                          cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCaaa                                                                                                                                   <  2:178528/110‑1 (MQ=255)

                                                          cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACCTGGc                                                                                           >  2:288357/1‑150 (MQ=255)

                                                          cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGc                                                                                           >  1:5596666/1‑150 (MQ=255)

                                                            tCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTg                                                                                         <  2:3346641/150‑1 (MQ=255)

                                                                cGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGt                                                                                     <  2:2282788/150‑1 (MQ=255)

                                                                  aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                   <  1:4745855/150‑1 (MQ=255)

                                                                  aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                   >  1:2947485/1‑150 (MQ=255)

                                                                  aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                   >  2:1071322/1‑150 (MQ=255)

                                                                       gTTCACCGAGGCAGCGCTGGGGGCGGTGATACCGGCATCCGAGTTGCGCTCTCCCAGCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGACTGCCTGCGACATGGCTGCCGTGGCTTGc                                                                              >  2:2386445/1‑150 (MQ=255)

                                                                       gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGc                                                                              <  2:5516715/150‑1 (MQ=255)

                                                                         tCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCg                                                                            <  1:2688845/150‑1 (MQ=255)

                                                                            ccGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGc                                                                         <  1:4491773/150‑1 (MQ=255)

                                                                            ccGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGc                                                                         >  1:3611566/1‑150 (MQ=255)

                                                                             cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCGGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  2:4322644/150‑1 (MQ=255)

                                                                             cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  1:3367906/150‑1 (MQ=255)

                                                                             cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  1:622436/150‑1 (MQ=255)

                                                                             cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  1:330471/150‑1 (MQ=255)

                                                                                 gCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                    >  2:2364413/1‑150 (MQ=255)

                                                                                 gCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                    <  2:3982495/150‑1 (MQ=255)

                                                                                    gCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCg                                                                 >  2:4404830/1‑150 (MQ=255)

                                                                                    gCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCg                                                                 >  2:757080/1‑150 (MQ=255)

                                                                                     cGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGt                                                                <  2:4822637/150‑1 (MQ=255)

                                                                                      gATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTc                                                               <  1:4840680/150‑1 (MQ=255)

                                                                                          ttGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGa                                                           <  2:1481860/150‑1 (MQ=255)

                                                                                           tGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAg                                                          >  2:1253469/1‑150 (MQ=255)

                                                                                            ggCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCAAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGAGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                         <  2:4158940/150‑1 (MQ=255)

                                                                                            ggCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCAAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGAGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                         <  2:1305555/150‑1 (MQ=255)

                                                                                             gCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCg                                                        <  1:3768620/150‑1 (MQ=255)

                                                                                                 tGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTaga                                                    <  1:171145/150‑1 (MQ=255)

                                                                                                  gATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTagag                                                   <  2:3163416/150‑1 (MQ=255)

                                                                                                    tACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                 <  1:1555036/150‑1 (MQ=255)

                                                                                                      ccGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                               <  1:1241657/150‑1 (MQ=255)

                                                                                                       ctggcttcGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                              >  2:4841016/8‑150 (MQ=255)

                                                                                                       cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCAGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                              >  2:5586108/1‑150 (MQ=255)

                                                                                                           aTCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                          <  1:3048137/150‑1 (MQ=255)

                                                                                                             ccGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                        <  1:4336044/150‑1 (MQ=255)

                                                                                                             ccGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                        <  1:1630934/150‑1 (MQ=255)

                                                                                                              cGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGcc                                       >  2:1684155/1‑150 (MQ=255)

                                                                                                               gAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCt                                      >  1:1473462/1‑150 (MQ=255)

                                                                                                                  ttGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGt                                   >  2:2058367/1‑150 (MQ=255)

                                                                                                                       ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                              <  2:306721/150‑1 (MQ=255)

                                                                                                                       ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                              <  2:5430820/150‑1 (MQ=255)

                                                                                                                       ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                              <  2:4693547/150‑1 (MQ=255)

                                                                                                                       ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                              >  2:1594984/1‑150 (MQ=255)

                                                                                                                           aCCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                          >  1:1737351/1‑150 (MQ=255)

                                                                                                                            ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcccc                         >  2:1446761/1‑150 (MQ=255)

                                                                                                                            ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcccc                         <  1:2308909/150‑1 (MQ=255)

                                                                                                                               tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                      <  1:3641934/150‑1 (MQ=255)

                                                                                                                                    tCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                 >  2:4069057/1‑150 (MQ=255)

                                                                                                                                    tCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                 >  2:2750171/1‑150 (MQ=255)

                                                                                                                                      tttgCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg               >  1:247476/4‑150 (MQ=255)

                                                                                                                                      gAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg               <  1:3332795/150‑1 (MQ=255)

                                                                                                                                      gAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg               <  1:5149500/150‑1 (MQ=255)

                                                                                                                                        aGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGc             <  1:1172326/150‑1 (MQ=255)

                                                                                                                                         gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt            >  2:4182696/1‑150 (MQ=255)

                                                                                                                                         gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt            <  2:909711/150‑1 (MQ=255)

                                                                                                                                           cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAcc                                                                  <  1:599804/92‑1 (MQ=38)

                                                                                                                                           cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAcc                                                                  >  2:599804/3‑94 (MQ=38)

                                                                                                                                            ggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGaa         >  2:1887643/2‑150 (MQ=255)

                                                                                                                                            aGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGaa         >  2:3526706/1‑150 (MQ=255)

                                                                                                                                                 tgtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCg    >  1:3136219/1‑150 (MQ=255)

                                                                                                                                                   tgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                 >  1:1639997/1‑103 (MQ=255)

                                                                                                                                                   tgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                 <  2:1639997/103‑1 (MQ=255)

TTGACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4654830‑4655126

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