BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,820,243:1	+T	5.6%	coding (928/1683 nt)	ligB ←	DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,820,243	1	.	T	5.6%	429.8 / 20.1	107	coding (928/1683 nt)	ligB	DNA ligase B
Reads supporting (aligned to +/- strand): ref base . (56/45); new base T (4/2); total (60/47)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 107 total aligned reads displayed.

ATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820097‑3820391

atCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGc                                                                                                                                                   <  2:748767/150‑1 (MQ=255)

atCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGc                                                                                                                                                   <  1:772707/150‑1 (MQ=255)

 tCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCt                                                                                                                                                  <  1:1520102/150‑1 (MQ=255)

   ccACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTAc                                                                                                                                                >  2:1412537/1‑150 (MQ=255)

    cACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACa                                                                                                                                               >  2:97209/1‑150 (MQ=255)

        ccTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGtt                                                                                                                                           <  2:745040/150‑1 (MQ=255)

          tGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGa                                                                                                                                         <  1:689405/150‑1 (MQ=255)

               gcgcCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccgaag                                                                                                                                    <  1:1188817/150‑7 (MQ=255)

                  ccTGACGTAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCa                                                                                                                                 <  2:572463/150‑1 (MQ=255)

                   cTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAg                                                                                                                                <  2:950831/150‑1 (MQ=255)

                       cGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGcca                                                                                                                            >  2:328748/1‑150 (MQ=255)

                          aaCCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACGTG‑AGCTACAGGTTGATATTTCCAGGccacca                                                                                                                         >  1:393498/1‑150 (MQ=255)

                            ccAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgc                                                                                                                                         <  1:920487/132‑2 (MQ=255)

                            ccAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgc                                                                                                                                         >  2:920487/1‑131 (MQ=255)

                             cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCCCCAGcc                                                                                                                      >  2:1503347/1‑150 (MQ=255)

                             cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                      >  1:1308612/1‑150 (MQ=255)

                             cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                      >  2:629098/1‑150 (MQ=255)

                              aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACAGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                     >  2:842361/1‑150 (MQ=255)

                                tatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAct                                                                                                                   >  2:973471/1‑150 (MQ=255)

                                tatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAct                                                                                                                   >  1:144414/1‑150 (MQ=255)

                                tatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAct                                                                                                                   >  1:150137/1‑150 (MQ=255)

                                        cGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTgg                                                                                                           >  1:21073/1‑150 (MQ=255)

                                          cTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGcc                                                                                                         >  1:1618910/1‑150 (MQ=255)

                                            ggACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCg                                                                                                       >  1:22396/1‑150 (MQ=255)

                                             gACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCgg                                                                                                      <  1:1444721/150‑1 (MQ=255)

                                                     tATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAAt                                                                                              <  2:104015/150‑1 (MQ=255)

                                                       tcatcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATgg                                                                                            >  2:48195/1‑150 (MQ=255)

                                                         atcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATg                                                                                             <  2:1092833/147‑1 (MQ=255)

                                                         atcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATg                                                                                             >  1:1092833/1‑147 (MQ=255)

                                                         atcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCg                                                                                          >  2:413807/1‑150 (MQ=255)

                                                         atcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCg                                                                                          >  2:407843/1‑150 (MQ=255)

                                                          tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         <  1:798395/150‑1 (MQ=255)

                                                          tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         <  2:457866/150‑1 (MQ=255)

                                                           ccAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCggg                                                                                        >  1:926928/1‑150 (MQ=255)

                                                              gCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGAtt                                                                                     >  2:38688/1‑150 (MQ=255)

                                                                 tGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTg                                                                                  <  1:560362/150‑1 (MQ=255)

                                                                 tGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTg                                                                                  <  2:1382137/150‑1 (MQ=255)

                                                                  gACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTgg                                                                                 >  2:1168698/1‑150 (MQ=255)

                                                                   aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                <  1:1261792/150‑1 (MQ=255)

                                                                   aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                <  2:816073/150‑1 (MQ=255)

                                                                    cAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTcc                                                                                                                                  <  1:602670/99‑1 (MQ=255)

                                                                      ggTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGctct                                                                             >  1:1426487/1‑150 (MQ=255)

                                                                      ggTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGctct                                                                             >  1:1319905/1‑150 (MQ=255)

                                                                          cGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                         <  2:926928/150‑1 (MQ=255)

                                                                           gAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgc                                                                        >  1:1200984/1‑150 (MQ=255)

                                                                            aGTGACGCAACCACCGATATTTTACCGCTCTTAACCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑ATCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGTTAGCCAATGGCGGGATTCTGGCTCTTTcgcc                                                                       <  2:169186/150‑1 (MQ=255)

                                                                               gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCCCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                    <  1:477633/150‑1 (MQ=255)

                                                                               gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                    >  2:1285289/1‑150 (MQ=255)

                                                                               gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                    >  2:943387/1‑150 (MQ=255)

                                                                               gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                    >  2:1408108/1‑150 (MQ=255)

                                                                                 cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                  <  2:1416339/150‑1 (MQ=255)

                                                                                   cAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTa                                                                <  1:97209/150‑1 (MQ=255)

                                                                                        accGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACt                                                           >  2:39749/1‑150 (MQ=255)

                                                                                           gATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTAcc                                                        >  1:694570/1‑149 (MQ=255)

                                                                                               ttttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCAt                                                    >  2:653625/1‑150 (MQ=255)

                                                                                               ttttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCAt                                                    >  2:452922/1‑150 (MQ=255)

                                                                                                tttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATc                                                   <  2:254359/150‑1 (MQ=255)

                                                                                                tttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATc                                                   <  2:612059/150‑1 (MQ=255)

                                                                                                  tACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCgg                                                 <  1:260767/150‑1 (MQ=255)

                                                                                                  tACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCCGGCTCTTTCGCCGCTCGTACAACTACGCCATCgg                                                 <  1:930001/150‑1 (MQ=255)

                                                                                                    ccGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTg                                               <  1:48195/150‑1 (MQ=255)

                                                                                                     cGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGa                                              <  2:150137/150‑1 (MQ=255)

                                                                                                     cGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGa                                              >  1:597448/1‑150 (MQ=255)

                                                                                                      gCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGAc                                             <  1:1503347/150‑1 (MQ=255)

                                                                                                         ctTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGaa                                          >  2:501734/1‑150 (MQ=255)

                                                                                                           tACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                        >  1:367340/1‑150 (MQ=255)

                                                                                                           tACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                        <  2:935886/150‑1 (MQ=255)

                                                                                                            aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                       >  2:365571/1‑150 (MQ=255)

                                                                                                            aCCCACCGCAAACTGAATTGCCTTAACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGACACTCTGCCTGGCCCGGTATCCAATGGCGGGATTCTGGCTCTTTCGCCGCTAGTACAACTACGCCATCGGTGACGAAggg                                       >  2:124033/1‑150 (MQ=255)

                                                                                                               cACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTaa                                    >  2:1576664/1‑150 (MQ=255)

                                                                                                                 ccGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAtt                                  >  2:1288945/1‑150 (MQ=255)

                                                                                                                   gCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                >  2:1551060/1‑150 (MQ=255)

                                                                                                                   gCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                <  2:1491116/150‑1 (MQ=255)

                                                                                                                     aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                              <  1:1457182/150‑1 (MQ=255)

                                                                                                                     aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                              >  2:1479130/1‑150 (MQ=255)

                                                                                                                        cTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCttt                           <  2:22396/150‑1 (MQ=255)

                                                                                                                         tGAATTGCCTTCACTTCGGCAACGTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTc                          <  2:393498/150‑1 (MQ=255)

                                                                                                                          gAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         <  2:147246/150‑1 (MQ=255)

                                                                                                                                 cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  1:1131899/1‑150 (MQ=255)

                                                                                                                                 cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  <  1:783563/150‑1 (MQ=255)

                                                                                                                                 cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  2:576568/1‑150 (MQ=255)

                                                                                                                                 cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  2:1351621/1‑150 (MQ=255)

                                                                                                                                 cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  2:1553420/1‑150 (MQ=255)

                                                                                                                                   tCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAt                >  1:650290/1‑150 (MQ=255)

                                                                                                                                   tCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAt                >  1:222087/1‑150 (MQ=255)

                                                                                                                                     aCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg              <  1:452922/150‑1 (MQ=255)

                                                                                                                                        tCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAATGGTAATTCCGCTTTCCACCACTCATTGCGt           <  1:842361/150‑1 (MQ=255)

                                                                                                                                        tCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt           >  2:1142585/1‑150 (MQ=255)

                                                                                                                                         cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          >  1:954802/1‑150 (MQ=255)

                                                                                                                                         cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          <  1:653625/150‑1 (MQ=255)

                                                                                                                                         cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          <  1:649015/150‑1 (MQ=255)

                                                                                                                                         cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          <  1:413807/150‑1 (MQ=255)

                                                                                                                                         cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          <  1:407843/150‑1 (MQ=255)

                                                                                                                                          atcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAc         >  1:177818/3‑150 (MQ=255)

                                                                                                                                             aaacTG‑AGCTACAGGTTGATATTTCCAGTCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGCGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcgcg      >  2:1540005/4‑150 (MQ=255)

                                                                                                                                              aCCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcgcgc     >  1:260100/1‑150 (MQ=255)

                                                                                                                                                cTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACgcgcgc   <  2:222087/150‑1 (MQ=255)

                                                                                                                                                cTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACgcgcgc   >  2:1468487/1‑150 (MQ=255)

                                                                                                                                                cTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACgcgcgc   <  2:489939/150‑1 (MQ=255)

                                                                                                                                                 tG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCa  >  1:1426904/1‑150 (MQ=255)

ATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820097‑3820391

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