BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,423,621	C→G	9.1%	noncoding (35/76 nt)	thrV ←	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,621	0	C	G	9.1%	337.6 / 19.7	108	noncoding (35/76 nt)	thrV	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base C (53/44); new base G (4/6); total (58/50)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 109 total aligned reads displayed.

GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGAC  >  NC_000913/3423472‑3423764

ggCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAc                                                                                                                                                 <  2:1222132/150‑1 (MQ=255)

  cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCa                                                                                                                                               <  1:1716897/150‑1 (MQ=255)

  cgcTACGGCGTTTCACTTCTGAGTTCGGCAGGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCTAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCa                                                                                                                                               <  1:122377/150‑1 (MQ=255)

     tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGGCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAgg                                                                                                                                            <  2:263295/150‑1 (MQ=255)

      atggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAATAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                           <  1:57560/148‑1 (MQ=255)

      aCGGCGTCTCACTTCTGAGTTCGGCATTGGGCCATGCGGGACCACCGCGCTACTGGCGCCAGACAAATTCTATGCTACTCTTCCTAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                           <  1:310204/150‑1 (MQ=38)

          cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgag                                                                                                                                       <  2:639631/150‑3 (MQ=255)

           gTTTCACTTCTGAGTGCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgagg                                                                                                                                      <  2:737473/150‑4 (MQ=37)

           gTTGCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgagg                                                                                                                                      <  2:1409215/150‑4 (MQ=37)

            tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgaggt                                                                                                                                     <  1:444765/150‑5 (MQ=255)

            tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgaggt                                                                                                                                     <  1:1009285/150‑5 (MQ=255)

             ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCGCACCCTTGGTAAGGGTgaggtc                                                                                                                                    <  1:1384617/150‑6 (MQ=37)

             ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGAAACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctc                                                                                                                                    <  1:616902/150‑1 (MQ=255)

              tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTTTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCCTGGTAAGGGTgaggtcc                                                                                                                                   <  1:11853/150‑7 (MQ=25)

               cACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTa                                                                                                                                  <  1:1272025/150‑1 (MQ=255)

                aCTTCTGAGTTCGGCATGGGGTCAGTTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGAGCTCACCCTTACCAAGGGTGCGCTCTAc                                                                                                                                 <  1:65954/150‑1 (MQ=255)

                aCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                                                                                                                 <  1:438331/150‑1 (MQ=255)

                aCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTATTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                                                                                                                 <  1:937930/150‑1 (MQ=255)

                aCTTCTGAGTGCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                                                                                                                 <  2:1507690/150‑1 (MQ=255)

                    ccgatctCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                                                                                                                             <  1:1635387/144‑1 (MQ=255)

                    cTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                                                                                                                             >  1:720432/1‑150 (MQ=255)

                        gTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACATTAACCTAAAAAGAGGTGCTGATACCCAGAGTCTAACTGGTGACCTCACCCTTACCAAGGGTGAGCTCTACCAACTGac                                                                                                                         >  1:1670004/1‑149 (MQ=38)

                           cGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCGGAACTGTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGGCCTCACGCTTAGCAAGGGTGCGCTCTGCGAACAGatcgg                                                                                                                      >  2:1635387/1‑146 (MQ=37)

                                  tagcgcggtggtCCCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGc                                                                                                               >  2:1682149/7‑150 (MQ=255)

                                  gggTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCACCTGAGCCATATCAGc                                                                                                               >  2:909388/1‑150 (MQ=255)

                                  gggTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGc                                                                                                               <  2:1185516/150‑1 (MQ=255)

                                    gTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCAc                                                                                                             <  2:670090/150‑1 (MQ=255)

                                     gcAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACg                                                                                                            <  1:940860/149‑1 (MQ=255)

                                     gcAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACg                                                                                                            <  2:1725451/149‑1 (MQ=255)

                                        ggTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTa                                                                                                         >  1:1416298/1‑150 (MQ=255)

                                              accaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTg                                                                                                   <  2:1629930/150‑1 (MQ=255)

                                               ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa                                                                                                  <  1:1037663/150‑1 (MQ=255)

                                               ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa                                                                                                  <  1:127522/150‑1 (MQ=255)

                                                caccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGAt                                                                                                 <  1:1353646/150‑1 (MQ=255)

                                                   cgcgCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGcc                                                                                              <  2:592786/150‑1 (MQ=255)

                                                    gcgcTACTGCCGCCAGACAAATTCTTTTCTAATATGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGAACTCCCCCTTACCAAGGGTGCGCTCTACCAACAGAGCCATATCAGCAAGCTAAATTTTATGCCt                                                                                             >  2:1333182/1‑150 (MQ=37)

                                                     ggcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                            <  1:582659/149‑1 (MQ=255)

                                                       cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGc                                                                                          >  2:1703119/1‑150 (MQ=255)

                                                        tACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCa                                                                                         <  1:655156/150‑1 (MQ=255)

                                                         aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                        <  2:1358532/150‑1 (MQ=255)

                                                         aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                        <  1:836784/150‑1 (MQ=255)

                                                            gccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                     >  2:541219/1‑150 (MQ=255)

                                                            gccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                     >  1:143296/1‑150 (MQ=255)

                                                            gccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAACAAGGGAGCGCTCTACAAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                     >  1:1393797/1‑150 (MQ=255)

                                                             ccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTcc                                                                                    >  2:493143/1‑150 (MQ=255)

                                                               gccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTATTTATTTTTTTTGATAGCCAGAGTCGAACTGGTTAACTAAAATTTACCAATGGTGATCTCTACGATCTTAGCATAATCATTTCGCTATATTTGATGCCTGGCAGTTCCCt                                                                                  <  2:334470/150‑1 (MQ=255)

                                                                 cAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                <  2:720432/150‑1 (MQ=255)

                                                                 cAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                >  2:1657662/1‑150 (MQ=255)

                                                                 cAGACAAATTCCTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                >  1:1128221/1‑150 (MQ=255)

                                                                      aaaTTCTTTTCCAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCCACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACACTCg                                                                           >  1:844147/1‑150 (MQ=255)

                                                                        aTTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCCCGCTAAATTTGATGCCTGGCACTTCCCTACTCTCGCa                                                                         >  2:346706/1‑150 (MQ=255)

                                                                         ttCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCAt                                                                        >  1:557037/1‑150 (MQ=255)

                                                                           cTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCCCCCTTACCAAGGGTGCGATCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgg                                                                      >  2:828935/1‑150 (MQ=255)

                                                                               tCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgag                                                                  >  2:1534571/1‑150 (MQ=255)

                                                                                    tCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                             <  1:1476347/150‑1 (MQ=255)

                                                                                         cGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                        <  2:143296/150‑1 (MQ=255)

                                                                                             cTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTcccg                                                    >  2:658746/1‑149 (MQ=255)

                                                                                                tAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTCCCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTGCCCGACTCTCGCATGGGGAGACACCACACTACCATCg                                                 >  1:863490/1‑150 (MQ=255)

                                                                                                  aCCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgc                                               >  1:1128043/1‑150 (MQ=255)

                                                                                                  aCCTAAAAAGTGGTGCTGATACCAAGAGTCGAACTTGGGACATAAAAATTACCAAGGGTGAGCTCAACCAACATAGCCATATAAGAACGCTAAATTTGATGCCTGGCAGTTCCCTACTCACGCATGGGTAGACCCAAAACTACCATCGgc                                               >  2:1544002/1‑150 (MQ=21)

                                                                                                   ccTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcg                                              >  1:194543/1‑150 (MQ=255)

                                                                                                    cTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                             >  2:1250257/1‑150 (MQ=255)

                                                                                                     tAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCCCACTACCATCGGCGCt                                            >  2:1365557/1‑150 (MQ=255)

                                                                                                      aaaaaGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTCCCATCGGCGCTa                                           >  1:1422113/1‑150 (MQ=255)

                                                                                                          aGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGc                                       >  2:2150/1‑150 (MQ=255)

                                                                                                             ggtgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGCGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                    >  1:1094776/1‑150 (MQ=255)

                                                                                                             ggtgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCCTCGGCGCTACTGCttt                                    >  1:229293/1‑150 (MQ=255)

                                                                                                             ggtgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                    >  1:817837/1‑150 (MQ=255)

                                                                                                               tgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTc                                  <  2:1695972/150‑1 (MQ=255)

                                                                                                                gCTGATACCCAGACTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGATCCATATCAGCACGCTAAATTTGATGCCTTTAAGTTCTCTACTCTCGCATGGGGATACCCAACACTACCATCGGCGCTACGGCGTTTCa                                 <  2:814858/150‑1 (MQ=37)

                                                                                                                 cTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                >  1:1625556/1‑150 (MQ=255)

                                                                                                                  tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                               <  2:1128221/150‑1 (MQ=255)

                                                                                                                  tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                               >  1:614287/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                          >  1:415576/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGACATATCAGCACGAGAAATTTGATGCCTGGCAGTTCCCTCCTCTCGCATGGGGAGACCCCACCCAACCATCGGCGCTACGGGGGTTCACGTCTgg                          >  1:940752/1‑149 (MQ=25)

                                                                                                                           gagTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                      <  2:1393797/150‑1 (MQ=255)

                                                                                                                               cGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                  >  1:865049/1‑150 (MQ=255)

                                                                                                                                gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                 >  2:809628/1‑150 (MQ=255)

                                                                                                                                 aaCTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTGCCATCGGCGCTCCGGCGTTTCAGTTCTGAGTTCGGCATg                >  1:506567/1‑150 (MQ=255)

                                                                                                                                  aCTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg               >  2:795230/1‑150 (MQ=255)

                                                                                                                                  aCTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg               <  2:1416298/150‑1 (MQ=255)

                                                                                                                                    tGGGGACCTCTCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTTTTGTTTGGTTGTTCCGTTATCTCGGATTGGGAGACCTCTATATACTATATGCGTTAAGGAGTTTTTCTTATGTGTTCGGCATgggg             <  2:1670004/150‑1 (MQ=12)

                                                                                                                                    tGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgggg             <  1:909388/150‑1 (MQ=255)

                                                                                                                                       ggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGAAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa          >  1:283956/1‑150 (MQ=255)

                                                                                                                                        gACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg         >  1:1618764/1‑150 (MQ=255)

                                                                                                                                         agcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAgg        >  2:370321/5‑150 (MQ=255)

                                                                                                                                          gcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCCTATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt       >  1:1596059/4‑150 (MQ=18)

                                                                                                                                          gcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt       >  1:1651121/4‑150 (MQ=255)

                                                                                                                                          ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCCGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt       <  1:1682149/150‑1 (MQ=255)

                                                                                                                                          ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt       >  2:1704014/1‑150 (MQ=255)

                                                                                                                                          ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATAAGCACGCTAAATTTGATGCCTGGCAGTTCCCAACTCTCGCATGGGGAGACACCACACTACCATCGGCGCTAAGGCGTTTCACATCTGAGTTCGGCATGGGGTCAGGt       >  1:369867/1‑150 (MQ=21)

                                                                                                                                           cTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg      >  2:912717/1‑150 (MQ=255)

                                                                                                                                             cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCCCACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg    >  2:1104531/1‑150 (MQ=255)

                                                                                                                                             cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCCTCGGCGCTACGGCGTTTCCCTTCTGAGTTCGGCATGGGGTCAGGTggg    >  2:328355/1‑150 (MQ=255)

                                                                                                                                             cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg    >  2:1381108/1‑150 (MQ=255)

                                                                                                                                              aCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa   >  1:225423/1‑150 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt            >  1:51906/1‑140 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt            <  2:51906/140‑1 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac  <  1:493143/150‑1 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCGCGCAGGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac  <  1:346706/150‑1 (MQ=255)

GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGAC  >  NC_000913/3423472‑3423764

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