BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,836,277	A→G	5.8%	noncoding (56/95 nt)	selC →	tRNA‑Sec

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,836,277	0	A	G	5.8%	344.2 / 9.4	123	noncoding (56/95 nt)	selC	tRNA‑Sec
Reads supporting (aligned to +/- strand): ref base A (60/53); new base G (4/3); total (66/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-02

Only 100 of 123 total aligned reads displayed.

GCGTGTATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGTTTTGTGA  >  NC_000913/3836129‑3836425

 cGTGTATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGc                                                                                                                                                    >  2:1339629/1‑150 (MQ=255)

 cGTGTATTAGGCGGAAAAAACTGAGCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGc                                                                                                                                                    <  1:683032/150‑1 (MQ=255)

  gtgtATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCg                                                                                                                                                   <  2:1312613/150‑1 (MQ=255)

   tgtATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCgg                                                                                                                                                  >  1:640476/1‑150 (MQ=255)

   tgtATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCgg                                                                                                                                                  >  1:646703/1‑150 (MQ=255)

   tgtATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCgg                                                                                                                                                  <  2:1269130/150‑1 (MQ=255)

        tAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCg                                                                                                                                             >  1:513833/1‑150 (MQ=255)

        tAGGCGGAAAAAACTGATCTGGCGGATGTAGAAACTCAAGGAAGTACCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGACCGCTGGCGGcccca                                                                                                                                             <  1:1042908/150‑2 (MQ=255)

         aGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTCGCAATCGAGGCGCGGAAGATCGTCGTATCCGGTGAGGCGGCTTGACTTAAAATCCAGTTGGGACCGCTGGCGGccccaa                                                                                                                                            <  2:1190026/150‑3 (MQ=255)

          ggCGTAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCggg                                                                                                                                           <  2:530371/150‑1 (MQ=255)

          ggCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCCGCGGTCCCggg                                                                                                                                           >  2:456752/1‑150 (MQ=255)

          ggCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGACCGCTGGCGGccccaac                                                                                                                                           <  1:281408/150‑4 (MQ=255)

             ggAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATACTGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAg                                                                                                                                        <  2:1467110/150‑1 (MQ=255)

                aaaaaCTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGtt                                                                                                                                     >  2:294209/1‑150 (MQ=255)

                   aaCTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCACCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGa                                                                                                                                  >  2:1518451/1‑150 (MQ=255)

                    ccTGCTCTGGGTGATGTAGAAACTCAAGTAAGTATCTTTATTGCGCCCCGCCTCCATGTAGAAATCGTGGCGCGGAAGATCGTCGTCTCCTGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGAc                                                                                                                                 <  1:775614/149‑1 (MQ=255)

                     cTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGAcg                                                                                                                                >  1:606500/1‑149 (MQ=255)

                          cTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCtgt                                                                                                                           >  1:1101976/1‑150 (MQ=255)

                           tGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCtgtg                                                                                                                          >  1:1002961/1‑150 (MQ=255)

                            gggggATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGa                                                                                                                         <  2:758225/150‑1 (MQ=255)

                                   gTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGAt                                                                                                                        >  2:60343/1‑144 (MQ=255)

                                   gTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGAt                                                                                                                        <  1:60343/144‑1 (MQ=255)

                                   gTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCg                                                                                                                  <  2:931187/150‑1 (MQ=255)

                                    tAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGc                                                                                                                 >  1:769136/1‑150 (MQ=255)

                                      gAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCa                                                                                                               >  2:160959/1‑150 (MQ=255)

                                       aaaCTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCaa                                                                                                              <  2:1034705/150‑1 (MQ=255)

                                          cTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGCCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAAt                                                                                                           >  1:1675667/1‑150 (MQ=255)

                                           tCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATg                                                                                                          <  1:859166/150‑1 (MQ=255)

                                           tCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATg                                                                                                          <  1:471536/150‑1 (MQ=255)

                                            cAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGc                                                                                                         >  1:933707/1‑150 (MQ=255)

                                            cAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGc                                                                                                         >  2:1411765/1‑150 (MQ=255)

                                             aaGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGcc                                                                                                        <  1:1362103/150‑1 (MQ=255)

                                                gAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCctct                                                                                                     <  1:294209/150‑1 (MQ=255)

                                                  aGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTcc                                                                                                   <  2:1451331/150‑1 (MQ=255)

                                                  aGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTcc                                                                                                   >  1:1428647/1‑150 (MQ=255)

                                                      gCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGAc                                                                                               >  1:829744/1‑150 (MQ=255)

                                                      gCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGACAAAATGCCTCTCCTGAc                                                                                               >  2:786363/1‑150 (MQ=255)

                                                      gCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGGTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGAc                                                                                               >  2:779792/1‑150 (MQ=255)

                                                            aTGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTtctc                                                                                         >  2:394234/1‑150 (MQ=255)

                                                             ttttttccGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTcc                                                                                        <  1:96657/144‑1 (MQ=255)

                                                             gccgcCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTcc                                                                                        <  2:384563/148‑1 (MQ=255)

                                                                  cccGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAAccc                                                                                   >  2:609315/1‑150 (MQ=255)

                                                                    cGCCTCCCTGTAGCAATCGAGGCGGGGAAGATCGGCGGCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCAAACCCCt                                                                                 <  2:235143/150‑1 (MQ=255)

                                                                    cGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCt                                                                                 >  1:53719/1‑150 (MQ=255)

                                                                    cGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCt                                                                                 <  2:769136/150‑1 (MQ=255)

                                                                     gCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCtt                                                                                >  1:160170/1‑150 (MQ=255)

                                                                      ccTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTg                                                                               >  2:544853/1‑150 (MQ=255)

                                                                       cTCCATGTCGCAATCGAGGCGCGGAAGAGCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGt                                                                              <  1:1681473/150‑1 (MQ=255)

                                                                        tCCATGTAGCAATCGAGGCGCGGAAGATCGTCGGCTCCGGGGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTa                                                                             <  1:1604169/150‑1 (MQ=255)

                                                                         ccATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGTCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTACGACAAAATGCCTCTCCTGACCTTCTCCACACACTTGTaa                                                                            >  2:482624/1‑150 (MQ=255)

                                                                           aTGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATa                                                                          <  2:1277897/150‑1 (MQ=255)

                                                                            tGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATAc                                                                         <  1:68193/150‑1 (MQ=255)

                                                                              tAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGCCGCCTGTGATCTTCCGCCACAATGCCTCTCCTGACGTTCTCCAACCCCGTGTAATACTg                                                                       >  1:1147666/1‑150 (MQ=255)

                                                                              tAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTg                                                                       >  1:113145/1‑150 (MQ=255)

                                                                                gCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTc                                                                     >  2:1016810/1‑150 (MQ=255)

                                                                                gCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTc                                                                     <  2:513833/150‑1 (MQ=255)

                                                                                gCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTc                                                                     <  2:512307/150‑1 (MQ=255)

                                                                                   aTCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAAt                                                                  <  1:524229/150‑1 (MQ=255)

                                                                                          cgcgGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGata                                                           >  1:1106623/1‑150 (MQ=255)

                                                                                                 gATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTa                                                    <  1:1588190/150‑1 (MQ=255)

                                                                                                  aTCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCCACCCCTTGTAATACTGTCAATCCAGAGATTTTCTaa                                                   >  1:494013/1‑149 (MQ=255)

                                                                                                  aTCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTAc                                                   <  1:160959/150‑1 (MQ=255)

                                                                                                     gtcgtcTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGc                                                >  1:454123/1‑150 (MQ=255)

                                                                                                      gcggctctGGTGAGGCGGCTGGACTTCAATTCCAGTTGGGGCCGCCTGCGGTCCCGGGCAGGTTCGACTCCTGTGTTCTTCCTCAAAAATGCCTCTCCTGACTTTCTCCAACCCCTTGTATTACTGTCAATCCAGATATTTTCTACTGCt                                               <  1:240879/146‑1 (MQ=255)

                                                                                                       cgtcTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTc                                              <  2:160170/150‑1 (MQ=255)

                                                                                                         tctcCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGGAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATCCTGTCAAGCCAGATATTTTCTACTGCTaaa                                            >  2:888277/1‑150 (MQ=255)

                                                                                                         tctcCGGTGAGGCGGATGGACTTCAAATCCAGTTGGTGCCGCCAGCGGTCCCGGGCAGGTTCGACTCATGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTGAATACAGAAATTTTCTAATGATCaa                                            >  2:894285/1‑150 (MQ=255)

                                                                                                          ctcCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCACCCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAat                                           >  1:1089408/1‑150 (MQ=255)

                                                                                                           tcCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAata                                          >  1:1510865/1‑150 (MQ=255)

                                                                                                            ccGGTGAGGCTTATGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAatat                                         <  2:646703/150‑1 (MQ=255)

                                                                                                            ccGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAatat                                         <  2:640476/150‑1 (MQ=255)

                                                                                                            ccGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAatat                                         >  2:353641/1‑150 (MQ=255)

                                                                                                              ggTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTc                                       <  1:544853/150‑1 (MQ=255)

                                                                                                                tGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCaa                                     >  2:1658978/1‑150 (MQ=255)

                                                                                                                   ggcggcTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCg                                  <  1:1408497/150‑1 (MQ=255)

                                                                                                                    gcggcTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCgg                                 <  1:1642107/150‑1 (MQ=255)

                                                                                                                     cggcTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTc                                       >  2:644873/1‑143 (MQ=255)

                                                                                                                      ggcTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCt                               <  2:1572593/150‑1 (MQ=255)

                                                                                                                      ggcTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCt                               >  1:1344568/1‑150 (MQ=255)

                                                                                                                         tGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCa                            <  2:560798/150‑1 (MQ=255)

                                                                                                                          ggACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAg                           <  2:829744/150‑1 (MQ=255)

                                                                                                                          ggACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAg                           >  1:1340136/1‑150 (MQ=255)

                                                                                                                           gACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAgg                          <  2:53719/150‑1 (MQ=255)

                                                                                                                            aCTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGa                         >  1:1249478/1‑150 (MQ=255)

                                                                                                                              ttCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATAATCAATCGGCTTCAGGAAt                       >  2:1116411/1‑150 (MQ=255)

                                                                                                                              ttCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGGCACTCCAGATATTTTCTACTGCTCAATAATCAATCGGCTTCAGGAAt                       >  1:283248/1‑150 (MQ=255)

                                                                                                                               tCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATa                      >  2:1704294/1‑150 (MQ=255)

                                                                                                                                cAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCGTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATAc                     >  1:1708893/1‑150 (MQ=255)

                                                                                                                                 aaaTCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACGGTCAATCCAGATATTTTCTACTGCTCAATATTCATTCGGCTACAGGAATAcc                    >  1:1591234/1‑150 (MQ=255)

                                                                                                                                  aaTCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCt                   >  1:521982/1‑150 (MQ=255)

                                                                                                                                    tCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTgg                 <  2:1101976/150‑1 (MQ=255)

                                                                                                                                        cccgggaCCGCTGGCGGCCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCtt             >  1:1360963/4‑149 (MQ=255)

                                                                                                                                        cccgggaCCGCTGGCGGCCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTa             >  1:1361235/4‑150 (MQ=255)

                                                                                                                                        cccgggaCCGCTGGCGGCCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTa             >  1:525834/4‑150 (MQ=255)

                                                                                                                                         ccgggaCCGCTGGCGGCCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCACCCCCTTGTAATACTGTCAATCCAGATATTTTCTACGGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTaa            >  2:1297931/3‑150 (MQ=255)

                                                                                                                                           gaacctgccAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAgg          >  1:1180267/7‑150 (MQ=255)

                                                                                                                                             ggCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGtt        >  1:574288/1‑150 (MQ=255)

                                                                                                                                                 gccAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGTTTtgt    <  1:1411765/150‑1 (MQ=255)

                                                                                                                                                  ccAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGTTTtgtg   <  2:1147666/150‑1 (MQ=255)

                                                                                                                                                   cAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGTTTggtgc  >  2:1326959/1‑149 (MQ=255)

GCGTGTATTAGGCGGAAAAAACTGATCTGGGGGATGTAGAAACTCAAGGAAGTAGCTATAATGCGCCCCGCCTCCATGTAGCAATCGAGGCGCGGAAGATCGTCGTCTCCGGTGAGGCGGCTGGACTTCAAATCCAGTTGGGGCCGCCAGCGGTCCCGGGCAGGTTCGACTCCTGTGATCTTCCGCCAAAATGCCTCTCCTGACGTTCTCCAACCCCTTGTAATACTGTCAATCCAGATATTTTCTACTGCTCAATATTCAATCGGCTTCAGGAATACCTGGTCTAAGGTTTTGTGA  >  NC_000913/3836129‑3836425

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