BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,508	T→C	E79G (GAA→GGA)	yidK ←	putative cotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,508	0	T	C	99.6%	696.0 / ‑4.7	223	E79G (GAA→GGA)	yidK	putative cotransporter
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (121/101); total (122/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 230 total aligned reads displayed.

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:3478484/1‑93 (MQ=255)

aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:1689018/1‑93 (MQ=255)

 gAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        <  2:1661293/93‑1 (MQ=255)

  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       >  2:4476225/1‑93 (MQ=255)

    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  1:1004715/1‑93 (MQ=255)

    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  2:2685130/1‑93 (MQ=255)

     aTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGcc                                                                                    >  1:3098505/1‑93 (MQ=255)

      tCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGccc                                                                                   >  1:2151122/1‑93 (MQ=255)

        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 >  1:658241/1‑93 (MQ=255)

        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 <  1:3373522/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  2:3689560/93‑1 (MQ=255)

             tggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            >  1:3087573/1‑93 (MQ=255)

               gtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGAc                                                                          >  1:2059568/1‑93 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:3514045/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:1379366/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:4002050/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:1423851/93‑1 (MQ=255)

                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        <  1:1832419/93‑1 (MQ=255)

                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        >  2:2044387/1‑93 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATGACCGACATGcc                                                                     >  1:4213119/1‑93 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     >  2:2179568/1‑93 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  2:2422075/1‑93 (MQ=255)

                      tGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGc                                                                   >  1:866405/1‑93 (MQ=255)

                       gCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCt                                                                  >  2:1275378/1‑93 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:3227078/93‑1 (MQ=255)

                           cgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTg                                                              <  1:3371548/93‑1 (MQ=255)

                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            <  2:2164979/93‑1 (MQ=255)

                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            <  2:1623203/93‑1 (MQ=255)

                              ttCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAg                                                           >  1:3045542/1‑93 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  2:915167/1‑93 (MQ=255)

                                    tAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGc                                                     >  2:1586534/1‑93 (MQ=255)

                                      gcgcGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGccc                                                   <  1:2272712/93‑1 (MQ=255)

                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGTCTGCCCgg                                                 >  1:1220799/1‑93 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  1:4478960/93‑1 (MQ=255)

                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             <  1:1000592/93‑1 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACGTGCCGCTTTTGTAGGACTGCCCGGAAAgg                                            >  2:1533727/1‑93 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:2281458/1‑93 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:4169260/1‑93 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:4210806/1‑93 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           >  1:4064979/1‑93 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           <  1:1662580/93‑1 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          >  2:486229/1‑93 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          <  1:3143322/93‑1 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:4568093/1‑93 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:2561434/1‑93 (MQ=255)

                                                  gATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAc                                       >  1:604887/1‑93 (MQ=255)

                                                  aaTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAc                                       <  2:3089266/92‑1 (MQ=255)

                                                   aTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAcc                                      <  2:3747564/93‑1 (MQ=255)

                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCCGCCCATCACCGACATGCAGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     >  1:2100422/1‑93 (MQ=255)

                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     <  1:1938035/93‑1 (MQ=255)

                                                     tAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAg                                    <  2:2859294/93‑1 (MQ=255)

                                                      aGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgt                                   >  1:3773266/1‑93 (MQ=255)

                                                       gcgcGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgtt                                  >  1:4670129/1‑93 (MQ=255)

                                                         gcgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgt                                >  1:2038400/1‑93 (MQ=255)

                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               <  2:3171027/93‑1 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              <  1:2173105/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  1:3607000/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             >  2:2486960/1‑93 (MQ=255)

                                                             ggAAGATCAGCGTAACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  2:384043/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  1:2456081/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  1:4528115/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  1:4299708/93‑1 (MQ=255)

                                                               aaGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTg                          <  1:1118334/93‑1 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  2:140980/1‑93 (MQ=255)

                                                                  aTCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaa                       >  2:1466219/1‑93 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  1:2879871/93‑1 (MQ=255)

                                                                   tCAGCATCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  2:4312215/93‑1 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    >  2:948744/1‑93 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    <  2:1857318/93‑1 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    >  2:1007729/1‑93 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  2:2607610/1‑93 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  1:1506154/1‑93 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  1:3905845/1‑93 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:1333322/1‑93 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  2:74922/1‑93 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  2:4351405/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:3688849/1‑93 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  2:3505257/1‑93 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTACGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:2516157/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCCACATGCCGCTTTTGTAGGACTGCCCGGGAACGCCGACCAGGTGGTCCGTGGAAAggttg                 >  1:3496112/1‑93 (MQ=255)

                                                                        gTCACCGCTCAAGTCGCTCCCCAGCCCATCCGCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:1746345/93‑1 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  1:3984409/1‑93 (MQ=255)

                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGCAAggctggt               >  2:3436928/1‑93 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             >  2:2373833/1‑93 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             <  1:1550071/93‑1 (MQ=255)

                                                                             cGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAg            >  2:627778/1‑93 (MQ=255)

                                                                              gCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGc           >  1:2994908/1‑93 (MQ=255)

                                                                               cTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCa          <  1:258707/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  2:289627/1‑93 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         <  1:4306469/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         <  2:4565694/93‑1 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        <  2:85419/93‑1 (MQ=255)

                                                                                    gTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  1:1353385/1‑93 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGTAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  2:2649901/93‑1 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  1:1867707/1‑93 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  1:797927/93‑1 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  2:3270175/1‑93 (MQ=255)

                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGa   <  1:3845702/93‑1 (MQ=255)

                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGa   <  2:3387807/93‑1 (MQ=255)

                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCAGTGGAAAGGTTGGTTAGCATTAACGaa  <  2:1130823/93‑1 (MQ=255)

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

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