BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,778,876	G→A	5.7%	G493D (GGC→GAC)	lldP →	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,778,876	0	G	A	5.7%	1395.1 / 46.3	417	G493D (GGC→GAC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand): ref base G (179/214); new base A (9/15); total (188/229)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 425 total aligned reads displayed.

TTCTTCTCGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTCA  >  NC_000913/3778728‑3779023

    tctcGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGt                                                                                                                                                  >  1:114372/1‑148 (MQ=255)

    tctcGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTc                                                                                                                                                 <  3:151013/149‑1 (MQ=255)

        gCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCAcc                                                                                                                                              <  2:875787/148‑1 (MQ=255)

           gTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGt                                                                                                                                           <  1:172595/148‑1 (MQ=255)

           gTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGt                                                                                                                                           <  2:571977/148‑1 (MQ=255)

            ttCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTaa                                                                                                                                         <  4:695986/149‑1 (MQ=255)

            ttCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGtg                                                                                                                                          <  2:311562/148‑2 (MQ=255)

              ccTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAg                                                                                                                                        <  1:414202/148‑1 (MQ=255)

                 cGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGtggtatt                                                                                                                                     <  2:193304/148‑7 (MQ=255)

                  ggctggctGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatgat                                                                                                                                   <  4:433491/149‑1 (MQ=255)

                  ggctggctGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatga                                                                                                                                    <  1:715327/148‑1 (MQ=255)

                      ggctggGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTcc                                                                                                                               <  3:239209/149‑1 (MQ=255)

                          gggggTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCg                                                                                                                            >  2:380165/1‑148 (MQ=255)

                           ggggTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCa                                                                                                                          <  3:275948/149‑1 (MQ=255)

                             ggTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCaa                                                                                                                         <  2:60547/148‑1 (MQ=255)

                             ggTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAAt                                                                                                                        <  4:104213/149‑1 (MQ=255)

                               tATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTa                                                                                                                                          <  2:458841/129‑1 (MQ=255)

                                 ttCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATcta                                                                                                                     <  1:305575/148‑1 (MQ=255)

                                  tCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGTGTCACCGGTa                                                                                                                                          <  1:477495/126‑1 (MQ=255)

                                     tGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatcgct                                                                                                                >  3:375194/1‑149 (MQ=255)

                                        ccGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatga                                                                                                                                    <  2:137288/126‑1 (MQ=255)

                                             tCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCtg                                                                                                        >  4:331411/1‑149 (MQ=255)

                                                     cTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTgcggcg                                                                                                 <  2:819107/148‑1 (MQ=255)

                                                     cTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTgcggcg                                                                                                 >  1:346581/1‑148 (MQ=255)

                                                      tCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatcgcta                                                                                                               <  2:680336/133‑1 (MQ=255)

                                                        aTCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatgat                                                                                                                                   >  4:18960/1‑111 (MQ=255)

                                                        aTCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatgat                                                                                                                                   <  3:18960/111‑1 (MQ=255)

                                                           tAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGcc                                                                                          <  3:823498/149‑1 (MQ=255)

                                                               gCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGTCTGTGCGGCGGTAGGCCtgg                                                                                       <  2:174664/148‑1 (MQ=255)

                                                                cccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCtgt                                                                                                       >  3:708720/1‑131 (MQ=255)

                                                                     ttCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCaa                                                                                >  4:728781/1‑149 (MQ=255)

                                                                      tCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCaa                                                                                <  1:511357/148‑1 (MQ=255)

                                                                       cgccgcGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCaaa                                                                               <  2:753896/148‑1 (MQ=255)

                                                                        gccgcGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAag                                                                              <  2:330017/148‑1 (MQ=255)

                                                                         ccgcGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAaga                                                                             >  2:695743/1‑148 (MQ=255)

                                                                          cgcgCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAagag                                                                            <  1:229710/148‑1 (MQ=255)

                                                                              cTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGa                                                                       <  4:244033/149‑1 (MQ=255)

                                                                                gCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGAtt                                                                     <  4:357416/149‑1 (MQ=255)

                                                                                 cAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGAtt                                                                     >  1:684967/1‑148 (MQ=255)

                                                                                  aaGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCa                                                                                                                                                >  1:612/1‑72 (MQ=255)

                                                                                   aGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTg                                                                   <  1:86938/148‑1 (MQ=255)

                                                                                    gCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGAt                                                                      <  2:426392/144‑1 (MQ=255)

                                                                                    gCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGtt                                                                 <  4:747418/149‑1 (MQ=255)

                                                                                     ccACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGtt                                                                 <  1:594466/148‑1 (MQ=255)

                                                                                     ccACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTc                                                                <  3:731162/149‑1 (MQ=255)

                                                                                     ccACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTc                                                                <  3:307095/149‑1 (MQ=255)

                                                                                     ccACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTc                                                                >  3:752186/1‑149 (MQ=255)

                                                                                      cACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTcc                                                               >  3:859553/1‑149 (MQ=255)

                                                                                         cGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAACATCATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCt                                                            <  4:784624/149‑1 (MQ=255)

                                                                                          gcagcaCAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCtt                                                           <  4:375194/149‑1 (MQ=255)

                                                                                            agcaCAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTa                                                         <  3:901172/149‑1 (MQ=255)

                                                                                            agcaCAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCATCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCtttc                                                         <  3:783794/149‑2 (MQ=255)

                                                                                              cacaACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTAc                                                        >  1:544911/1‑148 (MQ=255)

                                                                                              cacaACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACt                                                       <  3:118864/149‑1 (MQ=255)

                                                                                                  acaaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCAGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTGGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTc                                                    >  1:744347/1‑148 (MQ=255)

                                                                                                   caaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGATTCTGATTTGTTCCGCTTTACTGTCa                                                   >  2:450204/1‑148 (MQ=255)

                                                                                                   caaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCa                                                   >  1:32110/1‑148 (MQ=255)

                                                                                                     aaTTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtgg                                                                                       <  3:483852/110‑1 (MQ=255)

                                                                                                     aaTTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCaaa                                                 >  1:395361/1‑148 (MQ=255)

                                                                                                        tGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAacac                                              >  1:292474/1‑148 (MQ=255)

                                                                                                         ggCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAAcaca                                             <  2:292474/148‑1 (MQ=255)

                                                                                                           cGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGc                                           <  1:810450/148‑1 (MQ=255)

                                                                                                            gTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCt                                         <  3:530121/149‑1 (MQ=255)

                                                                                                             tctcTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTg                                        >  3:667880/1‑149 (MQ=255)

                                                                                                               tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTc                                                    >  3:443450/1‑135 (MQ=255)

                                                                                                               tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACCGCCTGa                                       <  2:497516/148‑1 (MQ=255)

                                                                                                               tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGa                                       <  1:914285/148‑1 (MQ=255)

                                                                                                               tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                      <  4:466806/149‑1 (MQ=255)

                                                                                                               tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                      >  3:808609/1‑149 (MQ=255)

                                                                                                                  gATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCt                                    >  1:654870/1‑148 (MQ=255)

                                                                                                                   aTCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTg                                                      >  1:892719/1‑129 (MQ=255)

                                                                                                                   aTCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTg                                                      <  2:892719/129‑1 (MQ=255)

                                                                                                                    tCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTc                                  >  1:126065/1‑148 (MQ=255)

                                                                                                                    tCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCa                                 >  4:219735/1‑149 (MQ=255)

                                                                                                                     cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTTCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCAc                                <  4:540521/149‑1 (MQ=255)

                                                                                                                     cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCa                                 <  2:73222/148‑1 (MQ=255)

                                                                                                                     cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCa                                 <  2:544911/148‑1 (MQ=255)

                                                                                                                     cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCAc                                <  3:341647/149‑1 (MQ=255)

                                                                                                                     cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCAc                                >  4:868446/1‑149 (MQ=255)

                                                                                                                         tGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTg                             >  1:637840/1‑148 (MQ=255)

                                                                                                                           cTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                      <  4:587763/137‑1 (MQ=255)

                                                                                                                           cTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGta                           >  1:453540/1‑148 (MQ=255)

                                                                                                                              gTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                    >  4:289989/1‑88 (MQ=255)

                                                                                                                              gTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                    <  3:289989/88‑1 (MQ=255)

                                                                                                                                tGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTg                      <  2:619491/148‑1 (MQ=255)

                                                                                                                                 gccgccAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTggg                    <  4:30787/149‑1 (MQ=255)

                                                                                                                                    gccAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTcc                                                               >  4:137213/1‑103 (MQ=255)

                                                                                                                                    gccAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTcc                                                               <  3:137213/103‑1 (MQ=255)

                                                                                                                                      tcatcttaccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGa               >  4:505332/8‑149 (MQ=255)

                                                                                                                                       aaTACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGAt              <  3:643314/149‑1 (MQ=255)

                                                                                                                                        aTACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGAt              >  1:130595/1‑148 (MQ=255)

                                                                                                                                            taccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcacc          >  2:122066/2‑148 (MQ=255)

                                                                                                                                            caccGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcacca         <  3:585053/149‑1 (MQ=255)

                                                                                                                                             accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGc                                                             >  4:664018/1‑96 (MQ=255)

                                                                                                                                             accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCt                                                            <  2:246557/97‑1 (MQ=255)

                                                                                                                                             accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCt                                                            >  1:246557/1‑97 (MQ=255)

                                                                                                                                             accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcacca         >  2:513716/1‑148 (MQ=255)

                                                                                                                                               cGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACg       <  1:734265/148‑1 (MQ=255)

                                                                                                                                                  tGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCtt    >  2:742415/1‑148 (MQ=255)

                                                                                                                                                    gCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCAcc                               >  4:843389/1‑119 (MQ=255)

TTCTTCTCGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTCA  >  NC_000913/3778728‑3779023

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