BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,820,248	Δ1 bp	5.7%	coding (923/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,248	0	A	.	5.7%	353.5 / 12.5	106	coding (923/1683 nt)	ligB	DNA ligase B
Reads supporting (aligned to +/- strand): ref base A (57/43); new base . (4/2); total (61/45)
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.30e-01

Only 100 of 106 total aligned reads displayed.

CCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTT  >  NC_000913/3820100‑3820395

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:150137/1‑150 (MQ=255)

                             tatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAct                                                                                                                       >  1:144414/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                                                                                                 >  1:1092833/1‑147 (MQ=255)

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

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

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

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

                                                       tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                             <  1:798395/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                                                                                                                                      >  2:602670/1‑99 (MQ=255)

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

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

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

                                                                   ggTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGctct                                                                                 >  1:489939/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‑AGCTACAGGTTGATATTTCCAGGCCCCCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                        >  1:927026/1‑150 (MQ=255)

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

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

                                                                            gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCt                                                                        >  2:943387/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:452922/1‑150 (MQ=255)

                                                                                            ttttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCAt                                                        >  2:653625/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                                                  >  1:597448/1‑150 (MQ=255)

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

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

                                                                                                   gCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGAc                                                 >  1:1577113/1‑150 (MQ=255)

                                                                                                   gCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGAc                                                 <  1:1503347/150‑1 (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)

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

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

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

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

                                                                                                                  aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                                  <  1:1457182/150‑1 (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                      >  2:576568/1‑150 (MQ=255)

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

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

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

                                                                                                                              cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                      <  1:783563/150‑1 (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)

                                                                                                                                      tatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa              >  2:296768/4‑150 (MQ=255)

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

                                                                                                                                      cGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa              <  1:413807/150‑1 (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)

                                                                                                                                          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:489939/150‑1 (MQ=255)

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

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

                                                                                                                                                  gCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACt   <  1:1551060/150‑1 (MQ=255)

                                                                                                                                                   cTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACtt  >  2:1591692/1‑150 (MQ=255)

CCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTT  >  NC_000913/3820100‑3820395

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