BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	672,245	C→A	5.6%	V847L (GTG→TTG) ‡	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	672,245	0	C	A	5.6%	376.7 / 10.1	125	V847L (GTG→TTG) ‡	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (53/65); new base A (4/3); total (57/68)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

Only 100 of 125 total aligned reads displayed.

AGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAG  >  NC_000913/672097‑672390

aggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt                                                                                                                                                  >  2:1441034/1‑150 (MQ=255)

   catcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGATTACCTGGTAAGTAAATCACATTa                                                                                                                                               >  1:677230/4‑150 (MQ=255)

   aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCTACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                               >  2:1839923/1‑150 (MQ=255)

   aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                               <  2:605899/150‑1 (MQ=255)

    atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc                                                                                                                                              >  2:1072713/3‑150 (MQ=255)

        gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc                                                                                                                                          <  2:846406/150‑1 (MQ=255)

        gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc                                                                                                                                          >  2:1571441/1‑150 (MQ=255)

         ggttgTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAcc                                                                                                                                         <  2:946177/146‑2 (MQ=255)

         cGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAcc                                                                                                                                         <  1:1082657/150‑2 (MQ=255)

         cGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAcc                                                                                                                                         <  1:371314/150‑2 (MQ=255)

            ggTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTATACGTACAgga                                                                                                                                      >  1:354227/1‑148 (MQ=255)

               aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACg                                                                                                                                   <  1:1843788/150‑1 (MQ=255)

                 cACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAACTCACTTTACGTACAGTAACGcc                                                                                                                                 >  1:240937/1‑150 (MQ=255)

                 cACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGcc                                                                                                                                 >  2:1855268/1‑150 (MQ=255)

                 cACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGACAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAACACTTTACGTACAGTAACGcc                                                                                                                                 >  1:350771/1‑150 (MQ=255)

                   cGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCAt                                                                                                                               >  2:748074/1‑150 (MQ=255)

                        aTGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGa                                                                                                                          <  1:661566/150‑1 (MQ=255)

                         tGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGat                                                                                                                         >  2:972593/1‑150 (MQ=255)

                          gccagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGata                                                                                                                        <  1:1126330/150‑1 (MQ=255)

                          gccagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGata                                                                                                                        >  2:1201616/1‑150 (MQ=255)

                            gcgccaCACCCGGCGGTGATTACCACCGCCAGAGATAACAACAATGGTGCCAGATATCGCACGCTTCCTCCCGCGCTTCGCCAACGCCCAGATTGAGGTGTTTACCTGGTACGTAAATCACTGTCCGTCCAGTCACGCCAGCAAGATAtt                                                                                                                      >  1:1477626/3‑150 (MQ=255)

                              gccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAtttt                                                                                                                    <  1:157499/150‑1 (MQ=255)

                               ccaCACCCGGCGGTGATTAACACCGCCAGAGATAAGAACAATGTTGAGAGATATCGCACGCTTCCTCCCGCGCTTACCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTg                                                                                                                   <  2:755028/150‑1 (MQ=255)

                               ccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTg                                                                                                                   <  2:1763210/150‑1 (MQ=255)

                                 acacCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCt                                                                                                                 <  2:1083049/150‑1 (MQ=255)

                                 aaacCCGGAGGTGATTAACACCGCCAGAGATAACAACAATGTTTCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTCAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCt                                                                                                                 <  2:1395787/148‑1 (MQ=255)

                                    cccGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTAcc                                                                                                              <  1:459415/150‑1 (MQ=255)

                                     ccGGCGGTTATTAACACCGCCAGAGATAACACCAATGTTGCCAGATATCGCACGCTTCCTCCCTCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGAAATTTTGCTACCa                                                                                                             <  2:1498175/150‑1 (MQ=255)

                                     ccGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGACAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCa                                                                                                             <  1:848707/150‑1 (MQ=255)

                                      ggggggTGTTTAACACCTCAAGAGATCACAACAATGTTGCGAGATATAGCACGCTTCCTCCCGCGGTTAGCCACCGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCCTCAAGTTATTTTGCTACCAg                                                                                                            <  2:330393/146‑1 (MQ=255)

                                          ggTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGt                                                                                                        <  2:1744369/150‑1 (MQ=255)

                                           gTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGtt                                                                                                       >  1:827453/1‑150 (MQ=255)

                                           gTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGtt                                                                                                       <  1:1750861/150‑1 (MQ=255)

                                             gATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTcc                                                                                                     <  2:1636238/150‑1 (MQ=255)

                                                tAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTgg                                                                                                  <  2:1119148/150‑1 (MQ=255)

                                                tAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTgg                                                                                                  <  1:1297237/150‑1 (MQ=255)

                                                 aaCACCGCCAGAGATAACAACAATGTTGACAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGCAGTTTACCTGCTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGc                                                                                                 >  1:269886/1‑150 (MQ=255)

                                                  acacCGCCAGAGATAACAACAATGTTTCCAGATATCGCACGCTTCCTCCCGCGCGTAGCCAACGACCAGATAGAGGAGATTACCTGGTACGTAAAGCACTTGACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGcc                                                                                                <  1:285630/150‑1 (MQ=255)

                                                  acacCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGcc                                                                                                >  1:468796/1‑150 (MQ=255)

                                                  acacCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGcc                                                                                                <  2:927056/150‑1 (MQ=255)

                                                   cacCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTACCGCGCTTAGCCAACGACCACATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCa                                                                                               <  2:662062/150‑1 (MQ=255)

                                                     ccGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACTTAAATCACTATACGTACAGTAACGCCATAAAGATATATTGCTACCAGATGTTCGTGGCCAGc                                                                                             >  2:774131/1‑150 (MQ=255)

                                                       gCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCAc                                                                                           <  1:1141602/150‑1 (MQ=255)

                                                             gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                     >  2:499123/1‑150 (MQ=255)

                                                             gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                     >  1:1191745/1‑150 (MQ=255)

                                                              aTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgcg                                                                                    <  2:179795/150‑1 (MQ=255)

                                                                 acaacaATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAc                                                                                 <  1:1441034/150‑1 (MQ=255)

                                                                  caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                                >  2:1290813/1‑150 (MQ=255)

                                                                  caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                                <  1:26088/150‑1 (MQ=255)

                                                                    acaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTg                                                                              >  2:694289/1‑150 (MQ=255)

                                                                       aTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttc                                                                           <  1:1424732/150‑1 (MQ=255)

                                                                           tGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTcc                                                                       <  2:743169/150‑1 (MQ=255)

                                                                           tGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTcc                                                                       <  1:1201616/150‑1 (MQ=255)

                                                                               aGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTg                                                                   <  2:498007/150‑1 (MQ=255)

                                                                               aGATATCGCACGCTGCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTg                                                                   <  2:1177488/150‑1 (MQ=255)

                                                                                 atatCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCg                                                                 <  1:1224478/150‑1 (MQ=255)

                                                                                  tatCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGt                                                                >  1:1671491/1‑150 (MQ=255)

                                                                                   atCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAGTGAGGAGTTTACCTGGTACGTCAATCACTTTACGTACAGTACCGCCATCAAGATAGTTTGCTACCAGATGTTCCTGGCCAGCCCGTTCGCGACCCAGTTCTGACGTTGCGTc                                                               >  2:324093/1‑150 (MQ=255)

                                                                                     cGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCa                                                             >  2:691968/1‑150 (MQ=255)

                                                                                      gCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAc                                                            >  1:668752/1‑150 (MQ=255)

                                                                                       cACGCTTCCTCCCGCGCTTAGCCAACGCCCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACGTGTTCTTCCGTTGCGTCCAcc                                                           >  2:937359/1‑150 (MQ=255)

                                                                                       cACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAcc                                                           <  1:1855268/150‑1 (MQ=255)

                                                                                       cACGCTTCCTCCCGCGCTTAGACAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGTCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAcc                                                           <  2:1624420/150‑1 (MQ=255)

                                                                                       cACGCTTCCTCACGCGCTTAGCCACCGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAcc                                                           <  2:1834147/150‑1 (MQ=255)

                                                                                         cGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCgg                                                         <  1:1855405/150‑1 (MQ=255)

                                                                                            ttCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTAGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAAc                                                      <  2:1797466/150‑1 (MQ=255)

                                                                                                tccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTg                                                  <  1:1811142/150‑1 (MQ=255)

                                                                                                  ccGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGAt                                                >  2:1073006/1‑150 (MQ=255)

                                                                                                   cgcgCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCGCTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGAtt                                               >  2:118376/1‑150 (MQ=255)

                                                                                                     ctattAGACTACGAACAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGAtttt                                             <  1:80072/147‑1 (MQ=255)

                                                                                                       cTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                           <  1:748074/150‑1 (MQ=255)

                                                                                                           gCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGACCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGAGTTTGGCACg                                       >  1:460552/1‑150 (MQ=255)

                                                                                                           gCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACg                                       >  1:467574/1‑150 (MQ=255)

                                                                                                             cAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGa                                     >  2:802655/1‑150 (MQ=255)

                                                                                                                cGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACtt                                  >  2:756684/1‑150 (MQ=255)

                                                                                                                   ccAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTAc                               >  2:683917/1‑150 (MQ=255)

                                                                                                                   ccAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTAc                               >  2:1618595/1‑150 (MQ=255)

                                                                                                                     aGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCg                             <  1:761149/150‑1 (MQ=255)

                                                                                                                       aTTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGtt                           <  1:1688379/150‑1 (MQ=255)

                                                                                                                       aTTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAAAGTAACGCCATCAAGATGTTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACTGTGATTTTGGCACTGACTTTACCGtt                           >  1:1187690/1‑150 (MQ=255)

                                                                                                                           aggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAcc                       <  1:771170/150‑1 (MQ=255)

                                                                                                                              agTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGc                    <  2:831894/150‑1 (MQ=255)

                                                                                                                                 cgttactgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacc                 >  1:1551008/8‑150 (MQ=255)

                                                                                                                                    ccTGGTACGTAAATCACTTTACGTACAGTAACGCCATCATGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACg              <  1:1571441/150‑1 (MQ=255)

                                                                                                                                    ccTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACg              <  1:1072713/150‑1 (MQ=255)

                                                                                                                                     cTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGa             <  1:1239172/150‑1 (MQ=255)

                                                                                                                                     cTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCCTTGCGTCAACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGa             >  1:133088/1‑150 (MQ=255)

                                                                                                                                      tGGTACTTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGAGAACCTGTTCTTCCGTTGAGTCCAACGGAACGGTGAATTTGTCACGTAATTTACCGTTAACCTGCACAAAGAc            >  2:1388458/1‑150 (MQ=255)

                                                                                                                                       tgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc           >  2:432690/2‑150 (MQ=255)

                                                                                                                                       ggTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc           >  1:641810/1‑150 (MQ=255)

                                                                                                                                        gTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCa          >  1:1198414/1‑150 (MQ=255)

                                                                                                                                           cGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCg       >  2:590244/1‑150 (MQ=255)

                                                                                                                                              aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGTCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTgg    <  1:972593/150‑1 (MQ=255)

                                                                                                                                              aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGCCTTTACCGTTAACCTGCACCACGACCCGCGTgg    >  1:367863/1‑150 (MQ=255)

                                                                                                                                              aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTgg    >  1:856168/1‑150 (MQ=255)

                                                                                                                                              aaaTCACTTAACCTACAGTTACGCCATCAAGATATTTTGCTACCAGAGGTTCCTGGCCAGCACGATCGCGAACCTGTTCTTCCGTTGCGGCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCgggg    <  2:269174/150‑1 (MQ=255)

                                                                                                                                               aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGa   >  1:1834855/1‑150 (MQ=255)

                                                                                                                                               aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGa   >  1:139694/1‑150 (MQ=255)

                                                                                                                                                aTCACTTTACGTTCAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAg  <  1:156601/150‑1 (MQ=255)

                                                                                                                                                aTCACTTTACGTACAGAACCGCCATCAAGATTTTTTGCTCCCCGATGTTCCTGGCCCGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGAATTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAg  <  2:354227/150‑1 (MQ=255)

AGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAG  >  NC_000913/672097‑672390

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