BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,928,802	G→T	7.5%	T13K (ACG→AAG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,928,802	0	G	T	7.5%	332.4 / 12.7	117	T13K (ACG→AAG)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base G (51/57); new base T (3/6); total (54/63)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 119 total aligned reads displayed.

TAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATC  >  NC_000913/1928654‑1928951

tAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGACCGCGTATCGTCCCTCAGCCAGTCGTAATTATCGATGCGCGTATCGCCATGAAGCGt                                                                                                                                                      >  1:323667/1‑150 (MQ=255)

tAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGAAGACCGTAACTATTTACTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCAGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCAACAAGCGt                                                                                                                                                      >  1:1146365/1‑150 (MQ=255)

     tGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCGTCAAgg                                                                                                                                                 >  1:333439/1‑150 (MQ=255)

     tGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATgg                                                                                                                                                 <  1:3999/150‑1 (MQ=255)

      gCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGc                                                                                                                                                <  1:387835/150‑1 (MQ=255)

      gCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGc                                                                                                                                                <  1:1497719/150‑1 (MQ=255)

       cGATCCAGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCTTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTGCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCg                                                                                                                                               <  1:365667/150‑1 (MQ=255)

         aTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgat                                                                                                                                             <  1:1783965/150‑3 (MQ=255)

           ccTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgatac                                                                                                                                           <  1:489485/150‑5 (MQ=255)

            cTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGTTAGGCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTgggg                                                                                                                                          <  1:147137/150‑1 (MQ=255)

            cTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTgggg                                                                                                                                          <  1:1695749/150‑1 (MQ=255)

             tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGa                                                                                                                                         >  1:1456384/1‑150 (MQ=255)

              gCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGaa                                                                                                                                        >  2:1622332/1‑150 (MQ=255)

                    cTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcggg                                                                                                                                  >  1:436895/1‑150 (MQ=255)

                    cTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcggg                                                                                                                                  <  2:1620286/150‑1 (MQ=255)

                      tgttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTTTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcg                                                                                                                                <  1:173418/150‑1 (MQ=255)

                      tgttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcg                                                                                                                                <  1:1452195/150‑1 (MQ=255)

                        ttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcggg                                                                                                                              >  1:206988/1‑149 (MQ=255)

                         tgtgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTTCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTAGCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCt                                                                                                                             <  1:759194/150‑1 (MQ=255)

                           tgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGAGCGGCttt                                                                                                                           >  2:1526181/1‑150 (MQ=255)

                            gAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTTTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCAGGCGGCtttt                                                                                                                          >  2:301440/1‑150 (MQ=255)

                            gAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                          <  2:1778341/150‑1 (MQ=255)

                            gAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                          <  1:1696036/150‑1 (MQ=255)

                              ggCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTgg                                                                                                                        >  1:683020/1‑150 (MQ=255)

                              ggCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGAGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGATAGCGGGCGGCTTTTgg                                                                                                                        >  1:424653/1‑150 (MQ=255)

                              ggCCATAACCCGATGCCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGAGCGTATCTCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTgg                                                                                                                        >  1:945928/1‑150 (MQ=255)

                               gCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGt                                                                                                                       >  1:601375/1‑150 (MQ=255)

                                ccATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTa                                                                                                                      >  2:1635125/1‑150 (MQ=255)

                                 cATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAg                                                                                                                     <  2:1818033/150‑1 (MQ=255)

                                   tCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCa                                                                                                                   <  1:483471/150‑1 (MQ=255)

                                        cGATGACCGTAACTATTTTCTTGGTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTa                                                                                                              <  2:1745261/150‑1 (MQ=255)

                                         gATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTAt                                                                                                             <  2:1024766/150‑1 (MQ=255)

                                            gACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGt                                                                                                          <  1:845566/150‑1 (MQ=255)

                                                   aCTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttctttct                                                                                                   <  1:661063/150‑1 (MQ=255)

                                                    cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                                  >  1:1405880/1‑150 (MQ=255)

                                                     tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTTGCTCAGAACGCGTATCGTCACGAAGCCAGTTGTAATTATCGATGCGCGTATCGCCATGAAGCGTCAGGGCGTCGGGTATCCGTGCGGCTTTTGGTAGCATGTTATTGTGCTTTCTgg                                                                                                 >  1:105010/1‑150 (MQ=255)

                                                     tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTgg                                                                                                 <  1:405324/150‑1 (MQ=255)

                                                     tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTgg                                                                                                 <  1:1301487/150‑1 (MQ=255)

                                                      aTTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGt                                                                                                <  2:1833525/150‑1 (MQ=255)

                                                       ttttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACTCGTATCGTCCCGCAGCCAGTAGGAATTATCGATGCGCGTATCGCCATGAACCGTCATGGCGTGGGTAATGCGGGCGGCATTTGGTAGCATGTTATTGTTCTTTCTGGtt                                                                                               >  2:1427204/1‑150 (MQ=255)

                                                         ttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                             <  2:1713146/150‑1 (MQ=255)

                                                            ttgttgCAGGTAGTCCAGGACTTCTGGCTGTGAACGCGTATCGTCCCTCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAAc                                                                                          <  1:1732885/150‑1 (MQ=255)

                                                             tgttgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACa                                                                                         >  1:142974/1‑150 (MQ=255)

                                                                tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                      >  1:194659/1‑150 (MQ=255)

                                                                tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                      >  1:1905713/1‑150 (MQ=255)

                                                                 gCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCtt                                                                                     <  2:1835701/150‑1 (MQ=255)

                                                                    ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                                  >  2:1064592/1‑150 (MQ=255)

                                                                    ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                                  >  2:1476801/1‑150 (MQ=255)

                                                                        gTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAggg                                                                              >  1:165985/1‑150 (MQ=255)

                                                                        gTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAggg                                                                              >  1:474732/1‑150 (MQ=255)

                                                                         tCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCTGTAGTAATTTTCGATGCGCGTATCGCCAGGAAGCGTCATGGCGGGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTAGAAgggg                                                                             >  1:1724540/1‑149 (MQ=255)

                                                                          ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                            <  2:763601/150‑1 (MQ=255)

                                                                          ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                            >  2:1682004/1‑150 (MQ=255)

                                                                          ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                            <  2:194659/150‑1 (MQ=255)

                                                                          ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                            <  1:1022080/150‑1 (MQ=255)

                                                                          ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                            <  1:1046429/150‑1 (MQ=255)

                                                                           cAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTgg                                                                           <  2:364039/150‑1 (MQ=255)

                                                                           cAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTgg                                                                           >  1:30115/1‑150 (MQ=255)

                                                                               aCTTCTGGCTGAGAACGCGTATCGGCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaa                                                                       <  2:729090/150‑1 (MQ=255)

                                                                                 ttCTGGCTGACAACGCGGAAGGTCCCCCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGATGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                                     <  2:316561/150‑1 (MQ=255)

                                                                                  tCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACt                                                                    >  2:269798/1‑150 (MQ=255)

                                                                                    tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCa                                                                  >  2:1682063/1‑150 (MQ=255)

                                                                                     ggCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTGTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCAc                                                                 >  2:1813883/1‑150 (MQ=255)

                                                                                         gagaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCggg                                                             <  2:1830368/150‑1 (MQ=255)

                                                                                          agaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                            >  1:655340/1‑150 (MQ=255)

                                                                                          agaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                            <  1:1756465/150‑1 (MQ=255)

                                                                                           gaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGa                                                           >  1:803702/1‑150 (MQ=255)

                                                                                           gaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGa                                                           <  1:1346619/150‑1 (MQ=255)

                                                                                             aCGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATg                                                         >  2:30432/1‑150 (MQ=255)

                                                                                             aCGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATg                                                         <  1:1558493/150‑1 (MQ=255)

                                                                                              cgcgTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGc                                                        >  1:1340332/1‑150 (MQ=255)

                                                                                              cgcgTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGc                                                        >  1:1690801/1‑150 (MQ=255)

                                                                                               gcgTATCGTCCCGCAGCCAGTAGTAAATATCGAGGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCCGGCGGCGTTTGGTAGCATGTGAGTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCa                                                       <  1:301440/150‑1 (MQ=255)

                                                                                                cgTATCGTCACGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATCAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTATGGTAGCATGTTATTTTTCTTTCTGGTTGAAACATCTTATAAGGGTGTCAAAACTAACCGGGGATGCaa                                                      >  2:469332/1‑150 (MQ=255)

                                                                                                   agcGTCCCGCAGCCAGTAGTAATAATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGTGAATGCGGGCGGCTTTTGGTAGCATGTTATTTTTCTTTCTGGTTGAAACATCTTATACGGGTGGCAAAACTCACCGGGGATGCAAGCg                                                   <  1:128900/148‑1 (MQ=255)

                                                                                                   aTCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCg                                                   <  2:311818/150‑1 (MQ=255)

                                                                                                        cccGCATCCAGTAGTAAATCTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTATCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACa                                              <  1:70697/150‑1 (MQ=255)

                                                                                                        cccGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACa                                              >  1:847887/1‑150 (MQ=255)

                                                                                                            cAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATTGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgggg                                          <  1:677361/150‑1 (MQ=255)

                                                                                                              gCCAGTAGTAATTATCGATGCGCGTATCTCGATGTAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                        <  1:744720/150‑1 (MQ=255)

                                                                                                              gCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                        <  1:1488505/150‑1 (MQ=255)

                                                                                                                cAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTc                                      >  1:335257/1‑150 (MQ=255)

                                                                                                                    agtaATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTg                                  <  1:442961/150‑1 (MQ=255)

                                                                                                                        aTTATCGATGCGCGTATCGCCATGAAGCGTCTTGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTa                              <  2:1897015/150‑1 (MQ=255)

                                                                                                                            tCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatg                          <  2:1522569/150‑1 (MQ=255)

                                                                                                                               aTGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGatgatg                       <  2:114367/150‑1 (MQ=255)

                                                                                                                                tGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGAtgatga                      <  2:1405880/150‑1 (MQ=255)

                                                                                                                                  cgcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACa                    >  1:1536191/1‑150 (MQ=255)

                                                                                                                                   gcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCAACGGGGATGCAAGCGAAACAGGGGAGTCATTTCTTAGATGATGACAg                   >  1:115906/1‑150 (MQ=255)

                                                                                                                                    cgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCGGGTTGCAACAGCTTATAAGGGTGGCAACACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCGTAGATGATGACAgg                  >  2:1776136/1‑150 (MQ=255)

                                                                                                                                      ccacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa                >  1:258573/4‑150 (MQ=255)

                                                                                                                                        tCGCCATGAAGCGTCATGGCGTGTGGAATGCGGGCGGCTTTTGGTAGCATTTTTTTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGAAAAACTCACCGGGGTTGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAAt              <  1:1897862/150‑1 (MQ=255)

                                                                                                                                        tCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCTAGCGAAACAGGGGAGTCTTTGCTTAGATGATGACAGGTAAt              <  1:1064592/150‑1 (MQ=255)

                                                                                                                                        tCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAAt              <  1:1476801/150‑1 (MQ=255)

                                                                                                                                        acGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAAt              >  2:708520/2‑150 (MQ=255)

                                                                                                                                          gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGGATGTTATTGTTCTTGCTGGTTGAAACATCTTATTAGGGTGGCAAAACTCACCGGGGATGAAAGCGAAACAGGGGAGTCATTGCTTAGATGTTGACAGGTATTgg            >  1:97230/1‑150 (MQ=255)

                                                                                                                                          gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTGTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGGGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg            <  1:1767787/150‑1 (MQ=255)

                                                                                                                                              tGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGcgcg        >  1:1327059/1‑150 (MQ=255)

                                                                                                                                                    gTCATGGCGTGGGGAATGCGTGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATc  <  1:1526181/150‑1 (MQ=255)

                                                                                                                                                    gTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATc  >  1:487409/1‑150 (MQ=255)

TAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATC  >  NC_000913/1928654‑1928951

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