BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,428,505	T→G	5.4%	D102A (GAC→GCC)	sohB ←	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,428,505	0	T	G	5.4%	309.4 / 5.7	109	D102A (GAC→GCC)	sohB	peptidase
Reads supporting (aligned to +/- strand): ref base T (43/60); new base G (3/3); total (46/63)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 112 total aligned reads displayed.

ACGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTCCTTG  >  NC_002947/4428357‑4428654

aCGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTc                                                                                                                                                      <  2:96604/150‑1 (MQ=255)

   cgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGt                                                                                                                                                       >  2:1519535/1‑146 (MQ=255)

   cgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGt                                                                                                                                                       <  1:1519535/146‑1 (MQ=255)

   cgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGcc                                                                                                                                                   >  1:2199096/1‑150 (MQ=255)

       aGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCa                                                                                                                                               >  2:2106585/1‑150 (MQ=255)

        gTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCaa                                                                                                                                              >  2:854333/1‑150 (MQ=255)

         ttGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCaaa                                                                                                                                             <  2:946450/150‑1 (MQ=255)

            tctcGCACCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                          <  2:703527/147‑1 (MQ=255)

                 cAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                                     >  1:1887342/1‑150 (MQ=255)

                 cAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTAACACGGCGTCGCGCGCGGTCGCGAGCTTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                                     <  2:1172491/150‑1 (MQ=255)

                   gCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACg                                                                                                                                   >  2:14777/1‑150 (MQ=255)

                    ccAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGa                                                                                                                                  <  2:1205021/150‑1 (MQ=255)

                     cAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGGTTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGaa                                                                                                                                 >  2:1929899/1‑150 (MQ=255)

                     cAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGaa                                                                                                                                 >  2:270339/1‑150 (MQ=255)

                        gCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCa                                                                                                                                                  <  2:2133636/130‑1 (MQ=255)

                        gCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCa                                                                                                                                                  >  1:2133636/1‑130 (MQ=255)

                         ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAcc                                                                                                                             >  1:328302/1‑150 (MQ=255)

                         ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAcc                                                                                                                             <  1:14777/150‑1 (MQ=255)

                          cGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAccc                                                                                                                            >  2:626632/1‑150 (MQ=255)

                              aCTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGCCATCAAAGTCCAGTACGAATACCCGGct                                                                                                                        <  2:1487191/150‑1 (MQ=255)

                                 gtgCACCAGGCCACCGGCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCtt                                                                                                                     <  2:120742/150‑1 (MQ=255)

                                 gtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCtt                                                                                                                     <  1:867582/150‑1 (MQ=255)

                                  tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTg                                                                                                                    <  1:2106585/150‑1 (MQ=255)

                                  tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTg                                                                                                                    <  1:2157278/150‑1 (MQ=255)

                                     aCCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGACATCaaa                                                                                                                                             >  1:993357/1‑122 (MQ=255)

                                       cAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTc                                                                                                               >  2:2375770/1‑150 (MQ=255)

                                          gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCAt                                                                                                                                                 >  1:1031892/1‑113 (MQ=255)

                                          gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCAt                                                                                                                                                 <  2:1031892/113‑1 (MQ=255)

                                          gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGc                                                                                                            <  2:1981978/150‑1 (MQ=255)

                                           ccACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGcc                                                                                                           >  1:800739/1‑150 (MQ=255)

                                            cACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCt                                                                                                          >  2:1664131/1‑150 (MQ=255)

                                             aCCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAg                                                                                                                                       >  1:2189346/1‑120 (MQ=255)

                                             aCCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAg                                                                                                                                       <  2:2189346/120‑1 (MQ=255)

                                             aCCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                       <  2:2356726/136‑1 (MQ=255)

                                             aCCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCtt                                                                                                         >  2:2249015/1‑150 (MQ=255)

                                              ccgccgCTTTCCAGGCGCAGTACCCATTCGTCACTCGCTGTAGCTGTTGTCAGCAGTGCGGTGATTTCGTCGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTg                                                                                                        <  1:539339/150‑1 (MQ=255)

                                               cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAccc                                                                                                                            >  1:1111794/1‑129 (MQ=255)

                                               cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAccc                                                                                                                            <  2:1111794/129‑1 (MQ=255)

                                                 ccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATCCCCGGCTCTTGTCCTCGGCCTTGccc                                                                                                     <  1:1032625/150‑1 (MQ=255)

                                                  cgcTTTCCAGGCGCAGCACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCct                                                                                                    >  1:2271201/1‑150 (MQ=255)

                                                    cTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCcttc                                                                                                  <  2:1401222/150‑1 (MQ=255)

                                                        ccAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGACGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTtcttct                                                                                              <  1:1374545/150‑1 (MQ=255)

                                                         cAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctg                                                                                             <  1:1818168/150‑1 (MQ=255)

                                                          aGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCATTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgc                                                                                            >  1:140230/1‑150 (MQ=255)

                                                           ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgct                                                                                           >  2:1470857/1‑150 (MQ=255)

                                                           ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgct                                                                                           <  1:626632/150‑1 (MQ=255)

                                                            gcgcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                                     >  2:559847/1‑107 (MQ=255)

                                                            gcgcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                                     <  1:559847/107‑1 (MQ=255)

                                                             cgcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGt                                                                                         <  2:2199096/150‑1 (MQ=255)

                                                              gcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGtt                                                                                        >  1:678851/1‑150 (MQ=255)

                                                                   aCCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTTCCACGCGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTcc                                                                                   <  1:760911/150‑1 (MQ=255)

                                                                   aCCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTcc                                                                                   <  1:1664131/150‑1 (MQ=255)

                                                                   aCCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTcc                                                                                   >  1:678283/1‑150 (MQ=255)

                                                                   aCCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCCGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCGCGGCCTTGCCCTGCTTCTGCTGTTTTTcc                                                                                   <  2:526306/150‑1 (MQ=255)

                                                                         tCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTGCTGCTGTTTTTCCGCCTTg                                                                             <  1:53924/150‑1 (MQ=255)

                                                                                gcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGt                                                                                         >  1:900534/1‑131 (MQ=255)

                                                                                gcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGt                                                                                         <  2:900534/131‑1 (MQ=255)

                                                                                gcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctgct                                                                      >  1:1290712/1‑150 (MQ=255)

                                                                                  gcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCtgctgc                                                                    >  1:580036/1‑150 (MQ=255)

                                                                                  gcgGTAGCGAGCGTCAGCAGTGCGGGGATTTCGTTGCGCAGGCTTTGGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCtgctgc                                                                    >  2:1090874/1‑150 (MQ=255)

                                                                                               tCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCt                                                       <  1:363031/150‑1 (MQ=255)

                                                                                                  gcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGc                                                                           <  1:1069656/127‑1 (MQ=255)

                                                                                                  gcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGc                                                                           >  2:1069656/1‑127 (MQ=255)

                                                                                                    agTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGc                                                  <  2:678851/150‑1 (MQ=255)

                                                                                                      tGCGTTTATTTTGTTGGGCAGGTCCGCGTTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTg                                                <  2:140230/150‑1 (MQ=255)

                                                                                                       gCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCTCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  1:270339/150‑1 (MQ=255)

                                                                                                       gCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  1:854333/150‑1 (MQ=255)

                                                                                                       gCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  2:2193349/150‑1 (MQ=255)

                                                                                                       gCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  1:2159385/150‑1 (MQ=255)

                                                                                                       gCGGTGATTTCGGTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  2:433527/150‑1 (MQ=255)

                                                                                                         ggTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTTTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGc                                             <  1:773810/150‑1 (MQ=255)

                                                                                                         ggTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGc                                             <  1:432179/150‑1 (MQ=255)

                                                                                                         ggTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGc                                             <  1:2375770/150‑1 (MQ=255)

                                                                                                          gTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCTTGCTGCTTGCGCAGGGCCTTGAGCTGGGcc                                            >  2:1985788/1‑150 (MQ=255)

                                                                                                            gATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCtt                                          <  2:580036/150‑1 (MQ=255)

                                                                                                             aTTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTg                                         >  1:743476/1‑150 (MQ=255)

                                                                                                               ttCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTc                                       >  1:2164156/1‑150 (MQ=255)

                                                                                                                 cGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGACTTGTAGa                                     <  2:1802463/150‑1 (MQ=255)

                                                                                                                  gTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               <  1:1441981/139‑1 (MQ=255)

                                                                                                                  gTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                               >  2:1441981/1‑139 (MQ=255)

                                                                                                                       gcAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGc                               >  2:1992750/1‑150 (MQ=255)

                                                                                                                        cAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGcc                              <  2:1887342/150‑1 (MQ=255)

                                                                                                                         aGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCTGCTCGTGGCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCg                             <  2:328302/150‑1 (MQ=255)

                                                                                                                         aGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCg                             <  2:1964895/150‑1 (MQ=255)

                                                                                                                          ggCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                                      >  1:423799/1‑124 (MQ=255)

                                                                                                                          ggCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                                      <  2:423799/124‑1 (MQ=255)

                                                                                                                           gCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGa                           <  2:800739/150‑1 (MQ=255)

                                                                                                                           gCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCACTTCTTCTGCTTTTTTTCCGCCTTGGCCAGCTGCTTGCGCAGGCCCTTGAGCTGGGGCTTGTCGAGCCGGCCGGa                           >  2:422582/1‑150 (MQ=255)

                                                                                                                                 ggTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAg                     <  1:1850358/150‑1 (MQ=255)

                                                                                                                                      ccGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTCGc                                                                           <  1:2407200/91‑1 (MQ=255)

                                                                                                                                       cGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGACTACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTc               <  2:1290712/150‑1 (MQ=255)

                                                                                                                                       cGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTc               >  2:1969914/1‑150 (MQ=255)

                                                                                                                                         aGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCAc             <  2:2271201/150‑1 (MQ=255)

                                                                                                                                          gacTTTGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg            >  1:56944/3‑150 (MQ=255)

                                                                                                                                             tttGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAg         <  1:1470857/150‑1 (MQ=255)

                                                                                                                                             tttGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAACGCACCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAg         >  1:24043/1‑150 (MQ=255)

                                                                                                                                             tttGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCtt                                                                              <  2:1152796/81‑1 (MQ=255)

                                                                                                                                             tttGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCtt                                                                              >  1:1152796/1‑81 (MQ=255)

                                                                                                                                                  tGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTCCt    >  1:184834/1‑150 (MQ=255)

                                                                                                                                                    tCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGGCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTTAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTCCTTg  <  1:281527/150‑1 (MQ=255)

ACGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTCCTTG  >  NC_002947/4428357‑4428654

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