BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	672,243	C→T	5.5%	V847V (GTG→GTA) ‡	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	672,243	0	C	T	5.5%	406.1 / 7.8	126	V847V (GTG→GTA) ‡	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (54/65); new base T (4/3); total (58/68)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 128 total aligned reads displayed.

GGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAG  >  NC_000913/672094‑672390

ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATc                                                                                                                                                     >  1:972357/1‑150 (MQ=255)

ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATc                                                                                                                                                     >  1:1797466/1‑150 (MQ=255)

 gaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa                                                                                                                                                    >  2:1043727/1‑150 (MQ=255)

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

    ggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                                                                                                                                                 <  2:1830765/150‑1 (MQ=255)

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

      aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                               <  2:605899/150‑1 (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)

                 tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc                                                                                                                                    >  2:106410/1‑150 (MQ=255)

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

                         aGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCaa                                                                                                                            >  2:771170/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)

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

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

                                  ccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTg                                                                                                                   <  2:280871/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)

                                          ggcggTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGa                                                                                                           <  1:810201/150‑1 (MQ=255)

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

                                              gTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGtt                                                                                                       >  1:827453/1‑150 (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)

                                                     acacCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTACATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGcc                                                                                                >  1:349591/1‑150 (MQ=255)

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

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

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

                                                         cGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGGTCCTGGCCAGCa                                                                                            <  2:1657956/150‑1 (MQ=255)

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

                                                                gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                     <  1:717532/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                           tggtGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTTGTACGGAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttct                                                                          <  2:1133911/147‑1 (MQ=255)

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

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

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

                                                                                  aGATATCGCACGCTTCCTACCGCTCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTATATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTg                                                                   <  2:721085/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)

                                                                                            cGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCgg                                                         <  1:1544819/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)

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

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

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

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

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

                                                                                                                   cGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACtt                                  >  2:756684/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)

                                                                                                                             gaggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAc                        >  1:1657618/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)

                                                                                                                                       ccTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAGCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGCCTTTACCGTTAACCTGCACCACg              >  1:225755/1‑150 (MQ=255)

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

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

                                                                                                                                        cTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGa             <  1:1239172/150‑1 (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)

                                                                                                                                           gTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATATTGCTACAAGATGTTCATGGCCAGCACGTTCGCGAACCTGATCTTCCGTAGCGTCCACCGGAACGGTGATTTTTGCACGGACTTTACCGTTAACCTGCACCACGACCa          >  2:64802/1‑150 (MQ=255)

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

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

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

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

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

GGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAG  >  NC_000913/672094‑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: `‑`