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	99.5%	1408.3 / ‑4.8	391	intergenic (‑82/‑206)	rhaB/rhaS	rhamnulokinase/RhaS
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (219/170); total (219/172)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 394 total aligned reads displayed.

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

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

gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc                                                                                        >  2:2604682‑M2/12‑93 (MQ=255)

gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc                                                                                        <  2:5380710‑M2/82‑1 (MQ=255)

 tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa                                                                                       >  2:2227569‑M2/11‑93 (MQ=255)

 tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa                                                                                       >  2:740849‑M2/11‑93 (MQ=255)

  cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt                                                                                      <  2:1266378‑M2/84‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  2:2889502‑M2/85‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  1:5322108‑M2/85‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     >  2:4253499‑M2/9‑93 (MQ=255)

     catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTgg                                                                                   <  1:298580‑M2/87‑1 (MQ=255)

      atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATtgga                                                                                  >  2:4872017‑M2/6‑92 (MQ=255)

         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc                                                                               >  1:1971275‑M2/3‑93 (MQ=255)

         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc                                                                               <  2:1483558‑M2/91‑1 (MQ=255)

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

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

             gagccTGCTGAATTTCATTACGTCCAGTCTAAAAAGCACCTGAATTCGCGCCCTTCTCGTAACGGACAGGAAAATGAGCCATTGGCAACCAgg                                                                           <  2:4742699/90‑1 (MQ=255)

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

              aTCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAggg                                                                          >  1:2569894/1‑93 (MQ=255)

              aTCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAggg                                                                          >  2:5756897/1‑93 (MQ=255)

                ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        <  1:6753188/93‑1 (MQ=255)

                ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        <  1:3140117/93‑1 (MQ=255)

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

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

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

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

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

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

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

                   gctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGa                                                                     >  1:7389175/1‑93 (MQ=255)

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

                    ctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGAt                                                                    <  2:7715762/93‑1 (MQ=255)

                      gAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGa                                                                  <  1:439033/93‑1 (MQ=255)

                      gAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGAAACCAGGGAAAGAtaa                                                                  >  1:2743175/1‑91 (MQ=255)

                       aaTTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGaa                                                                 <  2:3592994/93‑1 (MQ=255)

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

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

                         tttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACg                                                               >  2:5470010/1‑93 (MQ=255)

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

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

                            cATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtga                                                            >  2:7095744/1‑93 (MQ=255)

                               tACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgatga                                                         <  1:7890806/93‑1 (MQ=255)

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

                                  gACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      <  1:6281008/93‑1 (MQ=255)

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

                                   aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTGCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     <  1:965667/93‑1 (MQ=255)

                                    ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                    >  1:7555072/1‑93 (MQ=255)

                                    ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGAAACCAGGGAAAGATGAACGTGATGATGtta                                                    >  1:6629746/1‑92 (MQ=255)

                                     cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                   <  2:5887805/93‑1 (MQ=255)

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

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

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

                                           aaaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAttt                                             >  2:7164450/1‑93 (MQ=255)

                                            aaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtttt                                            >  2:4039085/1‑92 (MQ=255)

                                              aaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCt                                          <  2:3867610/93‑1 (MQ=255)

                                              aaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCt                                          <  2:4344997/93‑1 (MQ=255)

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

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

                                                  gcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAt                                      <  1:4210213/93‑1 (MQ=255)

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

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

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

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

                                                    cTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtg                                    >  2:4621940/1‑93 (MQ=255)

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

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

                                                       aaTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtgg                                 <  2:6130903/93‑1 (MQ=255)

                                                        aTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggt                                >  2:3031889/1‑93 (MQ=255)

                                                         ttCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggtg                               >  1:3550239/1‑93 (MQ=255)

                                                          tCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtga                              <  1:1525899/93‑1 (MQ=255)

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

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

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

                                                             cgACCTTCTCGTTACTGACAGGAAGATGAGCCATTGGCAACCAGGGAAAGATGAACGCGATGATGTTCACAATTTGCTGAATTGTGgtgatgt                           <  2:7798037/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                               aCCTTCTCGTTACAGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATAAACGTGATGATGTGCACAAGTGGCGGAATTGGggtgaggtga                         >  1:6021688/1‑93 (MQ=255)

                                                                ccTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtgat                        >  2:1668325/1‑93 (MQ=255)

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

                                                                  ttCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGc                      >  2:3394802/1‑93 (MQ=255)

                                                                     ctgTTGCTGACAGGAAAATGAGCCAGTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGCATTGTGGTGATGTGATGCTCa                   <  1:1361664/91‑1 (MQ=255)

                                                                      cGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAc                  >  1:6124530/1‑93 (MQ=255)

                                                                       gTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAcc                 <  1:7729699/93‑1 (MQ=255)

                                                                       gTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAcc                 >  2:5207252/1‑93 (MQ=255)

                                                                       gTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAcc                 >  2:3171611/1‑93 (MQ=255)

                                                                       gTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAcc                 <  1:233061/93‑1 (MQ=255)

                                                                        ttACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCg                <  1:5491781/93‑1 (MQ=255)

                                                                          aCTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCa              <  2:5105734/93‑1 (MQ=255)

                                                                          aCTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCa              <  2:5051297/93‑1 (MQ=255)

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

                                                                            tGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAtt            >  2:3767699/1‑93 (MQ=255)

                                                                            tGACAGGAAAAGGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAtt            <  2:5655473/93‑1 (MQ=255)

                                                                              aCAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttta          <  2:2645182/93‑2 (MQ=255)

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

                                                                                     aaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaaa   <  2:1922138/93‑1 (MQ=255)

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

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

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