BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,107,711	A→T	9.3%	T453S (ACC→TCC)	clsC →	cardiolipin synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,107,711	0	A	T	9.3%	275.5 / 13.5	100	T453S (ACC→TCC)	clsC	cardiolipin synthase C
Reads supporting (aligned to +/- strand): ref base A (47/43); new base T (2/8); total (49/51)
Fisher's exact test for biased strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 105 total aligned reads displayed.

GAGCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGG  >  NC_000913/1107563‑1107855

gagCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccac                                                                                                                                                 >  2:442560/1‑150 (MQ=255)

 agCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccacc                                                                                                                                                <  1:746184/150‑1 (MQ=255)

 agCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccacc                                                                                                                                                <  2:103962/150‑1 (MQ=255)

  gCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccacca                                                                                                                                               >  2:59705/1‑150 (MQ=255)

    gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGt                                                                                                                                             >  1:25003/1‑150 (MQ=255)

    gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGt                                                                                                                                             >  1:414243/1‑150 (MQ=255)

    gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATCATGCGGCCTCGCAGCTCCGTCTGGACAGGTGGGGACTGATCAACTCGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAgg                                                                                                                                             >  1:813375/1‑149 (MQ=255)

       aggCTGGCACAGTTTATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTc                                                                                                                                          <  2:859858/148‑1 (MQ=255)

        cGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCt                                                                                                                                         <  2:796564/150‑1 (MQ=255)

            ggCACAGTTAATTGATAAACGCTTTTTTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGaa                                                                                                                                     >  2:247319/1‑150 (MQ=255)

                 aGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                <  1:1028934/150‑1 (MQ=255)

                  gTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGt                                                                                                                               >  1:425814/1‑150 (MQ=255)

                    tAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTa                                                                                                                             <  2:836789/150‑1 (MQ=255)

                     aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                            <  2:839956/150‑3 (MQ=255)

                     aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACACGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                            <  2:1074675/150‑3 (MQ=255)

                      aTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttctt                                                                                                                           <  1:555425/150‑4 (MQ=255)

                         gATACACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:463792/150‑7 (MQ=255)

                         gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:468916/150‑7 (MQ=255)

                         gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:436751/150‑7 (MQ=255)

                         gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTc                                                                                                                        >  2:717751/1‑150 (MQ=255)

                                      ttCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGa                                                                                                           <  1:1151870/150‑1 (MQ=255)

                                       tCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGat                                                                                                          >  1:1100745/1‑150 (MQ=255)

                                        cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                         >  2:295534/1‑150 (MQ=255)

                                        cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                         >  2:406043/1‑150 (MQ=255)

                                        cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                         >  2:337650/1‑150 (MQ=255)

                                              cAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGccc                                                                                                   >  1:64696/1‑150 (MQ=255)

                                               aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                >  1:15404/1‑120 (MQ=255)

                                               aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                <  2:15404/120‑1 (MQ=255)

                                               aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGATGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCg                                                                                                  >  1:1193979/1‑150 (MQ=255)

                                                  atgatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTgg                                                                                               <  1:1526520/150‑1 (MQ=255)

                                                   tgatgCGGCCTGGCAGCTCCGTCTGGACAGGGGGGGACGGATCAACTGGGTTGATCGGCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGa                                                                                              >  2:926600/1‑150 (MQ=255)

                                                    gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                             >  2:1347181/1‑150 (MQ=255)

                                                    gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                             >  2:1389521/1‑150 (MQ=255)

                                                      tgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATg                                                                                           <  2:1524066/150‑1 (MQ=255)

                                                        cGGCCTGGCAGCTCGGTGTGGACAGGTGGGGACGTATCAAGTGTTTTGATCGTCATGCTATGAAAGTGATTATTCTCTAAAAAGAACCCGCGACAAGTTTCTGGAAGAGGGTTATTGTCAGACTGGCGTCGATATTGCCGGTGGAATTGt                                                                                         <  2:421598/150‑1 (MQ=255)

                                                        cGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGt                                                                                         <  2:181662/150‑1 (MQ=255)

                                                         ggCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGtt                                                                                        >  1:1168541/1‑150 (MQ=255)

                                                           ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                      >  2:24928/1‑150 (MQ=255)

                                                           ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                      <  1:510239/150‑1 (MQ=255)

                                                            cTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAtt                                                                                     <  1:59705/150‑1 (MQ=255)

                                                               gCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTa                                                                                  <  2:351816/150‑1 (MQ=255)

                                                                cAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaa                                                                                 >  1:674201/1‑150 (MQ=255)

                                                                 aGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaaa                                                                                <  2:425814/150‑1 (MQ=255)

                                                                  gCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAt                                                                               >  1:1161829/1‑150 (MQ=255)

                                                                     ccGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAtgtg                                                                            >  2:931305/1‑150 (MQ=255)

                                                                      cGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTgg                                                                           <  1:406043/150‑1 (MQ=255)

                                                                                ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                 >  1:451556/1‑150 (MQ=255)

                                                                                ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                 >  1:175826/1‑150 (MQ=255)

                                                                                   ggggACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTgg                                                              >  2:1530244/1‑150 (MQ=255)

                                                                                        cGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCtt                                                         >  2:312167/1‑150 (MQ=255)

                                                                                        cGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCtt                                                         <  1:1408705/150‑1 (MQ=255)

                                                                                         ggATCAACTGGGTTGATCGTGATGCAAAGTATTAGATTATTCTGTAAATAGTACCGGCCACCAGTTTCTGGAAGCGGGTTATTGTCAGACTGGCGTCGATATTGCGCGTGGATTGGTTATTGTAAATGTGTCTTAACGGGCTGGTGGttt                                                        <  2:918169/150‑1 (MQ=255)

                                                                                          gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                       <  1:24928/150‑1 (MQ=255)

                                                                                          gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                       >  2:510868/1‑150 (MQ=255)

                                                                                          gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                       >  2:121147/1‑150 (MQ=255)

                                                                                             cAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    <  2:969659/150‑1 (MQ=255)

                                                                                              aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                   >  1:768311/1‑150 (MQ=255)

                                                                                              aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                   >  1:254669/1‑150 (MQ=255)

                                                                                                   ggTTGATCGTCATGCAGAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCa                                              <  2:1091747/150‑1 (MQ=255)

                                                                                                    gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAATAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGGCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                             <  2:255058/150‑1 (MQ=255)

                                                                                                    gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                             >  1:1136980/1‑150 (MQ=255)

                                                                                                        aTCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCtt                                         >  1:1525443/1‑150 (MQ=255)

                                                                                                           gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                      <  1:447682/150‑1 (MQ=255)

                                                                                                           gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                      <  1:442560/150‑1 (MQ=255)

                                                                                                              aTGCAAAGAAAGTGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                   <  1:1390920/150‑1 (MQ=255)

                                                                                                              aTGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGTAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTTGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                   <  1:1396557/150‑1 (MQ=255)

                                                                                                               tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                  <  1:1389521/150‑1 (MQ=255)

                                                                                                               tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                  >  2:707308/1‑150 (MQ=255)

                                                                                                                gCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCg                                 <  1:717751/150‑1 (MQ=255)

                                                                                                                 cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    >  1:889709/1‑130 (MQ=255)

                                                                                                                 cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    <  2:889709/130‑1 (MQ=255)

                                                                                                                 cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                >  1:224893/1‑150 (MQ=255)

                                                                                                                 cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                <  1:295534/150‑1 (MQ=255)

                                                                                                                  aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                               >  2:1184448/1‑150 (MQ=255)

                                                                                                                  aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                               <  2:1168541/150‑1 (MQ=255)

                                                                                                                    agaaagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGttt                             <  2:674201/150‑1 (MQ=255)

                                                                                                                       aagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGGTTAACAACCGGTTTAcc                          >  2:1281358/1‑150 (MQ=255)

                                                                                                                       aagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTAcc                          <  1:312167/150‑1 (MQ=255)

                                                                                                                        agatATTATTCTCAACAAAGAACCCGCCTCCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTTGCGTACCGGGCTGGTGCTTTATCGTTCTGCCGCTTAACAACCGGTTTACCa                         <  2:813375/150‑1 (MQ=255)

                                                                                                                        agagATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCa                         >  2:1077285/1‑150 (MQ=255)

                                                                                                                         gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        <  1:1132095/150‑1 (MQ=255)

                                                                                                                         gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        >  1:783326/1‑150 (MQ=255)

                                                                                                                          agaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCaga                       <  1:677336/150‑1 (MQ=255)

                                                                                                                           tattattCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  2:974262/149‑1 (MQ=255)

                                                                                                                           gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  2:25003/150‑1 (MQ=255)

                                                                                                                           gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  1:121147/150‑1 (MQ=255)

                                                                                                                                ttCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaa                 >  1:1426347/1‑150 (MQ=255)

                                                                                                                                  ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa               >  1:791520/1‑150 (MQ=255)

                                                                                                                                  ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa               >  1:1269805/1‑150 (MQ=255)

                                                                                                                                    cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc             >  2:1245960/1‑150 (MQ=255)

                                                                                                                                    cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc             >  1:548383/1‑150 (MQ=255)

                                                                                                                                     aaaaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGATGGTGCTTTATCGTTCTCACGATTAACAACCGGTTTACCAGAGAACAAAAACt            <  1:707308/150‑1 (MQ=255)

                                                                                                                                       aaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTc          <  1:247319/150‑1 (MQ=255)

                                                                                                                                        ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa         >  1:1214353/5‑150 (MQ=255)

                                                                                                                                        ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa         >  1:1374988/5‑150 (MQ=255)

                                                                                                                                           gAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTa      >  1:1445145/1‑150 (MQ=255)

                                                                                                                                             aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        >  1:462739/1‑130 (MQ=255)

                                                                                                                                             aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        <  2:462739/130‑1 (MQ=255)

                                                                                                                                              cccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACg   >  1:205512/1‑150 (MQ=255)

                                                                                                                                               ccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACgg  >  2:524583/1‑150 (MQ=255)

GAGCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGG  >  NC_000913/1107563‑1107855

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