D3A :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,598,172	A→G	S110S (TCT→TCC)	SAUSA300_RS07880 ←	pseudouridine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,598,172	0	A	G	99.5%	533.7 / ‑5.6	202	S110S (TCT→TCC)	SAUSA300_RS07880	pseudouridine synthase B
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (2/199); total (2/200)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

Only 100 of 202 total aligned reads displayed.

AACCTTTTAATTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCAGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAGATAC  >  NC_007793/1598075‑1598267

   cTTTTAATTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGt                                                                                            <  1:1425040/100‑1 (MQ=255)

        aaTTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCAt                                                                                       <  1:3907019/100‑1 (MQ=255)

        aaTTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCAt                                                                                       <  1:4430462/100‑1 (MQ=255)

        aaTTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCAt                                                                                       <  1:397967/100‑1 (MQ=255)

        aaTTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCAt                                                                                       <  1:1286262/100‑1 (MQ=255)

        aaTTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATGTGTGAGTAATAACAATCCGGAAGTATCAt                                                                                       <  1:6316289/100‑1 (MQ=255)

             cGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCt                                                                                  <  1:3661311/100‑1 (MQ=255)

             cGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCt                                                                                  <  1:2234247/100‑1 (MQ=255)

             cGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCt                                                                                  <  1:136981/100‑1 (MQ=255)

             cGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCt                                                                                  <  1:1897036/100‑1 (MQ=255)

              gCATCATATTTATTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCGCCATCATTTGAGAGTAATAACAATCCGGAAGTATCATAGTCTa                                                                                 <  1:6605660/100‑1 (MQ=255)

              gCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCagatcgga                                                                                 >  1:3119864/1‑93 (MQ=255)

              gCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTa                                                                                 <  1:3038615/100‑1 (MQ=255)

              gCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTa                                                                                 <  1:703927/100‑1 (MQ=255)

                     atTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGAcc                                                                          <  1:2225440/100‑1 (MQ=255)

                      ttttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCa                                                                         <  1:664488/100‑1 (MQ=255)

                      ttttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCa                                                                         <  1:3838790/100‑1 (MQ=255)

                        ttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCaag                                                                       <  1:5392636/100‑2 (MQ=255)

                        ttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAAc                                                                       <  1:1533927/100‑1 (MQ=255)

                        ttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAAc                                                                       <  1:2094041/100‑1 (MQ=255)

                        ttttCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCAGTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAAc                                                                       <  1:2584883/100‑1 (MQ=255)

                           tCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAgg                                                                    <  1:6335967/100‑1 (MQ=255)

                           tCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAgg                                                                    <  1:1720699/100‑1 (MQ=255)

                           tCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAgg                                                                    <  1:1797232/100‑1 (MQ=255)

                           tCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAgg                                                                    <  1:2974150/100‑1 (MQ=255)

                           tCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAgg                                                                    <  1:2920636/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACATTCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:2382185/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGac                                                                  <  1:5238914/100‑2 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:3733836/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:5941316/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:80790/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:4106629/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGAt                                                                  <  1:52255/100‑1 (MQ=255)

                             ttAATTTGATATCTTGGATGTGTCATTAAATGAGTAAACTCTCCATCATGTGTGAGTAATAACAATCCGGAAGTATCATAGTCGAAACGACCAACAGGAt                                                                  <  1:5264544/100‑1 (MQ=255)

                              tAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATa                                                                 <  1:3718704/100‑1 (MQ=255)

                               aaTTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATaa                                                                <  1:6091839/100‑1 (MQ=255)

                                    gATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACg                                                           <  1:6494348/100‑1 (MQ=255)

                                       atCTTGGATGTGTCATTAAATTAGTAAATTCACCATGATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTgg                                                        <  1:5428281/100‑2 (MQ=255)

                                       atCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTgg                                                        <  1:5159750/100‑2 (MQ=255)

                                       atCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGc                                                        <  1:4593663/100‑1 (MQ=255)

                                       atCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGc                                                        <  1:494459/100‑1 (MQ=255)

                                        tCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATCCGTGCt                                                       <  1:6037959/100‑1 (MQ=255)

                                         cTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCtt                                                      <  1:1195588/100‑1 (MQ=255)

                                          ttGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTc                                                     <  1:4275005/100‑1 (MQ=255)

                                           tGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCt                                                    <  1:3323605/100‑1 (MQ=255)

                                           tGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCt                                                    <  1:4143798/100‑1 (MQ=255)

                                           tGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCt                                                    <  1:2672451/100‑1 (MQ=255)

                                            ggATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCta                                                   <  1:6014028/100‑1 (MQ=255)

                                              aTGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCtata                                                 <  1:555474/100‑1 (MQ=255)

                                                               tAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGt                                <  1:2371007/100‑1 (MQ=255)

                                                               tAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGt                                <  1:2464072/100‑1 (MQ=255)

                                                               tAAATTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGt                                <  1:3463624/100‑1 (MQ=255)

                                                                aaaTTCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTa                               <  1:4552247/100‑1 (MQ=255)

                                                                   ttCGCCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAAcc                            <  1:5739266/100‑1 (MQ=255)

                                                                   ttCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAACCGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAAcc                            <  1:2011600/100‑1 (MQ=255)

                                                                   ttCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAAcc                            <  1:2416176/100‑1 (MQ=255)

                                                                   ttCACCATCATTTGTGAGTAATAACAATCCGGAAGTACCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTAGATCTTTAAAGTAATCTGTAAcc                            <  1:5708302/100‑1 (MQ=255)

                                                                   ttCACCATCATTTGCGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAAcc                            <  1:4838503/100‑1 (MQ=255)

                                                                    tCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAAccc                           <  1:1233816/100‑2 (MQ=255)

                                                                    tCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCa                           <  1:3259861/100‑1 (MQ=255)

                                                                    tCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCa                           <  1:4751289/100‑1 (MQ=255)

                                                                    tCACCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCCTATATATCTTTAAAGTAATCTGTAACCa                           <  1:645629/100‑1 (MQ=255)

                                                                      gccATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACt                         >  1:3058580/2‑100 (MQ=255)

                                                                      aCCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACt                         <  1:5207143/100‑1 (MQ=255)

                                                                      aCCATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACt                         <  1:1053531/100‑1 (MQ=255)

                                                                       ccATCATTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTg                        <  1:283782/100‑1 (MQ=255)

                                                                        catcatTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGt                       <  1:5532395/100‑1 (MQ=255)

                                                                        catcatTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGt                       <  1:5354235/100‑1 (MQ=255)

                                                                         atcatTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCAcgtgc                      <  1:548263/100‑5 (MQ=255)

                                                                         atcatTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGtc                      <  1:5058762/100‑2 (MQ=255)

                                                                           catTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTAca                    <  1:1582220/100‑2 (MQ=255)

                                                                           catTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACg                    <  1:3750737/100‑1 (MQ=255)

                                                                            atTTGTGAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGt                   <  1:5572522/100‑1 (MQ=255)

                                                                                gtgAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCctt               <  1:296294/100‑2 (MQ=255)

                                                                                 tgAGTAATAACAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTCCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCctct              <  1:3682338/100‑1 (MQ=255)

                                                                                      aataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcat         <  1:1023375/100‑1 (MQ=255)

                                                                                      aataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcat         <  1:845980/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTGTCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:6128401/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGTCCAACAGGATAAATTCATGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:5498401/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:1092745/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:1326277/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:1601699/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:1672174/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:895200/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:1732040/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:2073586/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:2708312/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:5793191/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:5611239/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:5452399/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:3315397/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:5300674/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:4309596/100‑1 (MQ=255)

                                                                                       ataaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatc        <  1:4410399/100‑1 (MQ=255)

                                                                                         aaCAATCCGGAAGTGTCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAg      <  1:1169903/100‑1 (MQ=255)

                                                                                         aaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAg      <  1:2331888/100‑1 (MQ=255)

                                                                                         aaCAATCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAg      <  1:3624893/100‑1 (MQ=255)

                                                                                         aaCAATCCGGAAGTATCATAGTCTAAAAGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTACCATCAg      <  1:3222906/100‑1 (MQ=255)

                                                                                             aTCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTatcatcctctgt  <  1:6187490/100‑7 (MQ=255)

                                                                                             aTCCGGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAGATac  <  1:1415283/100‑1 (MQ=255)

AACCTTTTAATTTCGCAACATATTTTTTCTTAATTTGATATCTTGGATGTGTCATTAAATTAGTAAATTCACCATCATTTGTGAGTAATAACAATCCAGAAGTATCATAGTCTAAACGACCAACAGGATAAATACGTGCTTCTATATCTTTAAAGTAATCTGTAACCACTGTACGTCCTCTATCATCAGATAC  >  NC_007793/1598075‑1598267

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