BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,337,837:1	+T	6.0%	coding (998/1047 nt)	PP_2934 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,337,837	1	.	T	6.0%	471.6 / 22.9	117	coding (998/1047 nt)	PP_2934	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (56/54); new base T (5/2); total (61/56)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 128 total aligned reads displayed.

CGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAATGGCCAAG  >  NC_002947/3337688‑3337986

 gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGCGCTGGCGACCCCAGCGcctg                                                                                                                                                      <  1:534673/150‑3 (MQ=255)

  tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑gg                                                                                                                                                     <  2:564013/150‑1 (MQ=255)

  tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑gg                                                                                                                                                     >  1:1521198/1‑150 (MQ=255)

      cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGa                                                                                                                                                 <  2:3089/150‑1 (MQ=255)

          gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGt                                                                                                                                             >  2:1318422/1‑150 (MQ=255)

          gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGt                                                                                                                                             >  2:228855/1‑150 (MQ=255)

            ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAATGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTc                                                                                                                                           >  2:1069523/1‑150 (MQ=255)

             cACTGTCATGCCGAGCAGTACGCGCGAATGCCGGCGCTGGGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGTGCGTTCa                                                                                                                                          <  2:1344389/150‑1 (MQ=255)

              aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑g                                                                                                                                                      <  2:1599488/137‑1 (MQ=255)

              aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑g                                                                                                                                                      >  1:1599488/1‑137 (MQ=255)

              aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAg                                                                                                                                         <  1:1123934/150‑1 (MQ=255)

               cTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGc                                                                                                                                        <  2:2036110/150‑1 (MQ=255)

                 gCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTg                                                                                                                                      <  1:228855/150‑1 (MQ=255)

                     tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGAc                                                                                                                                  >  2:1817228/1‑150 (MQ=255)

                         gAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgc                                                                                                                              >  1:1089931/1‑150 (MQ=255)

                               tatgcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGGC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccg                                                                                                                        >  1:953352/4‑150 (MQ=255)

                                aCGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgc                                                                                                                       <  2:1089931/150‑1 (MQ=255)

                                  gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                                          >  2:43016/1‑129 (MQ=255)

                                  gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                                          <  1:43016/129‑1 (MQ=255)

                                  gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                                     >  2:899390/1‑150 (MQ=255)

                                  gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                                     >  1:1799052/1‑150 (MQ=255)

                                         tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                                              >  2:222648/1‑150 (MQ=255)

                                           ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCTATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                                            <  1:2068151/150‑1 (MQ=255)

                                           ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                                            >  2:1503000/1‑150 (MQ=255)

                                            cGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                                           >  2:1630088/1‑150 (MQ=255)

                                             ggCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACg                                                                                                          >  1:915163/1‑150 (MQ=255)

                                                    tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc                                                                                                   <  2:72298/150‑1 (MQ=255)

                                                         ggACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAccc                                                                                              >  2:1362265/1‑150 (MQ=255)

                                                               cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAAAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGTTAGGCCGACCCCGGCCa                                                                                        >  2:1239473/1‑150 (MQ=255)

                                                               cGACGGCGTCGACCGCCGGTACATACACGAATCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                                        >  1:65048/1‑150 (MQ=255)

                                                                gACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCaa                                                                                       >  1:2032291/1‑150 (MQ=255)

                                                                 aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCggactac                                                                                      <  1:1069523/150‑6 (MQ=255)

                                                                 aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                                      <  2:654210/150‑1 (MQ=255)

                                                                      gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                                                 >  2:1867665/1‑150 (MQ=255)

                                                                       tCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGACAAGGCGAGc                                                                                >  1:1844670/1‑150 (MQ=255)

                                                                          aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCaaa                                                                             <  2:1844670/150‑3 (MQ=255)

                                                                          aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                             <  2:2032291/150‑1 (MQ=255)

                                                                           ccgccgGTATATCCATACTTCGCTGCAGCGAGAACCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTGCCCAGCGGG‑GGCTGGGCGTACACCCGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAggg                                                                            <  2:2407781/150‑1 (MQ=255)

                                                                            cgccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGc                                                                           >  2:989937/1‑150 (MQ=255)

                                                                              ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAg                                                                                                      <  1:1324228/121‑1 (MQ=255)

                                                                              ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAg                                                                         >  2:1698564/1‑150 (MQ=255)

                                                                               cgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGc                                                                        <  2:915163/150‑1 (MQ=255)

                                                                                ggTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGcc                                                                       <  2:836347/150‑1 (MQ=255)

                                                                                  tataTCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCAc                                                                     <  2:410280/150‑1 (MQ=255)

                                                                                     atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                                  <  2:1521198/150‑1 (MQ=255)

                                                                                     atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                                  >  1:404464/1‑150 (MQ=255)

                                                                                      tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCt                                                                                                                                       >  2:1382383/1‑80 (MQ=255)

                                                                                      tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCt                                                                                                                                       <  1:1382383/80‑1 (MQ=255)

                                                                                        catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                                           <  1:76482/106‑1 (MQ=255)

                                                                                        catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                               >  1:2286183/1‑150 (MQ=255)

                                                                                        catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                               >  2:908780/1‑150 (MQ=255)

                                                                                        catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                                  <  1:1631759/147‑1 (MQ=255)

                                                                                        catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                                  >  2:1631759/1‑147 (MQ=255)

                                                                                         atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                                     <  2:908418/111‑1 (MQ=255)

                                                                                         atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                                     >  1:908418/1‑111 (MQ=255)

                                                                                         atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaag                                                              >  2:1825030/1‑150 (MQ=255)

                                                                                          tcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                                  >  1:588132/1‑97 (MQ=255)

                                                                                            atTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATc                                                           <  1:1698564/150‑1 (MQ=255)

                                                                                             ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                                               <  2:830791/81‑1 (MQ=255)

                                                                                             ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                                               >  1:830791/1‑81 (MQ=255)

                                                                                             ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                                      >  1:1836175/1‑122 (MQ=255)

                                                                                             ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAg                                                                              >  1:1864116/1‑130 (MQ=255)

                                                                                             ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAg                                                                              <  2:1864116/130‑1 (MQ=255)

                                                                                              tCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                                         <  1:205528/150‑1 (MQ=255)

                                                                                               cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                                        <  1:899390/150‑1 (MQ=255)

                                                                                               cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                                        <  1:1503000/150‑1 (MQ=255)

                                                                                               cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                                        <  1:156242/150‑1 (MQ=255)

                                                                                               cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCACTTGAGACGTGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                                        <  2:83296/150‑1 (MQ=255)

                                                                                                     aGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAc                                                  <  1:660959/150‑1 (MQ=255)

                                                                                                      tgtagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGTCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGCCAGCCACGCAGAAGATCCAGGCATAcc                                                 >  2:1150834/4‑150 (MQ=255)

                                                                                                      gCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcc                                                 <  2:1511829/150‑1 (MQ=255)

                                                                                                           aaaCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCT‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                                             >  1:2188945/1‑149 (MQ=255)

                                                                                                           aaaCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCT‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                                             <  2:2188945/149‑1 (MQ=255)

                                                                                                              cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGcc                                                                                                                                                        <  1:655008/40‑1 (MQ=255)

                                                                                                              cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGcc                                                                                                                                                        >  2:655008/1‑40 (MQ=255)

                                                                                                               ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                                              <  1:1218018/80‑1 (MQ=255)

                                                                                                               ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                                              >  2:1218018/1‑80 (MQ=255)

                                                                                                                  aGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCg                                     >  1:142830/1‑150 (MQ=255)

                                                                                                                   tgaacgctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                                    >  2:1976105/8‑150 (MQ=37)

                                                                                                                   gCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                                    <  1:1817228/150‑1 (MQ=255)

                                                                                                                     gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGAC‑GGCTGAGCCGAGAGGCTCAACAAGATCGTCTCTCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCa                                  >  1:1855359/1‑150 (MQ=255)

                                                                                                                        gctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACgg                               >  1:296251/3‑150 (MQ=38)

                                                                                                                           cAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgcgc                            <  1:1630088/150‑1 (MQ=255)

                                                                                                                             gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                                     <  1:2186076/59‑1 (MQ=255)

                                                                                                                             gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                                     >  2:2186076/1‑59 (MQ=255)

                                                                                                                              ccAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa                         >  1:809130/1‑150 (MQ=255)

                                                                                                                               ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAgg                                                       <  2:1875595/118‑1 (MQ=38)

                                                                                                                               ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAgg                                                       >  1:1875595/3‑120 (MQ=38)

                                                                                                                               ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa                         >  2:1703486/3‑150 (MQ=255)

                                                                                                                               cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                               <  2:94999/127‑1 (MQ=255)

                                                                                                                                   cgcTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCAcc                    <  1:475662/150‑1 (MQ=255)

                                                                                                                                           tCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                                      <  2:886259/124‑1 (MQ=255)

                                                                                                                                           accccAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCgg             <  2:296251/147‑1 (MQ=255)

                                                                                                                                            cGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaag                                                              >  2:315369/1‑99 (MQ=255)

                                                                                                                                            cGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaag                                                              <  1:315369/99‑1 (MQ=255)

                                                                                                                                              ccAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa                         <  2:1392375/134‑1 (MQ=255)

                                                                                                                                              ccAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa                         >  1:1392375/1‑134 (MQ=255)

                                                                                                                                                    cc‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAATGGCCaa   >  2:606790/1‑150 (MQ=255)

                                                                                                                                                    cc‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGACATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAATGGCCaa   <  1:908780/150‑1 (MQ=255)

                                                                                                                                                     c‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAATGGCCAAg  <  1:653149/150‑1 (MQ=255)

CGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAATGGCCAAG  >  NC_002947/3337688‑3337986

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