BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	672,245	C→A	5.1%	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.1%	411.3 / 7.1	137	V847L (GTG→TTG)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (66/63); new base A (6/2); total (72/65)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 139 total aligned reads displayed.

AAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGG  >  NC_000913/672096‑672388

aaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCAc                                                                                                                                                 <  1:1190033/150‑1 (MQ=255)

aaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCAc                                                                                                                                                 <  1:1295995/150‑1 (MQ=255)

 tggcatcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACTCCGCCAGAGATACCAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCGTAGCCACCGACCAGAGTGAGGCGTTTACCTGGTCCGTAAATCACt                                                                                                                                                >  1:304158/7‑150 (MQ=255)

 tggcatcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt                                                                                                                                                >  1:1978301/7‑150 (MQ=255)

    catcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                             >  1:1430087/4‑150 (MQ=255)

    aaCCTTCGTGTTATCTCGCAGATGCCAGCCACCCCAGGCGGTGATTAACACCGCCAGATATCACAACAATCTTGCCAGATATCGCACGCGTCCTCCCGCGCTTAGCAAACTACCAGATTGAGGATTTTACCTGGTACCTAAATCACTTTa                                                                                                                                             <  2:48756/150‑1 (MQ=255)

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

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

         gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGGGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc                                                                                                                                        <  2:1809531/150‑1 (MQ=255)

           gTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg                                                                                                                                      <  1:815860/150‑1 (MQ=255)

             ggTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAccagg                                                                                                                                    <  2:1969847/150‑5 (MQ=255)

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

                     gCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAAGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTATACGTACAGTAACGCCATc                                                                                                                            >  2:937126/1‑150 (MQ=255)

                       aGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCaa                                                                                                                          <  1:1799728/150‑1 (MQ=255)

                          tGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGat                                                                                                                       <  2:1285058/150‑1 (MQ=255)

                           gccagccaCACCCGGCGGTGATTTAAACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCGTAGCCCACGCCCAGATTGAGGAGTTTACCTGGTACGTTAATCAATTTACGTACAGTAACGCCATCAAGata                                                                                                                      >  2:1833848/1‑150 (MQ=255)

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

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

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

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

                            ccagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGatat                                                                                                                     <  2:674874/150‑1 (MQ=255)

                            ccagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGatat                                                                                                                     <  2:1879352/150‑1 (MQ=255)

                             cagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTACCGCCAACAAGAAAtt                                                                                                                    >  2:369205/1‑150 (MQ=255)

                             cagccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAtt                                                                                                                    <  1:1911296/150‑1 (MQ=255)

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

                                ccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTg                                                                                                                 >  2:881508/1‑150 (MQ=255)

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

                                   cacCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTa                                                                                                              >  2:1507957/1‑150 (MQ=255)

                                   cacCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTa                                                                                                              <  1:1105236/150‑1 (MQ=255)

                                              gATTAACACCGCCAGAGATAACATCAATGTTGCGAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTcc                                                                                                   <  1:221659/150‑1 (MQ=255)

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

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

                                               tttAACACCGCGAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGACGAGTTTACCTGGAACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCt                                                                                                  <  1:1076595/149‑1 (MQ=255)

                                               ttaatcacCGCCGGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCt                                                                                                  >  1:213427/6‑150 (MQ=255)

                                               aTTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCt                                                                                                  >  2:721038/1‑150 (MQ=255)

                                                ttAACACCGCCAGAGATAACAACAATGGTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTg                                                                                                 <  2:1221061/150‑1 (MQ=255)

                                                taaaCACCGCCAGAGATAACTACAATGTTGACAGATATCGCACGCTTCCTCACGCGCTTATCAAACGACAAGATTGAGTAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACTCCATCAATATATTTTGCTACCAGATGTTCCTg                                                                                                 <  1:1551978/148‑1 (MQ=255)

                                                  aaCACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGAGATTTTGCTACCAGATGTTCCTGGc                                                                                               >  2:256305/1‑150 (MQ=255)

                                                    cacCGCCAGAGATAACAACAATGTTCCCAGATAACGCACGCTTCCTCCCCCGCTTAGCAAACGACCAGATTGAGGAGTTTACCTGATACGTAAATCACTATACGTACAGTAACGCCATCAAGATATAGACCTACCAGATGTTCCTCGCCa                                                                                             >  1:1745666/1‑150 (MQ=255)

                                                      ccGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGc                                                                                           >  2:437642/1‑150 (MQ=255)

                                                       cGCCAGAGATAACACCAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAATTTGCTACCAGATGTTCCTGGCCAGCa                                                                                          <  2:841260/150‑1 (MQ=255)

                                                       cGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCa                                                                                          >  1:1504650/1‑150 (MQ=255)

                                                        gCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCAc                                                                                         >  1:479618/1‑150 (MQ=255)

                                                             agaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcg                                                                                    <  2:1868830/150‑1 (MQ=255)

                                                             agaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcg                                                                                    <  2:1430087/150‑1 (MQ=255)

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

                                                                tAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGa                                                                                 >  2:1391467/1‑150 (MQ=255)

                                                                  acaacaATGTTGCCAGATATCGCACGCTTCCTCCCGTGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAc                                                                               <  1:1016441/150‑1 (MQ=255)

                                                                   caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCCGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                              <  1:1885300/150‑1 (MQ=255)

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

                                                                   caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCAACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                              >  1:1238011/1‑150 (MQ=255)

                                                                       aaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGtt                                                                          <  1:773365/150‑1 (MQ=255)

                                                                          gTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttctt                                                                       >  2:894076/1‑150 (MQ=255)

                                                                          gTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttctt                                                                       <  1:418893/150‑1 (MQ=255)

                                                                          gTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttctt                                                                       >  2:1933452/1‑150 (MQ=255)

                                                                          gTTGCCAGATATCGCACGCTTCCTAACGCGCTTAGCCAACGACCAGATTGAGGAGTTTACATGGTACGTAAATCACTTTACGTACAGATACGCGATCAAGATATTTTGCTACAAGATGTTCCTGGCGAGCACGTTCGCGAACCTGttctt                                                                       >  2:1806323/1‑150 (MQ=255)

                                                                            tGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTAACCTGGTACGAAAATCACTTTACGTACAGTAACGGCATCAAGATATTTTGCTACCAGATGTACCTGGCCAGAACGTTCGCGAACCTGTTCTTcc                                                                     >  1:124789/1‑150 (MQ=255)

                                                                             gCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCg                                                                    <  1:412632/150‑1 (MQ=255)

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

                                                                                 gATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGc                                                                >  1:749139/1‑150 (MQ=255)

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

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

                                                                                   tatCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTTCCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGt                                                              >  1:1288808/1‑150 (MQ=255)

                                                                                   tatCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGt                                                              <  1:665689/150‑1 (MQ=255)

                                                                                    atCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTc                                                             <  2:1787027/150‑1 (MQ=255)

                                                                                    atCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTc                                                             >  1:679021/1‑150 (MQ=255)

                                                                                     tCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTcc                                                            >  1:458889/1‑150 (MQ=255)

                                                                                      cGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCa                                                           <  1:353785/150‑1 (MQ=255)

                                                                                       gCACGCTTCCTCCCGCGCTTAGCCAACGCCCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAc                                                          >  2:1161972/1‑150 (MQ=255)

                                                                                          cGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCgg                                                       >  1:932779/1‑150 (MQ=255)

                                                                                           gCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGa                                                      >  2:1240256/1‑150 (MQ=255)

                                                                                             ttCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCGTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAAc                                                    >  1:533984/1‑150 (MQ=255)

                                                                                                ctccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGt                                                 >  2:1319/1‑150 (MQ=255)

                                                                                                ctccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGt                                                 >  1:1222394/1‑150 (MQ=255)

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

                                                                                                              cAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATTACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCATATGTTCCTGGCCAGCACGTTCGCGAACCTGTGCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGa                                   <  1:1341262/150‑1 (MQ=255)

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

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

                                                                                                                          tGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACATTAACGCCATCAAGATATTTTTCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTaa                       <  1:2005123/150‑1 (MQ=255)

                                                                                                                          tGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTaa                       <  2:1556800/150‑1 (MQ=255)

                                                                                                                           gaggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAc                      >  1:1142195/1‑150 (MQ=255)

                                                                                                                           gaggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAc                      >  2:1017009/1‑150 (MQ=255)

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

                                                                                                                                   tACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacca              <  2:1978301/150‑1 (MQ=255)

                                                                                                                                   tACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTACCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacca              >  1:516443/1‑150 (MQ=255)

                                                                                                                                   tACCGGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATCTTTTGCTCCCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacca              <  2:304158/150‑1 (MQ=255)

                                                                                                                                      cTGGTACTTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTAGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGTTGATTTTGGCACGGACTTTACCGTTAACCTGCACAACGa           >  2:588459/1‑150 (MQ=255)

                                                                                                                                      cTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCCCGa           >  2:1530609/1‑150 (MQ=255)

                                                                                                                                       tGGTACGTAAATCACTTTACGTACCGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAc          <  2:1319363/150‑1 (MQ=255)

                                                                                                                                       ctgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAc          >  1:567458/3‑150 (MQ=255)

                                                                                                                                       ctgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAc          >  1:1985846/3‑150 (MQ=255)

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

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

                                                                                                                                        ggTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGGTCCTGGCCAGCACGTTCGCGAACCGGATCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc         <  1:275820/150‑1 (MQ=255)

                                                                                                                                           aCGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTATGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGc      <  1:516362/150‑1 (MQ=255)

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

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

                                                                                                                                               aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAATCTGCACCACGACCAGCGTgg  >  2:459026/1‑150 (MQ=255)

                                                                                                                                               aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTgg  <  1:1514028/150‑1 (MQ=255)

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

AAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGG  >  NC_000913/672096‑672388

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: `‑`