BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,836,771	T→C	9.0%	L220P (CTG→CCG)	yfiP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,836,771	0	T	C	9.0%	829.5 / 49.3	277	L220P (CTG→CCG)	yfiP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (105/147); new base C (8/17); total (113/164)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 280 total aligned reads displayed.

CCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGT  >  NC_002947/1836622‑1836920

ccTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCt                                                                                                                                                       <  2:52179/150‑1 (MQ=255)

  tGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGc                                                                                                                                                     <  1:6312611/150‑1 (MQ=255)

  tGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGc                                                                                                                                                     >  2:670299/1‑150 (MQ=255)

       cgcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGTCAAGCGGCAGCCGCCGCtt                                                                                                                                                <  1:6083019/150‑1 (MQ=255)

       cgcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                                                <  2:26465/150‑1 (MQ=255)

       cgcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                                                <  2:2786614/150‑1 (MQ=255)

       cgcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCtt                                                                                                                                                <  2:5741462/150‑1 (MQ=255)

        gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                                               <  1:1766864/150‑1 (MQ=255)

        gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                                               <  2:2224632/150‑1 (MQ=255)

        gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                                               <  2:3050808/150‑1 (MQ=255)

          tCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                             <  1:5273205/150‑1 (MQ=255)

             aaGCGCGCTGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATGACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                          <  2:631197/150‑1 (MQ=255)

             aaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                          <  2:2875295/150‑1 (MQ=255)

             aaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                          <  2:6528949/150‑1 (MQ=255)

             aaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                          >  1:2452218/1‑150 (MQ=255)

               gcgcGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                                        <  1:3539174/150‑1 (MQ=255)

                cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGCCGCTTGCCgctcag                                                                                                                                       <  2:2093265/150‑6 (MQ=38)

                cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAg                                                                                                                                       <  2:5570583/150‑1 (MQ=255)

                 gcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAgg                                                                                                                                      >  1:2647100/1‑150 (MQ=255)

                  cgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                     <  1:1865593/150‑1 (MQ=255)

                  cgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                     <  1:3535706/150‑1 (MQ=255)

                   gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                    >  1:898911/1‑150 (MQ=255)

                   gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                    >  1:2493713/1‑150 (MQ=255)

                   gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCACTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGCCGCTTGCCGCTCAggta                                                                                                                                    <  2:4371063/150‑3 (MQ=38)

                      gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacgac                                                                                                                                 >  2:4432035/1‑150 (MQ=255)

                               tGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAg                                                                                                                        <  2:2362557/150‑1 (MQ=255)

                               tGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAg                                                                                                                        >  1:1811481/1‑150 (MQ=255)

                                    ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGc                                                                                                                                                     <  1:2359627/116‑1 (MQ=255)

                                    ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                   <  2:61592/150‑1 (MQ=255)

                                    ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                   <  1:6449167/150‑1 (MQ=255)

                                    ggCAGAAGGTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGCTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                   <  1:738242/150‑1 (MQ=255)

                                     gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                                  <  2:4249034/150‑1 (MQ=255)

                                     gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                                  <  1:5637966/150‑1 (MQ=255)

                                     gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                                  <  2:2231778/150‑1 (MQ=255)

                                      cAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGc                                                                                                                 <  1:3998838/150‑1 (MQ=255)

                                         aaGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATa                                                                                                              <  1:5456508/150‑1 (MQ=255)

                                          aGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATAc                                                                                                             <  2:41818/150‑1 (MQ=255)

                                              gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCacctggtc                                                                                                         >  1:1211823/1‑143 (MQ=255)

                                                   gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCTGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCtatc                                                                                                    <  1:3699893/150‑2 (MQ=255)

                                                   gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCtata                                                                                                    <  1:6759105/150‑1 (MQ=255)

                                                      gtgCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTc                                                                                                 >  1:3256225/1‑150 (MQ=255)

                                                          cTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAg                                                                                             >  2:1551614/1‑150 (MQ=255)

                                                           tGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                                            <  2:5950155/150‑1 (MQ=255)

                                                           tGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                                            >  1:456468/1‑150 (MQ=255)

                                                                cTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGcac                                                                                       >  2:5909150/1‑150 (MQ=255)

                                                                cTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGcac                                                                                       >  2:4662270/1‑150 (MQ=255)

                                                                cTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGcac                                                                                       >  2:3429620/1‑150 (MQ=255)

                                                                cTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGcac                                                                                       >  1:3565181/1‑150 (MQ=255)

                                                                       gCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCt                                                                                <  2:4139109/150‑1 (MQ=255)

                                                                             gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACgctca                                                                                                                                        >  1:2674824/1‑84 (MQ=255)

                                                                              ccAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTgg                                                                         <  2:534512/150‑1 (MQ=255)

                                                                               cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                        >  2:6223377/1‑150 (MQ=255)

                                                                               cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                        <  2:2514708/150‑1 (MQ=255)

                                                                                  gCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCa                                                                     <  2:6737196/150‑1 (MQ=255)

                                                                                     tCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                                  <  1:2540508/150‑1 (MQ=255)

                                                                                     tCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                                  <  1:430635/150‑1 (MQ=255)

                                                                                       gCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATg                                                                <  2:5083790/150‑1 (MQ=255)

                                                                                        cAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGa                                                               >  1:1899958/1‑150 (MQ=255)

                                                                                            cgcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGa                                                           <  2:5380557/150‑1 (MQ=255)

                                                                                            cgcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGa                                                           <  2:150596/150‑1 (MQ=255)

                                                                                            cgcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGa                                                           <  2:1158609/150‑1 (MQ=255)

                                                                                                ggATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTg                                                       >  2:1525794/1‑150 (MQ=255)

                                                                                                 gATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTgg                                                      <  2:5428874/150‑1 (MQ=255)

                                                                                                   tGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCa                                                    >  2:3689264/1‑150 (MQ=255)

                                                                                                          cTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCCAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAggg                                             <  1:5413580/150‑1 (MQ=255)

                                                                                                          cTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAggg                                             <  1:4676998/150‑1 (MQ=255)

                                                                                                           tGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGt                                            <  2:4006473/150‑1 (MQ=255)

                                                                                                               tgtgTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                                        <  1:4387750/150‑1 (MQ=255)

                                                                                                                 tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCga                                      >  2:6432770/1‑150 (MQ=255)

                                                                                                                 tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCga                                      <  2:1864553/150‑1 (MQ=255)

                                                                                                                   ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCgaga                                    >  2:4801383/1‑150 (MQ=255)

                                                                                                                   ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACCACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCAAAATGACGGAGCAGGCAGCCAGGGTGGGCgaca                                    >  1:228322/1‑148 (MQ=255)

                                                                                                                      aGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACgg                                 <  1:1551614/150‑1 (MQ=255)

                                                                                                                        ccTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                               >  1:987378/1‑150 (MQ=255)

                                                                                                                        ccTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                               >  1:3101316/1‑150 (MQ=255)

                                                                                                                        ccTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                               >  1:1291601/1‑150 (MQ=255)

                                                                                                                          tGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg                             >  1:1151901/1‑150 (MQ=255)

                                                                                                                          tGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg                             <  1:2013686/150‑1 (MQ=255)

                                                                                                                              ttACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGa                         >  1:4624758/1‑150 (MQ=255)

                                                                                                                                   tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                    >  1:3608827/1‑150 (MQ=255)

                                                                                                                                   tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                    >  1:6660154/1‑150 (MQ=255)

                                                                                                                                    gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                                   <  1:5921667/102‑1 (MQ=255)

                                                                                                                                    gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCt                   <  1:5197979/150‑1 (MQ=255)

                                                                                                                                     cgCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCtt                  <  1:4060205/149‑1 (MQ=255)

                                                                                                                                     aGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCtt                  <  1:4939310/150‑1 (MQ=255)

                                                                                                                                       cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGa                <  1:3383808/150‑1 (MQ=255)

                                                                                                                                        gtcAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt               >  2:1943496/3‑150 (MQ=255)

                                                                                                                                        gtcAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt               >  1:4462304/3‑150 (MQ=255)

                                                                                                                                           aaGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtctc            >  2:481761/1‑150 (MQ=255)

                                                                                                                                           aaGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtctc            <  1:3801160/150‑1 (MQ=255)

                                                                                                                                             gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc          <  1:666810/150‑1 (MQ=255)

                                                                                                                                             gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc          <  1:481761/150‑1 (MQ=255)

                                                                                                                                               ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCt                                                        <  2:688018/102‑1 (MQ=255)

                                                                                                                                               ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCtt        <  2:130290/150‑1 (MQ=255)

                                                                                                                                               ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTGTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCtt        >  2:4900548/1‑150 (MQ=255)

                                                                                                                                                 cggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGa      <  2:4065453/150‑1 (MQ=255)

                                                                                                                                                   gcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCCAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGAcc    <  2:4514042/150‑1 (MQ=255)

                                                                                                                                                    cTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCg   >  2:335487/1‑150 (MQ=255)

                                                                                                                                                    cTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCg   >  1:954245/1‑150 (MQ=255)

                                                                                                                                                     tGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGt  >  1:3059566/1‑150 (MQ=255)

CCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGT  >  NC_002947/1836622‑1836920

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