BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,122	C→T	27.7%	R34H (CGC→CAC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,122	0	C	T	27.7%	427.2 / 196.2	263	R34H (CGC→CAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (96/93); new base T (31/41); total (128/135)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 271 total aligned reads displayed.

CTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAA  >  NC_000913/4116974‑4117269

cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc                                                                                                                                                      <  3:483966/149‑1 (MQ=255)

             ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt                                                                                                                                         <  2:13924/149‑1 (MQ=255)

               gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCCGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                                                                                       >  2:2342847/1‑149 (MQ=255)

               gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTaa                                                                                                                                       <  1:1429049/149‑1 (MQ=255)

                aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                                                                                                      <  4:1138246/149‑1 (MQ=255)

                    cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata                                                                                                                                  >  4:789612/1‑149 (MQ=255)

                      gaaaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTATACCAGCCTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGACGTGTTCTGCCCAACCTGGTGTTGGGTAGATTTGATCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                                                                                >  2:59278/3‑149 (MQ=255)

                      gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                                                                                                                                            <  1:1837815/137‑1 (MQ=255)

                      gCAATTTGAACGGAACTTATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATGGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                                                                                <  4:262183/149‑1 (MQ=255)

                       cAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg                                                                                                                               <  4:1538310/149‑1 (MQ=255)

                          tttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa                                                                                                                            <  3:1374958/149‑1 (MQ=255)

                            tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATc                                                                                                                          <  3:1017948/149‑1 (MQ=255)

                            tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc                                                                                                                          <  4:1039931/149‑1 (MQ=255)

                              aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                                                                                        >  2:1006329/1‑149 (MQ=255)

                                cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                                      <  3:13465/149‑1 (MQ=255)

                                cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                                      <  2:998972/149‑1 (MQ=255)

                                  gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                                                                                                                    >  4:868973/1‑149 (MQ=255)

                                    aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCa                                                                                                                  <  3:1275906/149‑1 (MQ=255)

                                    aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCa                                                                                                                  <  4:2005532/149‑1 (MQ=255)

                                     cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCAt                                                                                                                 <  4:938513/149‑1 (MQ=255)

                                       gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTT‑GGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                                              >  2:2145411/1‑149 (MQ=255)

                                         tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACg                                                                                                             >  1:1161763/1‑149 (MQ=255)

                                           tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGAc                                                                                                           >  3:706727/1‑149 (MQ=255)

                                                  tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAACCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                                                                                    <  2:847277/149‑1 (MQ=255)

                                                     ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                                 >  3:760324/1‑149 (MQ=255)

                                                      cAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                                                >  4:35763/1‑149 (MQ=255)

                                                          aCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGCCACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggt                                                                                           >  1:1792687/1‑149 (MQ=255)

                                                          aCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggt                                                                                           <  4:635164/149‑1 (MQ=255)

                                                            ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggtgg                                                                                         >  2:1818222/1‑149 (MQ=255)

                                                            ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggtgg                                                                                         >  2:737420/1‑149 (MQ=255)

                                                            ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGATCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgcg                                                                                                     >  3:1519333/1‑138 (MQ=255)

                                                               tACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGc                                                                                      <  4:760324/149‑1 (MQ=255)

                                                                aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGcc                                                                                     <  4:1047609/149‑1 (MQ=255)

                                                                aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGcc                                                                                     >  4:1008966/1‑149 (MQ=255)

                                                                 ccAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGccc                                                                                    >  3:2231027/1‑149 (MQ=255)

                                                                 ccAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGccc                                                                                    >  3:172684/1‑149 (MQ=255)

                                                                  cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggt                                                                                           <  2:1883023/141‑1 (MQ=255)

                                                                  cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCt                                                                                   >  4:1293800/1‑149 (MQ=255)

                                                                        ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGa                                                                             <  1:1580115/149‑1 (MQ=255)

                                                                          aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAgc                                                                           <  1:638185/149‑1 (MQ=255)

                                                                          aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAgc                                                                           <  1:2449149/149‑1 (MQ=255)

                                                                           gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAgcg                                                                          >  2:2483925/1‑149 (MQ=255)

                                                                             tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                                                                                <  4:2108497/94‑1 (MQ=255)

                                                                                  gTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGat                                                                   >  1:1705612/1‑149 (MQ=255)

                                                                                      cccAGTTTTCCATTGGGTCGTGTTCTACTAAAACTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGGCCTGGTCGAGCGCAACGatatat                                                               >  2:776486/1‑149 (MQ=255)

                                                                                       ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                              >  3:1163633/1‑149 (MQ=255)

                                                                                       ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                                                                                 <  4:2415593/67‑1 (MQ=255)

                                                                                           aTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATAtttttt                                                          <  1:393154/149‑1 (MQ=255)

                                                                                               ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCTTTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                      <  1:2059643/149‑1 (MQ=255)

                                                                                               ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                      <  4:1657290/149‑1 (MQ=255)

                                                                                               ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                      <  2:1290320/149‑1 (MQ=255)

                                                                                               ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                      <  1:1268852/149‑1 (MQ=255)

                                                                                               ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                      <  3:1909206/149‑1 (MQ=255)

                                                                                                 aTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGt                                                    <  3:238298/149‑1 (MQ=255)

                                                                                                    gggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCAt                                                 <  4:1622474/149‑1 (MQ=255)

                                                                                                       tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGtt                                              <  3:2160136/149‑1 (MQ=255)

                                                                                                         gtgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTa                                            >  2:1144065/1‑149 (MQ=255)

                                                                                                         gtgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTa                                            >  3:445913/1‑149 (MQ=255)

                                                                                                          tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                                           <  3:366796/149‑1 (MQ=255)

                                                                                                          tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                                           <  1:816205/149‑1 (MQ=255)

                                                                                                          tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                                           <  1:1084200/149‑1 (MQ=255)

                                                                                                          tgtTCTACCCAACCTGGTTTTGGGTAGAGTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCC‑GGTGGTGCCCTGGTAGACCTCAACGATTTCCCTTTTCGCAGTCATGTCTaa                                           <  2:1185299/149‑1 (MQ=255)

                                                                                                               taaccAAGCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACAGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCTGTCATGTTTAATTGTc                                      <  4:48703/146‑1 (MQ=255)

                                                                                                               tACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc                                      <  4:238058/149‑1 (MQ=255)

                                                                                                                 cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc                                      <  1:1984649/147‑1 (MQ=255)

                                                                                                                 cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc                                    <  4:1506805/149‑1 (MQ=255)

                                                                                                                    aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                                 <  1:977949/149‑1 (MQ=255)

                                                                                                                     aCCTGGTTTTGGGTAGATTTGCTCAAGTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAg                                >  1:2368557/1‑149 (MQ=255)

                                                                                                                       cTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTc                              >  4:1341360/1‑149 (MQ=255)

                                                                                                                       cTGGTTTTGGGTAGATTTGCTCAAATTCGGGTGCCGCAACGCAATTGAAATTGGCATCGTGATCAATTACAACGGCGC‑GGAGCTGGGTGGTGCCCTGGTCAGCCGCAACATTATTTTTTTTTTCATTAATGTTTATTTGTCCGTTAGTc                              >  1:2035289/1‑149 (MQ=255)

                                                                                                                        tGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCa                             >  1:498928/1‑149 (MQ=255)

                                                                                                                          gTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCata                           >  1:887286/1‑149 (MQ=255)

                                                                                                                           ttttGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGTGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat                          >  3:1685582/1‑149 (MQ=255)

                                                                                                                           ttttGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat                          <  3:1842092/149‑1 (MQ=255)

                                                                                                                           ttttGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat                          <  3:1641646/149‑1 (MQ=255)

                                                                                                                               gggTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACa                      <  1:2256489/149‑1 (MQ=255)

                                                                                                                                 gTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGa                                                                    <  1:321021/101‑1 (MQ=255)

                                                                                                                                 gTAGATTTGCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGa                                                                    >  2:321021/1‑101 (MQ=255)

                                                                                                                                 gTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTa                                            <  3:2482634/125‑1 (MQ=255)

                                                                                                                                  tAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCg                   >  1:2519712/1‑149 (MQ=255)

                                                                                                                                    gATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCTTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGaa                 <  4:172684/149‑1 (MQ=255)

                                                                                                                                    gATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGaa                 <  1:1006329/149‑1 (MQ=255)

                                                                                                                                    gATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGaa                 <  3:1607617/149‑1 (MQ=255)

                                                                                                                                     aTTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAg                >  2:2310444/1‑149 (MQ=255)

                                                                                                                                         gCTCAAATTCGTGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttt            <  4:1885258/149‑1 (MQ=255)

                                                                                                                                         gCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttt            >  2:2086878/1‑149 (MQ=255)

                                                                                                                                           tCAAATTCTCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          >  3:2288359/1‑149 (MQ=255)

                                                                                                                                           tCAAATTCGGGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          <  2:1561387/149‑1 (MQ=255)

                                                                                                                                           tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          <  2:2035289/149‑1 (MQ=255)

                                                                                                                                           tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          <  4:2241621/149‑1 (MQ=255)

                                                                                                                                           tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          >  2:557800/1‑149 (MQ=255)

                                                                                                                                           tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt          >  2:2250636/1‑149 (MQ=255)

                                                                                                                                             aaaTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                                                <  3:800734/46‑1 (MQ=255)

                                                                                                                                             aaaTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                                                >  4:800734/1‑46 (MQ=255)

                                                                                                                                             aaaTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt‑ggtg                                                                                          >  3:944673/1‑67 (MQ=255)

                                                                                                                                             aaaTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGc               >  3:2069559/1‑142 (MQ=255)

                                                                                                                                              aaTTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGtt       <  3:1100115/149‑1 (MQ=255)

                                                                                                                                                ttCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCt     >  3:2433265/1‑149 (MQ=255)

                                                                                                                                                   gcgcTGCGACACGCTAATGATATTGGCATCGTGATCGATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGaa  <  2:1246358/149‑1 (MQ=255)

                                                                                                                                                   gcgcTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGaa  <  2:1619188/149‑1 (MQ=255)

CTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCT‑GGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAA  >  NC_000913/4116974‑4117269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`