BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,095,111	G→A	intergenic (‑82/‑206)	rhaB ← / → rhaS	rhamnulokinase/RhaS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,111	0	G	A	98.9%	982.0 / ‑2.6	277	intergenic (‑82/‑206)	rhaB/rhaS	rhamnulokinase/RhaS
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (180/94); total (182/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 280 total aligned reads displayed.

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc                                                                                        >  1:296444‑M2/12‑93 (MQ=255)

  cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt                                                                                      >  1:655650‑M2/10‑93 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCatg                                                                                     >  2:1084219‑M2/9‑92 (MQ=255)

   gacatatgGTGAACCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  1:3489429‑M2/85‑1 (MQ=255)

      atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGc                                                                                  >  1:3273779‑M2/6‑93 (MQ=255)

        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCaa                                                                                >  1:3791236‑M2/4‑93 (MQ=255)

         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATGGGCAAc                                                                               >  1:1660504‑M2/3‑93 (MQ=255)

          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc                                                                              >  2:4549461‑M2/2‑93 (MQ=255)

          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc                                                                              >  1:4525058‑M2/2‑93 (MQ=255)

           gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa                                                                             >  2:1991508/1‑93 (MQ=255)

            tgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAg                                                                            >  1:2190698/1‑93 (MQ=255)

            tgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAg                                                                            <  2:48004/93‑1 (MQ=255)

             gATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAgg                                                                           <  1:1955616/93‑1 (MQ=255)

                 ctgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaaa                                                                       <  1:185593/93‑1 (MQ=255)

                 ctgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaaa                                                                       >  2:1451154/1‑93 (MQ=255)

                  tgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGCACATGAGCCCTTGGCAACCAGGGAAAg                                                                      >  1:2471258/1‑93 (MQ=255)

                  tgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAg                                                                      >  2:4425910/1‑93 (MQ=255)

                   gctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGa                                                                     <  2:674536/93‑1 (MQ=255)

                   gctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTAGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGa                                                                     >  2:2955304/1‑93 (MQ=255)

                   gctgAACTTCATTACGACCAGTCTAAAAAGCGCCTGATTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGa                                                                     <  2:255867/93‑1 (MQ=255)

                    ctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGAt                                                                    >  1:155694/1‑93 (MQ=255)

                     tgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATg                                                                   >  1:780405/1‑93 (MQ=255)

                     tgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATg                                                                   >  1:2423808/1‑93 (MQ=255)

                     tgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATg                                                                   <  2:1207878/93‑1 (MQ=255)

                      gAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGa                                                                  >  2:104229/1‑93 (MQ=255)

                        aTTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAAc                                                                >  2:3735443/1‑93 (MQ=255)

                        aTTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAAc                                                                <  2:424439/93‑1 (MQ=255)

                         tttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACg                                                               <  1:3413065/93‑1 (MQ=255)

                          ttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGt                                                              >  2:3375115/1‑93 (MQ=255)

                          ttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGt                                                              <  2:2114378/93‑1 (MQ=255)

                           tCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtg                                                             <  2:1356466/93‑1 (MQ=255)

                             aTTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgat                                                           <  2:3219966/93‑1 (MQ=255)

                                aCGACCAGCCTAAAAAGCGCCTGAATTCGCCACCTTCTCGTTACTGACAGGAAACTGAGCCATTGGCAACCAGGGAAAGATGAACGTgatgat                                                        <  1:1596271/93‑1 (MQ=255)

                                 cGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGatgatg                                                       <  2:4195054/93‑1 (MQ=255)

                                  gACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      >  2:1491100/1‑93 (MQ=255)

                                  gACCAGTCTAAAAAGCGCCTGAACTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      <  2:2600045/93‑1 (MQ=255)

                                   aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     >  1:4285307/1‑93 (MQ=255)

                                    ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAAGGTGATGATGTTc                                                    >  2:3854655/1‑93 (MQ=255)

                                     cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                   >  2:1839722/1‑93 (MQ=255)

                                     cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                   >  1:3599095/1‑93 (MQ=255)

                                     cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                   >  1:1435890/1‑93 (MQ=255)

                                      aGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTcac                                                  <  1:4612441/93‑1 (MQ=255)

                                      aGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTcac                                                  >  1:2562182/1‑93 (MQ=255)

                                       gTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTcaca                                                 >  2:4268839/1‑93 (MQ=255)

                                       gTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTcaca                                                 >  1:3708528/1‑93 (MQ=255)

                                        tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  1:4458363/1‑93 (MQ=255)

                                         cTAAAAAGCGCCTGAATTCGCGACCTTCTTGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAt                                               >  1:419249/1‑93 (MQ=255)

                                         cTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAt                                               >  1:4361389/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              >  2:1019494/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              >  1:1129807/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              >  1:2261394/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              <  1:2859698/93‑1 (MQ=255)

                                           aaaaaGCTCCTGAATTCGCGACCTTCTCGTGACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAttt                                             <  2:551052/93‑1 (MQ=255)

                                             aaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGc                                           <  1:173491/93‑1 (MQ=255)

                                             aaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGc                                           <  2:4376279/93‑1 (MQ=255)

                                             aaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGc                                           >  2:2969974/1‑93 (MQ=255)

                                               aGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTg                                         <  1:4687784/93‑1 (MQ=255)

                                                  gcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAt                                      >  2:1403833/1‑93 (MQ=255)

                                                  gcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAt                                      >  2:3583416/1‑93 (MQ=255)

                                                   ccTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAtt                                     <  1:1726382/93‑1 (MQ=255)

                                                   ccTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAtt                                     >  1:504356/1‑93 (MQ=255)

                                                   ccTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAtt                                     >  1:24019/1‑93 (MQ=255)

                                                   ccTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAtt                                     >  1:3769488/1‑93 (MQ=255)

                                                    cTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtg                                    <  1:4232647/93‑1 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  2:285059/1‑93 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  2:3759325/1‑93 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  1:266807/1‑93 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   <  1:4059386/93‑1 (MQ=255)

                                                      gAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgtg                                  <  1:1631279/93‑1 (MQ=255)

                                                      gAATTCGCCACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgtg                                  <  1:879698/93‑1 (MQ=255)

                                                       aaTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtgg                                 >  1:101201/1‑93 (MQ=255)

                                                        aTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggt                                <  1:4101872/93‑1 (MQ=255)

                                                          tCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtga                              >  1:385532/1‑93 (MQ=255)

                                                           cgcgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgat                             <  2:2812014/93‑1 (MQ=255)

                                                             cgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgt                           <  1:2829487/93‑1 (MQ=255)

                                                              gACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtg                          >  2:327524/1‑93 (MQ=255)

                                                               aCCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         >  2:1572142/1‑93 (MQ=255)

                                                               aCCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         <  1:212987/93‑1 (MQ=255)

                                                                 cTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGtgatgtgatg                       >  2:839799/1‑93 (MQ=255)

                                                                   tctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCt                     >  1:2454348/1‑93 (MQ=255)

                                                                   tctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCt                     <  2:3074114/93‑1 (MQ=255)

                                                                   tctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCt                     >  2:3449832/1‑93 (MQ=255)

                                                                    ctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTc                    >  2:3962595/1‑93 (MQ=255)

                                                                    ctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTc                    >  1:3153282/1‑93 (MQ=255)

                                                                     tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   <  2:87385/93‑1 (MQ=255)

                                                                     tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   >  2:1893530/1‑93 (MQ=255)

                                                                      cGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAc                  <  1:1263313/93‑1 (MQ=255)

                                                                         tACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGc               >  2:1650503/1‑93 (MQ=255)

                                                                           cTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAt             <  1:472327/93‑1 (MQ=255)

                                                                           cTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACGGCAt             >  1:3385808/1‑93 (MQ=255)

                                                                            tGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAtt            <  2:4122493/93‑1 (MQ=255)

                                                                             gACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt           >  2:2800070/1‑93 (MQ=255)

                                                                             gACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttc           >  1:1513961/1‑92 (MQ=255)

                                                                                aGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCt        >  2:2461484/1‑93 (MQ=255)

                                                                                aGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCt        <  2:2039085/93‑1 (MQ=255)

                                                                                 ggAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTg       <  1:48973/93‑1 (MQ=255)

                                                                                  gAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGa      <  1:1791218/93‑1 (MQ=255)

                                                                                  gAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGa      >  1:105093/1‑93 (MQ=255)

                                                                                   aaaaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaa     <  2:4623548/93‑1 (MQ=255)

                                                                                     aaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaaa   >  1:2080559/1‑93 (MQ=255)

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

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: `‑`