BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,458,272	(G)_5→4	coding (218/1311 nt)	b2342 ←	putative acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,458,268	0	G	.	99.1%	1007.1 / ‑1.0	225	coding (222/1311 nt)	b2342	putative acyltransferase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base . (115/108); total (115/110)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 243 total aligned reads displayed.

GCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458176‑2458356

gCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTg                                                                                          <  2:5892144/93‑1 (MQ=255)

gCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTg                                                                                          >  1:5421902/1‑93 (MQ=255)

  cTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        >  2:4982973/1‑93 (MQ=255)

   tGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGAATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTgggg                                                                                       >  2:2422905/1‑93 (MQ=255)

    ggAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggc                                                                                      >  2:2236384/1‑92 (MQ=255)

      aaaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggctt                                                                                    >  1:2943911/1‑90 (MQ=255)

          tGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAgg                                                                               <  2:3303241/93‑1 (MQ=255)

            gTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCa                                                                             >  1:1968918/1‑93 (MQ=255)

            gTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCa                                                                             <  1:4633039/93‑1 (MQ=255)

             tAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAt                                                                            <  1:7091722/93‑1 (MQ=255)

              aGCGCAAGCGCGGCTGACGGTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAtt                                                                           <  2:2376130/93‑1 (MQ=255)

               gcgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAttt                                                                          >  1:6077852/1‑93 (MQ=255)

               gcgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAttt                                                                          <  1:2149065/93‑1 (MQ=255)

                cgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTg                                                                         <  2:878041/93‑1 (MQ=255)

                 gcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGt                                                                        <  1:2600840/93‑1 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       <  1:4969415/93‑1 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       >  1:6372783/1‑93 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       <  2:1034969/93‑1 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       <  2:7097280/93‑1 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       <  2:1520268/93‑1 (MQ=255)

                   aaGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  1:5834328/93‑1 (MQ=255)

                   aaGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  2:3180704/93‑1 (MQ=255)

                    aGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACg                                                                     >  1:1810256/1‑93 (MQ=255)

                    aGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACg                                                                     >  1:4817823/1‑93 (MQ=255)

                     gcgcGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGa                                                                    <  2:6912122/93‑1 (MQ=255)

                      cgcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAc                                                                   >  1:6499225/1‑93 (MQ=255)

                        cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCt                                                                 >  1:2800824/1‑93 (MQ=255)

                        cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCt                                                                 >  2:5569621/1‑93 (MQ=255)

                            tGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGAcc                                                             >  1:6691265/1‑93 (MQ=255)

                            tGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGAcc                                                             >  1:5790289/1‑93 (MQ=255)

                            tGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGAcc                                                             >  1:3428449/1‑93 (MQ=255)

                              aCGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCaa                                                           <  1:3289144/93‑1 (MQ=255)

                                gCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         >  2:5479238/1‑93 (MQ=255)

                                gCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         <  2:3482596/93‑1 (MQ=255)

                                 cTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGa                                                        >  1:7681726/1‑93 (MQ=255)

                                  tGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGAc                                                       >  2:1796588/1‑93 (MQ=255)

                                   gTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGAcc                                                      >  2:3923205/1‑93 (MQ=255)

                                    tAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCa                                                     <  2:1253273/93‑1 (MQ=255)

                                     aaGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAg                                                    <  1:403933/93‑1 (MQ=255)

                                     aaGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAg                                                    >  2:6771211/1‑93 (MQ=255)

                                         aTCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgt                                                >  1:3527924/1‑93 (MQ=255)

                                         aTCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgt                                                >  1:718590/1‑93 (MQ=255)

                                          tCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgtt                                               >  2:7858052/1‑93 (MQ=255)

                                           cGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTc                                              <  2:800953/93‑1 (MQ=255)

                                            ggTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCa                                             <  1:807780/93‑1 (MQ=255)

                                            ggTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCa                                             >  2:5072480/1‑93 (MQ=255)

                                            cgTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCa                                             <  1:3238806/92‑1 (MQ=255)

                                              tATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAAt                                           <  1:3779984/93‑1 (MQ=255)

                                              tATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAAt                                           <  1:4740417/93‑1 (MQ=255)

                                              tATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAAt                                           <  2:7377114/93‑1 (MQ=255)

                                               aTGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATc                                          <  1:7205557/93‑1 (MQ=255)

                                                tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAAtga                                         >  2:1540067/1‑91 (MQ=255)

                                                tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCa                                         >  1:5572524/1‑93 (MQ=255)

                                                tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCa                                         >  1:3853924/1‑93 (MQ=255)

                                                  tccattcgttCCCGTGCCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACt                                       <  1:3668857/91‑1 (MQ=255)

                                                   aCATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACtt                                      >  1:4210063/1‑93 (MQ=255)

                                                    cattcattCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTc                                     <  2:3007433/93‑1 (MQ=255)

                                                    cattcattCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTc                                     >  2:2800288/1‑93 (MQ=255)

                                                     attcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcg                                    >  1:6524956/1‑93 (MQ=255)

                                                      ttcattcCCGTACCGAGAACAATTTCACTCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   >  1:4399997/1‑93 (MQ=255)

                                                      ttcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   <  1:7932768/93‑1 (MQ=255)

                                                         attcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcggcgg                                >  2:3046047/1‑93 (MQ=255)

                                                         attcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcggcgg                                >  1:6285823/1‑93 (MQ=255)

                                                          ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               >  2:2088006/1‑93 (MQ=255)

                                                           tcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCgggg                              >  2:6019253/1‑93 (MQ=255)

                                                             ccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAt                            >  2:2460297/1‑93 (MQ=255)

                                                             ccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAt                            >  1:2743894/1‑93 (MQ=255)

                                                              cGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtc                           <  2:5153936/93‑1 (MQ=255)

                                                               gTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtct                          <  2:6465056/93‑1 (MQ=255)

                                                                  ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       <  1:1207225/93‑1 (MQ=255)

                                                                  ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       <  2:1482877/93‑1 (MQ=255)

                                                                  ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       >  2:1218521/1‑93 (MQ=255)

                                                                   cGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCt                      >  1:6229064/1‑93 (MQ=255)

                                                                    gagaACAATTTCAAGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTg                     <  2:2703971/93‑1 (MQ=255)

                                                                     agaACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTg                     <  1:228782/93‑1 (MQ=255)

                                                                      gaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCg                   >  2:846106/1‑93 (MQ=255)

                                                                        aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  2:6683620/1‑93 (MQ=255)

                                                                        aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  1:5292150/1‑93 (MQ=255)

                                                                           aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGGTCTCGCTGCGTGccg              >  1:2548235/1‑92 (MQ=255)

                                                                           aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGccc              >  2:3947368/1‑92 (MQ=255)

                                                                           aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCca              <  2:806374/93‑1 (MQ=255)

                                                                           aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCca              <  2:7843018/93‑1 (MQ=255)

                                                                           aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCca              <  1:1617968/93‑1 (MQ=255)

                                                                            tttCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcccg             >  2:1939897/1‑91 (MQ=255)

                                                                            tttCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcag             >  2:1141196/1‑93 (MQ=255)

                                                                             ttCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagc            <  2:231427/93‑1 (MQ=255)

                                                                              tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           >  2:4576322/1‑93 (MQ=255)

                                                                              tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           >  2:4175910/1‑93 (MQ=255)

                                                                               taCGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagcag          <  1:6817918/92‑1 (MQ=255)

                                                                               cACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagccg          >  1:7469177/1‑91 (MQ=255)

                                                                               cACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagcag          <  1:410007/93‑1 (MQ=255)

                                                                                aCGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGt         >  1:6147233/1‑93 (MQ=255)

                                                                                aCGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGt         >  2:728297/1‑93 (MQ=255)

                                                                                    gcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGc     >  2:6355934/1‑93 (MQ=255)

                                                                                    gcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGc     <  2:640854/93‑1 (MQ=255)

                                                                                     cAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGcc    <  1:5086381/93‑1 (MQ=255)

                                                                                      aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  2:7653004/1‑93 (MQ=255)

                                                                                       aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCtc  >  1:2010841/1‑92 (MQ=255)

                                                                                       aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTa  >  2:1347226/1‑93 (MQ=255)

                                                                                       aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTa  >  2:2604311/1‑93 (MQ=255)

GCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458176‑2458356

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