BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,394,501	T→A	6.5%	T411T (ACA→ACT)	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,394,501	0	T	A	6.5%	730.6 / 7.4	253	T411T (ACA→ACT)	PP_3865	tail protein
Reads supporting (aligned to +/- strand): ref base T (73/148); new base A (18/2); total (101/152)
Fisher's exact test for biased strand distribution p-value = 6.96e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

TAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAG  >  NC_002947/4394352‑4394650

tAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACt                                                                                                                                                       <  1:1761531/150‑1 (MQ=255)

    ccTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                                   <  1:1172097/150‑1 (MQ=255)

    ccTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                                   <  2:3826585/150‑1 (MQ=255)

    ccTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATAGGCCTGGCCGCGAAGGCCACACCGAGCACCCCCCTGGCCACTGTgg                                                                                                                                                   >  2:3431145/1‑150 (MQ=255)

     cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCACCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGc                                                                                                                                                  >  1:672371/1‑150 (MQ=255)

      ttcaTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCCAGCCCCCGCGTGGCCCCAGTGGcc                                                                                                                                                 >  1:5050052/4‑150 (MQ=255)

        gATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGTACCCGCGTGGCCACTGTGGCcgc                                                                                                                                               <  2:613780/150‑1 (MQ=255)

        gATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACAGTGGCcgc                                                                                                                                               >  2:2958503/1‑150 (MQ=255)

        gATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCCATCCCCCCGGCGATCAGGCCAAAAGGCGTGGCCCCGCAGGCCCCACCCAGCAACCCCGTGGCCACAGTGGCcgc                                                                                                                                               >  1:4461144/1‑150 (MQ=255)

          tCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCGCGGCGGTCAGGCCAATAGGCGTGGCCGCGAAGGCCCCACCGCGCAACCCCGTGGCCCCAGTgggcccgc                                                                                                                                             >  1:2251341/1‑144 (MQ=255)

          tCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgcgc                                                                                                                                             <  2:3550305/150‑2 (MQ=255)

           cAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATAGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGaa                                                                                                                                            <  2:478376/150‑1 (MQ=255)

            aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                           <  1:1123251/150‑4 (MQ=255)

            aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                           <  2:2627924/150‑4 (MQ=255)

             gCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCCCGGCGATCAGGCCCAAAGGCGTGGCCGCGAAGGCCACACCCCGCACCCCCGTGGCCACAGTGGCCGCGaccc                                                                                                                                          >  1:3884097/1‑147 (MQ=255)

               tCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCCCCGCCATCACGCCAAAAGGCGTGGCCGCGAAGGCCACACCGAGCCCCCCCCTGGCCACGGTGGCCGCGAACAcc                                                                                                                                        >  1:5408843/1‑150 (MQ=255)

               tCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGAACCCCGCGGCGATAAGGCCAAACGGCGTGGCCGCGAAGGCCACACCCAGCACCCCCGTGGCCACAGTGGCCCCGAACAcc                                                                                                                                        >  1:1351770/1‑150 (MQ=255)

                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccgc                                                                                                                                      <  1:3375191/150‑1 (MQ=255)

                  cacaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCGGGCGCGATCCCCCCGGCGGTCAGGCCAAAAGGCGTGGGGGCGAATGCCACAACGAGCACCCGCGTGGGCCCTGTGGGCGCGAAAAacgcc                                                                                                                                     >  1:1711691/1‑150 (MQ=255)

                  cacaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccgcc                                                                                                                                     <  1:3679133/150‑1 (MQ=255)

                   acaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCCCCGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCCAGCACCCGCGGGGCCCCGGTGGCCGCGAACAcccccc                                                                                                                                    >  1:1791301/1‑149 (MQ=255)

                       ggCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCGGCCGCGCGCCCCCCCGCGATCAGGCCAAAAGGCGTGGCCGCGCAGGCCCCCCCCAGCACCCCCGGGGCCCCAGTGGCCCCGCACACCCCcccacc                                                                                                                                >  1:5348581/1‑149 (MQ=255)

                         cgcTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCAATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCTCCCGCGTGGCCACTGTGGCCGCGAACACCGCcgcccgca                                                                                                                              <  2:2542883/150‑1 (MQ=255)

                           cTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAgc                                                                                                                            <  2:430282/150‑1 (MQ=255)

                              aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGCTCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCCCGGTGGCCGCGAACACCGCCGCACGCAGcgcg                                                                                                                         >  2:2700089/1‑150 (MQ=255)

                              aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACCGTGGCCGCGAACACCGCCGCACGCAGcgcg                                                                                                                         >  2:2685756/1‑150 (MQ=255)

                              aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCCCGGCGATAAGGCCCAACGGCGTGGCCGCGAAGGCCCCACCGAGCACCCCCGTGGCCACAGTGGCCGCGAACACCGCCCCACGCAGcgcg                                                                                                                         >  1:1782187/1‑150 (MQ=255)

                               agaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCgg                                                                                                                        <  2:2584145/150‑1 (MQ=255)

                                   cGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                    <  2:4057274/150‑1 (MQ=255)

                                   cGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                    <  1:667746/150‑1 (MQ=255)

                                    ggCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                                   <  1:2315205/150‑1 (MQ=255)

                                    ggCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                                   <  1:4719056/150‑1 (MQ=255)

                                         gACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAAACGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACGGTGGCCGCGAACACCCCCGCACGCAGCGCGGGGGAGgccgc                                                                                                              >  2:521096/1‑150 (MQ=255)

                                             ttCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGAACAGGCCAATAGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGGGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCCCGGTGAAGGCCGccccc                                                                                                          >  2:4365149/1‑149 (MQ=255)

                                             ttCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCCCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCCAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCg                                                                                                          >  2:2107121/1‑150 (MQ=255)

                                                  ccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAt                                                                                                     <  1:5070285/150‑1 (MQ=255)

                                                    aGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAtta                                                                                                   >  2:869919/1‑148 (MQ=255)

                                                    aGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGa                                                                                                   <  2:1139304/150‑1 (MQ=255)

                                                     gTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGAc                                                                                                  <  2:5406937/150‑1 (MQ=255)

                                                        ggCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGAt                                                                                               <  1:1087791/150‑1 (MQ=255)

                                                         gCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATg                                                                                              <  1:131854/150‑1 (MQ=255)

                                                            acgaTCAGACCGGCCGCCGCCGCGATCCCCCCGGCGATCAGGCCCATAGGCGTGGCCGCGAAGGCCACACCCAGCACCCCCGTGGCCACAGTGGCCGCGAACACCCCCGCACGCAGCGCGGTGAATGCCGCCCCCGCAAAGACGATGCcg                                                                                           >  1:3203780/1‑150 (MQ=255)

                                                                       ggCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTg                                                                                <  2:3840309/150‑1 (MQ=255)

                                                                        gccgccgccGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                                                               <  1:2794810/150‑1 (MQ=255)

                                                                         ccgccgccgCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCg                                                                              <  1:3769832/150‑1 (MQ=255)

                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  1:2188098/150‑1 (MQ=255)

                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  1:4984805/150‑1 (MQ=255)

                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  2:5032069/150‑1 (MQ=255)

                                                                           gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCggg                                                                            <  2:4224436/150‑1 (MQ=255)

                                                                             cgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                          >  1:2999019/1‑150 (MQ=255)

                                                                                cgcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTc                                                                       <  1:5580161/150‑1 (MQ=255)

                                                                                cgcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTc                                                                       <  2:4928852/150‑1 (MQ=255)

                                                                                     tCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                                  >  2:2079945/1‑150 (MQ=255)

                                                                                         cgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatcat                                                              <  1:2290343/150‑1 (MQ=255)

                                                                                                    ggCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                                   >  1:4901264/1‑150 (MQ=255)

                                                                                                     gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATg                                                  >  1:2018120/1‑150 (MQ=255)

                                                                                                       cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                                <  2:5166892/150‑1 (MQ=255)

                                                                                                        aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCt                                               <  2:3973663/150‑1 (MQ=255)

                                                                                                        aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCt                                               <  2:1782187/150‑1 (MQ=255)

                                                                                                             gCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATcc                                          >  1:4561327/1‑150 (MQ=255)

                                                                                                               gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                        >  2:4857518/1‑150 (MQ=255)

                                                                                                               gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                        <  2:3050512/150‑1 (MQ=255)

                                                                                                               gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                        <  2:3249512/150‑1 (MQ=255)

                                                                                                               gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                        <  1:3279309/150‑1 (MQ=255)

                                                                                                                tGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                       >  2:3601709/1‑150 (MQ=255)

                                                                                                                tGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                       <  2:3734255/150‑1 (MQ=255)

                                                                                                                tGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                       >  2:779432/1‑150 (MQ=255)

                                                                                                                 ggCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACgg                                      <  2:1608502/150‑1 (MQ=255)

                                                                                                                   ccGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGcc                                    >  1:134528/1‑150 (MQ=255)

                                                                                                                       gAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAg                                <  1:1013520/150‑1 (MQ=255)

                                                                                                                         aGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGtt                              <  1:4259722/150‑1 (MQ=255)

                                                                                                                          ggCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTg                             >  1:5521153/1‑150 (MQ=255)

                                                                                                                            ccACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCt                           >  2:1394179/1‑150 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAgg                                                                                                                  >  2:5311334/1‑62 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCt                           <  2:4318749/149‑1 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCt                           >  1:4318749/1‑149 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  1:2593660/150‑1 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  2:5443650/150‑1 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  2:5513831/150‑1 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  1:2647110/150‑1 (MQ=255)

                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          >  2:2384678/1‑150 (MQ=255)

                                                                                                                               cacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                        >  2:1228388/1‑150 (MQ=255)

                                                                                                                               cacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                        <  1:1559046/150‑1 (MQ=255)

                                                                                                                                   gAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCAc                    <  1:5453174/150‑1 (MQ=255)

                                                                                                                                       aCCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCa                >  1:2082184/1‑150 (MQ=255)

                                                                                                                                       aCCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCa                >  1:867482/1‑150 (MQ=255)

                                                                                                                                       aCCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCa                <  2:4561327/150‑1 (MQ=255)

                                                                                                                                         ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAtt              <  2:2251341/150‑1 (MQ=255)

                                                                                                                                           gcgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGa            <  2:2005766/150‑1 (MQ=255)

                                                                                                                                            cgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGAt           <  1:3604750/150‑1 (MQ=255)

                                                                                                                                               ggCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGcc        <  2:134528/150‑1 (MQ=255)

                                                                                                                                                gCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg       >  2:3881787/1‑150 (MQ=255)

                                                                                                                                                gCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg       <  2:368500/150‑1 (MQ=255)

                                                                                                                                                 ccACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCgg      <  1:2016577/150‑1 (MQ=255)

                                                                                                                                                  cacaGTGGCCACGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCa                                                                                    <  1:5173967/71‑1 (MQ=255)

                                                                                                                                                  cACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCggg     >  1:793498/1‑150 (MQ=255)

                                                                                                                                                   aCTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgc    <  2:482158/150‑1 (MQ=255)

                                                                                                                                                    cTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgca   <  2:346897/150‑1 (MQ=255)

                                                                                                                                                    cTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgca   <  1:4942573/150‑1 (MQ=255)

                                                                                                                                                     tgtgGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgcag  >  2:5100006/1‑150 (MQ=255)

TAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAG  >  NC_002947/4394352‑4394650

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