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.6%	860.4 / ‑5.4	232	P83S (CCG→TCG)	ribD	pyrimidine deaminase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (123/108); total (123/109)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 242 total aligned reads displayed.

GTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTT  >  NC_000913/432835‑433014

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

gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCg                                                                                         <  1:3300786/93‑1 (MQ=255)

gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCg                                                                                         <  1:3955950/93‑1 (MQ=255)

gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCg                                                                                         <  1:1782034/93‑1 (MQ=255)

gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCg                                                                                         <  1:2018498/93‑1 (MQ=255)

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

  acacGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTg                                                                                       <  1:6587564/93‑1 (MQ=255)

   cacGCGTTGCGTATGGCGGGCGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATCGTCGTACGCCATCGTGc                                                                                      >  1:3736772/1‑93 (MQ=255)

    acGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCt                                                                                     >  1:1138002/1‑93 (MQ=255)

      gcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCtgt                                                                                   <  1:5363198/93‑1 (MQ=255)

      gcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCtgt                                                                                   <  1:465711/93‑1 (MQ=255)

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

       cgTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGAACGCCTATGTGAAACTCGAAACCCGTAGCCATCATGGTCGTACGCCATCGTGCtgtg                                                                                  >  2:1445502/1‑93 (MQ=255)

        gTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGa                                                                                 <  1:461124/93‑1 (MQ=255)

        gTTGCGTATGGCGGGTGAAAAAGCCAAAGGCGCGACCGCCTATGTCACACTCGGACCCTGTCGCCATCATGGGCTTACGCCATTGTGCtgttc                                                                                 >  2:2699333/1‑91 (MQ=255)

        gTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGa                                                                                 <  2:4316336/93‑1 (MQ=255)

          tGCGTATGGCGGGTGAAAAAGCCAAAGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACg                                                                               >  2:3388314/1‑93 (MQ=255)

             gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGCTGCAc                                                                            >  2:5645465/1‑93 (MQ=255)

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

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

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

                tGGCGGGTGAAAAAGCCAAAGGTCCGACCACCTATGTCACACTCGAACGCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCa                                                                         <  1:5697328/93‑1 (MQ=255)

                 ggCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCAt                                                                        >  2:5034243/1‑93 (MQ=255)

                 ggCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCAt                                                                        >  1:1892204/1‑93 (MQ=255)

                   cGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcg                                                                      <  1:1592729/93‑1 (MQ=255)

                    gggTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgc                                                                     >  1:6332970/1‑93 (MQ=255)

                    gggTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgc                                                                     <  2:6527913/93‑1 (MQ=255)

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

                     ggTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAGCCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATcgcc                                                                    <  1:540222/93‑1 (MQ=255)

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

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

                         aaaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTg                                                                <  2:319951/93‑1 (MQ=255)

                         aaaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTg                                                                <  2:4538933/93‑1 (MQ=255)

                          aaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTgg                                                               <  1:6039424/93‑1 (MQ=255)

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

                            aaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTTTGACGCACTCATCGCCGCTGGCg                                                             >  2:2065372/1‑93 (MQ=255)

                            aaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCg                                                             >  2:2954373/1‑93 (MQ=255)

                             aGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGt                                                            <  2:1826116/93‑1 (MQ=255)

                             aGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGt                                                            <  1:6394733/93‑1 (MQ=255)

                              gCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTa                                                           >  2:5860697/1‑93 (MQ=255)

                              gCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTa                                                           >  2:3854323/1‑93 (MQ=255)

                               ccAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAg                                                          <  1:5020154/93‑1 (MQ=255)

                               ccAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAg                                                          <  1:4882048/93‑1 (MQ=255)

                               ccAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAg                                                          <  2:1550563/93‑1 (MQ=255)

                               ccAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAg                                                          >  2:1141587/1‑93 (MQ=255)

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

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

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

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

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

                                     gTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGcgcgcg                                                    <  1:1677104/93‑1 (MQ=255)

                                         gACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGGGCGTGGt                                                <  1:4197279/93‑1 (MQ=255)

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

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

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

                                         gACCGCCTATGTCACACTCGAACCCTGGAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGt                                                <  1:2982275/93‑1 (MQ=255)

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

                                           ccGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTg                                              >  1:805907/1‑93 (MQ=255)

                                            cGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGc                                             >  2:1668370/1‑93 (MQ=255)

                                             gCCTATGTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGcc                                            <  1:2144260/93‑1 (MQ=255)

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

                                                 aTGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTTGCGTAGCGCGCGTGGTTGCCTCGa                                        >  1:4944377/1‑93 (MQ=255)

                                                 aTGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGTCTCGa                                        >  1:4413574/1‑93 (MQ=255)

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

                                                  tGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAt                                       >  2:659905/1‑93 (MQ=255)

                                                    tCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGc                                     <  2:2791855/93‑1 (MQ=255)

                                                    tCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGc                                     >  1:2763674/1‑93 (MQ=255)

                                                      acacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCaa                                   >  1:3063964/1‑93 (MQ=255)

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

                                                       cacTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                  <  2:2532920/93‑1 (MQ=255)

                                                        acTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                 >  2:341456/1‑93 (MQ=255)

                                                        acTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                 <  2:6022215/93‑1 (MQ=255)

                                                        acTCGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                 >  2:5578178/1‑93 (MQ=255)

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

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

                                                           cGAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATcc                              >  1:2382521/1‑93 (MQ=255)

                                                            gAACCCTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCt                             <  1:6670907/93‑1 (MQ=255)

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

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

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

                                                                 cTGTAGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAccc                        >  2:4048633/1‑93 (MQ=255)

                                                                     aGCCATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAg                    >  2:4671996/1‑93 (MQ=255)

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

                                                                       ccATCATGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGt                  <  1:129622/93‑1 (MQ=255)

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

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

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

                                                                          tcatGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc               >  1:4712961/1‑93 (MQ=255)

                                                                           catGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCt              <  2:5149786/93‑1 (MQ=255)

                                                                           catGGTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCt              <  2:5513464/93‑1 (MQ=255)

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

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

                                                                              ggTCGTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTggg           <  1:2452764/93‑1 (MQ=255)

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

                                                                                 cGTACGCCAGCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt        <  1:3453816/93‑1 (MQ=255)

                                                                                  gTACGCCATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTg       <  1:4683927/93‑1 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  >  2:1508626/1‑93 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  >  1:300237/1‑93 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  >  2:6457724/1‑93 (MQ=255)

                                                                                       ccATCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACtt  >  1:3592564/1‑93 (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` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`