BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,858,764	A→G	5.3%	D238G (GAC→GGC)	PP_5133 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,858,764	0	A	G	5.3%	1004.0 / 26.5	304	D238G (GAC→GGC)	PP_5133	membrane protein
Reads supporting (aligned to +/- strand): ref base A (163/124); new base G (8/8); total (172/132)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 313 total aligned reads displayed.

GCTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGCGCC  >  NC_002947/5858615‑5858913

 cTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgac                                                                                                                                                      >  1:3170480/1‑150 (MQ=255)

  tCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgacg                                                                                                                                                     >  1:3236604/1‑150 (MQ=255)

    aCCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGaa                                                                                                                                                   <  2:1970538/150‑1 (MQ=255)

     ccGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGTCGGGCACGGCGAAg                                                                                                                                                  <  1:2043645/150‑1 (MQ=255)

      cgcgCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGc                                                                                                                                                 >  2:1157303/1‑150 (MQ=255)

        cgcTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCa                                                                                                                                               <  1:2963996/150‑1 (MQ=255)

           tGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGa                                                                                                                                            >  2:3744490/1‑150 (MQ=255)

                aTGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                                       <  2:3380512/150‑1 (MQ=255)

                 tGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgcg                                                                                                                                      >  1:584088/1‑150 (MQ=255)

                    tACTGGTCATGAAGGAGACCCGGCTCTATGGCTTAGTCTGGTAAACCAACTTGCTCGCCTCCGACCCGTTCATACACCTGACCAAAGCACTGGGCGCCCTCCCGTCGTTGCTGTGATTCTCCGTGCCCGACGAAGCCATGATCAGCGCCa                                                                                                                                   >  2:4174804/1‑150 (MQ=255)

                        ggTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                               >  2:2209652/1‑150 (MQ=255)

                        ggTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                               <  1:2784471/150‑1 (MQ=255)

                         gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                              >  2:2055678/1‑150 (MQ=255)

                         gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                              >  1:3920350/1‑150 (MQ=255)

                         gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                              >  2:1750715/1‑150 (MQ=255)

                         gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCCGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                              >  1:3641053/1‑150 (MQ=255)

                          tCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggcg                                                                                                                             >  1:3049559/1‑150 (MQ=255)

                            ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCCGCGGCGAc                                                                                                                           >  2:914168/1‑150 (MQ=255)

                            ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGAc                                                                                                                           >  2:3937519/1‑150 (MQ=255)

                             tGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACa                                                                                                                          <  1:289887/150‑1 (MQ=255)

                                   cGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCAACAGCCTGc                                                                                                                    <  2:4416604/150‑1 (MQ=255)

                                   cGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCAACAGCCTGc                                                                                                                    <  2:3149155/150‑1 (MQ=255)

                                     aCCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCa                                                                                                                  <  1:2537434/150‑1 (MQ=255)

                                         ggCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGccc                                                                                                              >  2:3601100/1‑150 (MQ=255)

                                         ggCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGGCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGccc                                                                                                              <  2:2798827/150‑1 (MQ=255)

                                          gcgcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCt                                                                                                             <  2:1784525/150‑1 (MQ=255)

                                               aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACg                                                                                                        >  1:1494067/1‑150 (MQ=255)

                                                tGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCAGGACGt                                                                                                       >  1:4247527/1‑150 (MQ=255)

                                                  gCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTgg                                                                                                     <  2:2862247/150‑1 (MQ=255)

                                                    tttGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGcc                                                                                                   >  2:1446721/1‑150 (MQ=255)

                                                    tttGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGcc                                                                                                   >  2:330533/1‑150 (MQ=255)

                                                       gTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                                <  1:67760/150‑1 (MQ=255)

                                                         cTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGccctggc                                                                                              <  1:1641381/150‑5 (MQ=255)

                                                         cTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCa                                                                                              >  1:1849991/1‑150 (MQ=255)

                                                          tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAg                                                                                             >  2:4495832/1‑150 (MQ=255)

                                                          tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAg                                                                                             >  2:2262393/1‑150 (MQ=255)

                                                           gggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                            >  2:1431561/1‑150 (MQ=255)

                                                           gggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                            >  2:3645442/1‑150 (MQ=255)

                                                             gAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCCTGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGcc                                                                                          >  1:646483/1‑150 (MQ=255)

                                                              aaaCCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCt                                                                                         <  1:2635874/150‑1 (MQ=255)

                                                                  caccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGc                                                                                     >  2:1618307/1‑150 (MQ=255)

                                                                   accTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGcc                                                                                    >  2:1521723/1‑150 (MQ=255)

                                                                    ccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCa                                                                                   <  1:4191749/150‑1 (MQ=255)

                                                                       tGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgct                                                                                <  2:2570300/150‑1 (MQ=255)

                                                                       tGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgct                                                                                <  1:510751/150‑1 (MQ=255)

                                                                         cTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctgg                                                                              <  2:4014547/150‑1 (MQ=255)

                                                                         cTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctgg                                                                              <  2:3262870/150‑1 (MQ=255)

                                                                         cTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACGGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctgg                                                                              <  2:3499473/150‑1 (MQ=255)

                                                                         cTCGCCGCCGACCCGTTCAGTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctgg                                                                              <  2:2668990/150‑1 (MQ=255)

                                                                         cTCGCCGCCGACCCGTTCAGTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctgg                                                                              <  2:4239746/150‑1 (MQ=255)

                                                                          tCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggc                                                                             >  2:2837091/1‑150 (MQ=255)

                                                                           cgccgcCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCACCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggct                                                                            <  2:3126141/150‑1 (MQ=255)

                                                                            gccgccGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggctg                                                                           >  1:4518096/1‑150 (MQ=255)

                                                                                  gACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGc                                                                     >  1:174620/1‑150 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTTGGCTTAGCCGTGACCGACGAAGCCATGATCCGCGCCAGCGGCGACACCCTGCCACCCCAGGACGTGGCCCCCCAGGCCAGGGCCAGCTGGCTGCTCGGCgtg                                                                  >  1:410602/1‑150 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                                  <  2:3422196/150‑1 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                                  >  1:1415178/1‑150 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                                  <  2:2632738/150‑1 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                                  >  1:3500437/1‑150 (MQ=255)

                                                                                     ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                                  >  1:3825167/1‑150 (MQ=255)

                                                                                       gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtggt                                                                >  2:3886800/1‑150 (MQ=255)

                                                                                          cattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                             <  1:147453/150‑1 (MQ=255)

                                                                                          cattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                             <  1:1234196/150‑1 (MQ=255)

                                                                                           attcaCCTGACCCACGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTg                                                            <  2:1921550/150‑1 (MQ=255)

                                                                                           attcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTg                                                            >  1:1363489/1‑150 (MQ=255)

                                                                                            ttcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                                           >  1:2452802/1‑150 (MQ=255)

                                                                                             tcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                                          <  1:1069803/150‑1 (MQ=255)

                                                                                               aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGtat                                                        >  1:4393701/1‑150 (MQ=255)

                                                                                                  tGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCGGGTATACg                                                     >  1:2120492/1‑150 (MQ=255)

                                                                                                      ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCGGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCct                                                 >  1:2259676/1‑150 (MQ=255)

                                                                                                        aaGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctgc                                               >  2:3846104/1‑150 (MQ=255)

                                                                                                        aaGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctgc                                               >  1:215043/1‑150 (MQ=255)

                                                                                                        aaGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctgc                                               >  1:702166/1‑150 (MQ=255)

                                                                                                         aGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctgct                                              <  2:2121146/150‑1 (MQ=255)

                                                                                                              tGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcccc                                         <  2:3109303/150‑1 (MQ=255)

                                                                                                                ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGc                                       >  1:864124/1‑150 (MQ=255)

                                                                                                                 gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGcc                                      >  2:2831778/1‑150 (MQ=255)

                                                                                                                   gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctg                                    >  1:3699395/1‑150 (MQ=255)

                                                                                                                   gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCCACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctg                                    >  1:1309553/1‑150 (MQ=255)

                                                                                                                     ccTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCCGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgct                                  <  1:1786237/150‑1 (MQ=255)

                                                                                                                      cTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctg                                 >  2:2551486/1‑150 (MQ=255)

                                                                                                                      cTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctg                                 <  2:1394764/150‑1 (MQ=255)

                                                                                                                       tCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgc                                >  2:739888/1‑150 (MQ=255)

                                                                                                                             cgtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCg                          <  2:584088/150‑1 (MQ=255)

                                                                                                                              gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                         >  1:899680/1‑150 (MQ=255)

                                                                                                                              gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                         <  2:4369031/150‑1 (MQ=255)

                                                                                                                              gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                         >  1:2715766/1‑150 (MQ=255)

                                                                                                                                  cTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTa                     >  1:2199955/1‑150 (MQ=255)

                                                                                                                                   tgtgCTTCGCCGTGCCCGACGAAGCCATTATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGACATGCCTGTGCCAGCTGGCTGCTCGGCGTGTTGCTTGTATACTGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTAt                    <  1:3356539/147‑1 (MQ=255)

                                                                                                                                   tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCC‑GGTTATg                   >  2:2352796/1‑150 (MQ=255)

                                                                                                                                    gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                   >  1:3664106/1‑150 (MQ=255)

                                                                                                                                    gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                   <  2:1288100/150‑1 (MQ=255)

                                                                                                                                      gCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGcc                 >  1:3864224/1‑150 (MQ=255)

                                                                                                                                        ttCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtg               >  1:3736385/1‑150 (MQ=255)

                                                                                                                                          cGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgtg             <  1:2611761/150‑1 (MQ=255)

                                                                                                                                            ccGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggc           <  1:3133742/150‑1 (MQ=255)

                                                                                                                                              gTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgc         <  1:954284/150‑1 (MQ=255)

                                                                                                                                                  ccGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgctggc     <  1:1157303/150‑1 (MQ=255)

                                                                                                                                                     acgaAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGCGcc  >  2:519544/1‑150 (MQ=255)

                                                                                                                                                     acgaAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGCGcc  <  2:2986422/150‑1 (MQ=255)

GCTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGCGCC  >  NC_002947/5858615‑5858913

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