BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,724,413	A→T	7.2%	S96C (AGC→TGC)	yfiM →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,724,413	0	A	T	7.2%	573.1 / 26.0	166	S96C (AGC→TGC)	yfiM	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (61/93); new base T (5/7); total (66/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 172 total aligned reads displayed.

ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAG  >  NC_000913/2724265‑2724561

acaGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGTTAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCaa                                                                                                                                                      <  1:142490/149‑1 (MQ=255)

acaGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCaa                                                                                                                                                      >  2:727810/1‑149 (MQ=255)

 caGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAg                                                                                                                                                     >  3:177136/1‑149 (MQ=255)

     aTCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCAcc                                                                                                                                                 <  4:224561/149‑1 (MQ=255)

       cAGGGGCTGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCgg                                                                                                                                               <  3:373774/149‑1 (MQ=255)

        aGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGc                                                                                                                                              <  4:6881/149‑1 (MQ=255)

         ggggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCGTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcg                                                                                                                                             <  2:239616/149‑2 (MQ=255)

          gggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCAc                                                                                                                                                  >  1:506529/1‑143 (MQ=255)

          gggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCta                                                                                                                                            <  3:1404454/149‑1 (MQ=255)

                  ccGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATATCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    <  3:729571/149‑1 (MQ=255)

                  ccGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    <  2:951743/149‑1 (MQ=255)

                   cGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTg                                                                                                                                   >  3:532573/1‑149 (MQ=255)

                    gggATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTgg                                                                                                                                  <  1:341868/149‑1 (MQ=255)

                       aTCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCaa                                                                                                                               <  4:1008024/149‑1 (MQ=255)

                        tCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAAc                                                                                                                              <  1:387450/149‑1 (MQ=255)

                          gCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTg                                                                                                                            >  1:914108/1‑149 (MQ=255)

                           cAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGa                                                                                                                           <  3:1382794/149‑1 (MQ=255)

                           cAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGa                                                                                                                           <  1:369683/149‑1 (MQ=255)

                                  aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCa                                                                                                                    <  2:1273146/149‑1 (MQ=255)

                                    gTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACt                                                                                                                  >  3:32614/1‑149 (MQ=255)

                                       tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTaaa                                                                                                               <  4:532573/149‑1 (MQ=255)

                                         gATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGa                                                                                                             >  1:127562/1‑149 (MQ=255)

                                           ttGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACg                                                                                                           >  3:942612/1‑149 (MQ=255)

                                               tGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAGGCATc                                                                                                       >  2:1548678/1‑149 (MQ=255)

                                                gTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATcc                                                                                                      <  3:1413299/149‑1 (MQ=255)

                                                gTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATcc                                                                                                      <  2:682328/149‑1 (MQ=255)

                                                  tctcTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATcccc                                                                                                    >  2:1651974/1‑149 (MQ=255)

                                                   ctctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCccct                                                                                                   <  4:335642/149‑1 (MQ=255)

                                                    tctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcga                                                                                                                                            <  2:882257/107‑3 (MQ=255)

                                                    tctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcga                                                                                                                                            >  1:882257/1‑105 (MQ=255)

                                                      tGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCa                                                                                                         <  1:1510087/140‑1 (MQ=255)

                                                           gTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    >  3:841582/1‑108 (MQ=255)

                                                           gTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    <  4:841582/108‑1 (MQ=255)

                                                            tttGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGc                                                                                          <  1:416222/149‑1 (MQ=255)

                                                               gggggCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGAt                                                                                       <  2:1302931/149‑1 (MQ=255)

                                                                  ggCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTa                                                                                    <  1:968407/149‑1 (MQ=255)

                                                                   gCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAg                                                                                   <  4:27421/149‑1 (MQ=255)

                                                                     gTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    >  1:9982/1‑98 (MQ=255)

                                                                     gTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCa                                                                                 <  4:1624818/149‑1 (MQ=255)

                                                                      tCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCAt                                                                                >  3:591620/1‑149 (MQ=255)

                                                                        aaaaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                            >  4:400105/1‑119 (MQ=255)

                                                                        aaaaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                            <  3:400105/119‑1 (MQ=255)

                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  4:742071/149‑1 (MQ=255)

                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  4:1609324/149‑1 (MQ=255)

                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  2:864421/149‑1 (MQ=255)

                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCGGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:28326/149‑1 (MQ=255)

                                                                               cTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAc                                                                       <  3:442981/149‑1 (MQ=255)

                                                                               cTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAc                                                                       >  4:466158/1‑149 (MQ=255)

                                                                               cTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAc                                                                       <  1:1599374/149‑1 (MQ=255)

                                                                                tttGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAcc                                                                      >  2:1205927/1‑149 (MQ=255)

                                                                                    ggATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAAc                                                                                                                              <  2:386994/89‑1 (MQ=255)

                                                                                      aTAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAAc                                                                <  4:389812/149‑1 (MQ=255)

                                                                                      aTAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAAc                                                                <  2:817120/149‑1 (MQ=255)

                                                                                      aTAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAAc                                                                <  2:127562/149‑1 (MQ=255)

                                                                                         gccgccCCGAAGGGGGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCCATGTAACATCAAAGAAACCAGAAACgac                                                             <  4:1157679/149‑3 (MQ=255)

                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACg                                                       >  1:1004792/1‑149 (MQ=255)

                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACg                                                       >  3:238108/1‑149 (MQ=255)

                                                                                                cGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCcaccac                                                        >  1:830479/1‑147 (MQ=255)

                                                                                                cGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTG‑AAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGc                                                     >  1:1455703/1‑149 (MQ=255)

                                                                                                  aaGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTGTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGcc                                                    <  4:177136/149‑1 (MQ=255)

                                                                                                        gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAg                                              <  3:466158/149‑1 (MQ=255)

                                                                                                           gctggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                           >  2:1380133/1‑149 (MQ=255)

                                                                                                           gctggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACACGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAAc                                                                >  1:635629/1‑128 (MQ=255)

                                                                                                            ctggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGt                                                                                                                      <  2:1468012/73‑1 (MQ=255)

                                                                                                              ggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTcaca                                        <  1:1694215/149‑1 (MQ=255)

                                                                                                                agctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATa                                      <  4:739138/149‑1 (MQ=255)

                                                                                                                     gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttccct                                                                                              >  3:1122483/1‑88 (MQ=255)

                                                                                                                     gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttccct                                                                                              <  4:1122483/88‑1 (MQ=255)

                                                                                                                        ggATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTAACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaaga                              >  1:1104657/1‑149 (MQ=255)

                                                                                                                         gATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAgaagaa                             <  2:594921/149‑1 (MQ=255)

                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            <  4:238108/149‑1 (MQ=255)

                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            <  2:328722/149‑1 (MQ=255)

                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            <  3:744128/149‑1 (MQ=255)

                                                                                                                           tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                           <  3:1621547/149‑1 (MQ=255)

                                                                                                                           tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                           <  4:1327669/149‑1 (MQ=255)

                                                                                                                             tGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCt                         <  4:526514/149‑1 (MQ=255)

                                                                                                                             tGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCt                         <  4:76016/149‑1 (MQ=255)

                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaag                               <  1:1058146/141‑1 (MQ=255)

                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaag                               >  2:1058146/1‑141 (MQ=255)

                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGt                       >  1:1698228/1‑149 (MQ=255)

                                                                                                                                 cTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGttt                     <  2:1698228/149‑1 (MQ=255)

                                                                                                                                 cTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGt                       >  1:127613/1‑147 (MQ=255)

                                                                                                                                  tGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGc                                                              >  3:1205773/1‑107 (MQ=255)

                                                                                                                                  tGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCGCATACCAGAAGAAGGCTGTTTc                    >  3:1036259/1‑149 (MQ=255)

                                                                                                                                   gggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACgg                                                      <  4:1565560/114‑1 (MQ=255)

                                                                                                                                   gggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTcc                   >  1:1627515/1‑149 (MQ=255)

                                                                                                                                      gtatagCCGGTGCTTGCACCGGCTATACCTTCTGGCAACTGACCCGTCACGAAAGACGCATCCCCTTCCCTTCGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCTGCAAGGTCACATACCAGAAGAATGCTGGTTCCAtt                >  4:651045/6‑149 (MQ=255)

                                                                                                                                        gTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTcc              <  2:408571/149‑1 (MQ=255)

                                                                                                                                           gCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATc                                                                               <  2:1581158/81‑1 (MQ=255)

                                                                                                                                           gCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATc                                                                               >  1:1581158/1‑81 (MQ=255)

                                                                                                                                           gCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            >  4:1213052/1‑132 (MQ=255)

                                                                                                                                           gCCGGTGCAAGCACCGGCTAGACCGTCTGGCAACTGACCCGTCACTAAAGTCGCATCCCCATCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt           <  1:1392740/149‑1 (MQ=255)

                                                                                                                                            ccagTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCATAACCGCCACCACGGCCATCAAGGTCACATATCAGAAGAAGGCTGTTTCCACTCCTAtt          <  4:1308667/146‑1 (MQ=255)

                                                                                                                                            ccGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAtt          <  3:651045/149‑1 (MQ=255)

                                                                                                                                            ccGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAtt          <  1:495449/149‑1 (MQ=255)

                                                                                                                                               gTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAt       <  3:1535342/149‑1 (MQ=255)

                                                                                                                                               gTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAt       <  3:1281744/149‑1 (MQ=255)

                                                                                                                                               gTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAt       <  1:164480/149‑1 (MQ=255)

                                                                                                                                                tGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAtt      >  1:786669/1‑149 (MQ=255)

                                                                                                                                                    aGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAg  >  1:1315893/1‑149 (MQ=255)

ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAG  >  NC_000913/2724265‑2724561

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