BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	432,925	C→T	P83S (CCG→TCG)	ribD →	pyrimidine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	432,925	0	C	T	99.2%	458.9 / ‑5.3	133	P83S (CCG→TCG)	ribD	pyrimidine deaminase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (69/63); total (69/64)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 139 total aligned reads displayed.

GTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTT  >  NC_000913/432835‑433014

gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCg                                                                                         >  2:597937/1‑93 (MQ=255)

 tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGt                                                                                        <  1:4047568/93‑1 (MQ=255)

     cgcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCtg                                                                                    <  1:2604708/93‑1 (MQ=255)

      gcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCtgt                                                                                   >  1:3292903/1‑93 (MQ=255)

       cgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCtgtg                                                                                  <  2:259040/93‑1 (MQ=255)

         ttGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGAc                                                                                >  1:130316/1‑93 (MQ=255)

         ttGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGAc                                                                                >  2:1775550/1‑93 (MQ=255)

            cGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCa                                                                             <  1:636942/93‑1 (MQ=255)

            cGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCa                                                                             <  1:4449123/93‑1 (MQ=255)

             gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCAc                                                                            <  2:3288582/93‑1 (MQ=255)

             gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCAc                                                                            >  2:1874782/1‑93 (MQ=255)

             gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACCCAc                                                                            >  1:3355598/1‑93 (MQ=255)

                tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCa                                                                         >  1:676290/1‑93 (MQ=255)

                tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCa                                                                         >  2:371547/1‑93 (MQ=255)

                tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCa                                                                         <  1:1493798/93‑1 (MQ=255)

                 ggCGGGTGAAAAAGCCAAAGGTGCCACCGCCTATGTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCAt                                                                        <  2:3414605/93‑1 (MQ=255)

                 ggCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCCCTCAt                                                                        >  2:967160/1‑93 (MQ=255)

                  gCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATc                                                                       <  2:2826739/93‑1 (MQ=255)

                  gCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATc                                                                       <  1:715935/93‑1 (MQ=255)

                    gagTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgc                                                                     <  2:2749678/91‑1 (MQ=255)

                     ggTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgcc                                                                    <  2:4175985/93‑1 (MQ=255)

                     ggTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgcc                                                                    >  1:2532370/1‑93 (MQ=255)

                       tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgccgc                                                                  <  2:4552573/93‑1 (MQ=255)

                       tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgccgc                                                                  >  1:4141947/1‑93 (MQ=255)

                        gAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGTCGCACTCATCGCCGCt                                                                 <  1:1693193/93‑1 (MQ=255)

                        gAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCt                                                                 >  1:2260744/1‑93 (MQ=255)

                           aaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGc                                                              <  2:3900426/93‑1 (MQ=255)

                           aaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGc                                                              >  2:4736338/1‑93 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTTGCCATCATGGTCGGACGCCATTGTGCTGTGACGCACTCATCGCCGCCGGCGTAgc                                                         >  2:4225302/1‑93 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         >  2:4189598/1‑93 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         >  2:1273521/1‑93 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         <  2:3965698/93‑1 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         >  2:402314/1‑93 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         <  2:622001/93‑1 (MQ=255)

                                cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                         <  2:1177599/93‑1 (MQ=255)

                                 aaaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcg                                                        <  1:244557/93‑1 (MQ=255)

                                  aaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgc                                                       >  1:2626157/1‑93 (MQ=255)

                                  aaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgc                                                       <  2:4701729/93‑1 (MQ=255)

                                    ggTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcgc                                                     >  2:1251669/1‑93 (MQ=255)

                                       gCGACCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGATCGTACGCCCTCGCGCTGTGACGCACTCCTCGACGCTCGTCTAGCGCGAGTg                                                  >  1:191344/1‑93 (MQ=255)

                                        cGACCGCCTATGTCACATTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTgg                                                 <  2:3049612/93‑1 (MQ=255)

                                         gACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGt                                                <  2:353925/93‑1 (MQ=255)

                                          aCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGtt                                               >  2:231486/1‑93 (MQ=255)

                                             gCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGcc                                            >  2:3491158/1‑93 (MQ=255)

                                             gCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGcc                                            <  2:3375252/93‑1 (MQ=255)

                                             gCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGcc                                            >  1:3652766/1‑93 (MQ=255)

                                              ccTATGTCACACTCGACCCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                           <  2:964948/93‑1 (MQ=255)

                                              ccTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTGCGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGGGTGGTTGCCt                                           <  2:2575441/93‑1 (MQ=255)

                                              ccTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                           >  2:459533/1‑93 (MQ=255)

                                                tATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                         >  1:875427/1‑93 (MQ=255)

                                                tATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                         >  2:40000/1‑93 (MQ=255)

                                                 aTGTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGa                                        <  1:3329031/93‑1 (MQ=255)

                                                 aTGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGa                                        <  2:3020565/93‑1 (MQ=255)

                                                 aTGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGa                                        >  1:4466095/1‑93 (MQ=255)

                                                 aTGTCACACTCGAACCCTGTAGCCATCAGGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGAGTAGCACGCGTGGTTGCCCCGa                                        <  2:2475145/93‑1 (MQ=255)

                                                  tGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAt                                       <  2:2382939/93‑1 (MQ=255)

                                                   gTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACCCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATg                                      <  1:1138337/93‑1 (MQ=255)

                                                   gTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATg                                      <  1:4036504/93‑1 (MQ=255)

                                                      acacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCaa                                   <  1:2862835/93‑1 (MQ=255)

                                                      acacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCaa                                   <  2:665496/93‑1 (MQ=255)

                                                      acacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCaa                                   <  1:2862836/93‑1 (MQ=255)

                                                       cacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAt                                       <  2:259922/88‑1 (MQ=255)

                                                       cacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                  >  2:1939387/1‑93 (MQ=255)

                                                       cacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                  >  2:1921317/1‑93 (MQ=255)

                                                         cTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGAt                                <  1:750446/93‑1 (MQ=255)

                                                          tCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATc                               <  2:3960829/93‑1 (MQ=255)

                                                          tCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATc                               >  1:578818/1‑93 (MQ=255)

                                                          tCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATc                               >  1:2553496/1‑93 (MQ=255)

                                                           cGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATcc                              <  1:3423190/93‑1 (MQ=255)

                                                             aaCCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTa                            >  1:2802856/1‑93 (MQ=255)

                                                              aCCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTaa                           <  1:3123696/93‑1 (MQ=255)

                                                               cccTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAc                          >  2:1328079/1‑93 (MQ=255)

                                                                ccTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAcc                         <  2:1878126/93‑1 (MQ=255)

                                                                  tGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCg                       >  1:3767187/1‑93 (MQ=255)

                                                                  tGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCg                       >  1:1851375/1‑93 (MQ=255)

                                                                    tAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCa                     <  2:4530738/93‑1 (MQ=255)

                                                                     aGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAg                    <  2:4606194/93‑1 (MQ=255)

                                                                      gCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAgg                   <  1:2390891/93‑1 (MQ=255)

                                                                      gCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAgg                   >  2:2206141/1‑93 (MQ=255)

                                                                      gCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAgg                   >  1:3032167/1‑93 (MQ=255)

                                                                        catcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTc                 >  2:2997212/1‑93 (MQ=255)

                                                                        catcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTc                 <  1:4283263/93‑1 (MQ=255)

                                                                         atcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCg                >  2:2463200/1‑93 (MQ=255)

                                                                          tcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc               >  2:2043374/1‑93 (MQ=255)

                                                                          tcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc               >  2:3796867/1‑93 (MQ=255)

                                                                          tcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc               <  1:4684711/93‑1 (MQ=255)

                                                                             tGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTgg            <  1:3349635/93‑1 (MQ=255)

                                                                              ggTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTcgg           <  1:3104252/93‑1 (MQ=255)

                                                                               gTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc          >  2:4374352/1‑93 (MQ=255)

                                                                                tCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCg         <  1:3808718/93‑1 (MQ=255)

                                                                                 cGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTGGATGCAAGATCCTAACCCGCAGGTCGCTGGGCgg        >  1:871671/1‑92 (MQ=255)

                                                                                 cGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCgg        >  2:3462352/1‑92 (MQ=255)

                                                                                 cGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCgg        >  1:871808/1‑92 (MQ=255)

                                                                                 cGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt        <  1:4397512/93‑1 (MQ=255)

                                                                                 cGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt        >  2:2970455/1‑93 (MQ=255)

                                                                                  gTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTg       >  1:2480289/1‑93 (MQ=255)

                                                                                   tACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTgg      >  1:1769076/1‑93 (MQ=255)

                                                                                     cGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAc    >  2:3621288/1‑93 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGGGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  <  1:3642734/93‑1 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  <  1:3642602/93‑1 (MQ=255)

GTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTT  >  NC_000913/432835‑433014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`