BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999,772	T→C	H4R (CAC→CGC)	b1922 ←	sigma factor 28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,772	0	T	C	98.6%	453.8 / ‑3.8	145	H4R (CAC→CGC)	b1922	sigma factor 28
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (67/76); total (68/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 152 total aligned reads displayed.

GATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999683‑1999861

gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAGCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGtt                                                                                        <  1:4028221/93‑1 (MQ=255)

gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGtt                                                                                        <  1:3310195/93‑1 (MQ=255)

 aTCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGttt                                                                                       <  2:4321346/93‑1 (MQ=255)

  tCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTa                                                                                      <  2:2326617/93‑1 (MQ=255)

       aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATCACGCTGCCACAGCGAGCGTTTATCCAt                                                                                 >  2:346252/1‑93 (MQ=255)

       aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAt                                                                                 >  2:346253/1‑93 (MQ=255)

       aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAt                                                                                 <  2:4172021/93‑1 (MQ=255)

        aGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAtt                                                                                >  2:263939/1‑93 (MQ=255)

         gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTa                                                                               <  1:2421650/93‑1 (MQ=255)

                gCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCTGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        <  2:287608/93‑1 (MQ=255)

                gCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        <  1:4425036/93‑1 (MQ=255)

                gCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        >  1:4720621/1‑93 (MQ=255)

                gCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        <  2:2491786/93‑1 (MQ=255)

                 cgcgcgGGCAGTCGAACGTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCa                                                                       <  1:1411876/91‑1 (MQ=255)

                 cTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCa                                                                       >  2:714400/1‑93 (MQ=255)

                    gcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCg                                                                    <  1:833116/93‑1 (MQ=255)

                      gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:638354/1‑93 (MQ=255)

                      gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  <  2:3818248/93‑1 (MQ=255)

                      gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  2:2598112/1‑93 (MQ=255)

                      gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:1734736/1‑93 (MQ=255)

                         cAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGtata                                                               <  2:1053298/93‑1 (MQ=255)

                          aGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATag                                                              <  2:3773062/93‑1 (MQ=255)

                           gTCGAACCTGCAGGCGCAATGCTTCGTGGCGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  2:3182279/93‑1 (MQ=255)

                           gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             >  2:2413491/1‑93 (MQ=255)

                           gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  1:2544323/93‑1 (MQ=255)

                             cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           <  1:640015/93‑1 (MQ=255)

                             cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           >  2:2026265/1‑93 (MQ=255)

                                aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTgca                                                        >  1:1367728/1‑91 (MQ=255)

                                aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:3872489/1‑93 (MQ=255)

                                aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        >  2:2971003/1‑93 (MQ=255)

                                aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:4647833/1‑93 (MQ=255)

                                 ccTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAt                                                       <  1:3352509/93‑1 (MQ=255)

                                 ccTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAt                                                       <  2:650795/93‑1 (MQ=255)

                                  cTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                                      >  1:1729271/1‑93 (MQ=255)

                                   tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                                     >  1:1621696/1‑93 (MQ=255)

                                     cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATGCATTACACCTTCAGCGGTATAGAGTGAATTCag                                                   >  1:2663551/1‑92 (MQ=255)

                                      aGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACg                                                  >  1:1557929/1‑93 (MQ=255)

                                      aGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACg                                                  >  1:4532055/1‑93 (MQ=255)

                                       ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                                                 <  1:533401/93‑1 (MQ=255)

                                        gcgcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGAt                                                <  1:3719642/93‑1 (MQ=255)

                                          gcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaa                                              >  2:75645/1‑93 (MQ=255)

                                            aaTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                                            >  2:3867925/1‑93 (MQ=255)

                                             aTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  1:2475607/93‑1 (MQ=255)

                                             aTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  2:2408794/93‑1 (MQ=255)

                                              tGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAg                                          <  2:532045/93‑1 (MQ=255)

                                              tGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAg                                          >  2:673407/1‑93 (MQ=255)

                                              tGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAg                                          >  1:3775190/1‑93 (MQ=255)

                                               gCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGc                                         >  1:4746230/1‑93 (MQ=255)

                                                cTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGcc                                        <  1:3519816/93‑1 (MQ=255)

                                                 ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  2:2278217/1‑93 (MQ=255)

                                                 ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  1:404904/1‑93 (MQ=255)

                                                 ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       <  1:635455/93‑1 (MQ=255)

                                                 cgcttGACCCCCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       <  1:635454/89‑1 (MQ=255)

                                                  tCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCt                                      >  1:2848028/1‑93 (MQ=255)

                                                    gTGACGCACCAGCTGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGc                                    <  1:306996/93‑1 (MQ=255)

                                                    gTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGc                                    <  2:3350037/93‑1 (MQ=255)

                                                    gTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAAAGCCCTGc                                    >  1:2329208/1‑93 (MQ=255)

                                                     tGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCg                                   <  2:2023032/93‑1 (MQ=255)

                                                      gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  >  2:1265076/1‑93 (MQ=255)

                                                        cGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTta                                <  2:1676227/93‑1 (MQ=255)

                                                           aCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTatat                             <  2:4543247/93‑1 (MQ=255)

                                                            ccAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATg                            >  1:1327335/1‑93 (MQ=255)

                                                             cAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGa                           <  1:2455505/93‑1 (MQ=255)

                                                             aaGCGGGCAAAAACGCTGTCACAGGGAGCTTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATACACAGCCCTGCGTTATATGa                           <  1:3535029/92‑1 (MQ=255)

                                                               gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         >  1:3330299/1‑93 (MQ=255)

                                                               gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         <  2:1884062/93‑1 (MQ=255)

                                                               gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         >  1:2777131/1‑93 (MQ=255)

                                                               gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         >  1:3710923/1‑93 (MQ=255)

                                                                 gggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTa                       >  2:1663792/1‑93 (MQ=255)

                                                                  gtaCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      <  1:301872/91‑1 (MQ=255)

                                                                  ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  2:3271319/1‑93 (MQ=255)

                                                                  ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      <  1:1116084/93‑1 (MQ=255)

                                                                  ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      <  2:796973/93‑1 (MQ=255)

                                                                   gACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATc                     <  2:136199/93‑1 (MQ=255)

                                                                       tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 >  2:4175625/1‑93 (MQ=255)

                                                                       tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 >  1:89833/1‑93 (MQ=255)

                                                                       tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGGATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 >  2:303663/1‑93 (MQ=255)

                                                                        aaCGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCAt                <  2:4522779/93‑1 (MQ=255)

                                                                         aCGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATg               <  1:528018/93‑1 (MQ=255)

                                                                           gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  1:3432666/93‑1 (MQ=255)

                                                                           gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  1:3432801/93‑1 (MQ=255)

                                                                             tGCGACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTa           <  2:2973054/93‑1 (MQ=255)

                                                                             tGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTa           >  2:2487385/1‑93 (MQ=255)

                                                                              gCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTAt          >  1:95092/1‑93 (MQ=255)

                                                                              gCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTAt          >  1:3631699/1‑93 (MQ=255)

                                                                               ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  2:4388480/93‑1 (MQ=255)

                                                                               ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  1:1360682/93‑1 (MQ=255)

                                                                                cacaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATcc        >  1:3680060/1‑93 (MQ=255)

                                                                                 acaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATcct       >  2:1853823/1‑92 (MQ=255)

                                                                                  cggCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGt      <  1:4570023/91‑1 (MQ=255)

                                                                                   agcgagcgTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGtt     <  2:754024/93‑1 (MQ=255)

                                                                                    gCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGttt    >  2:39454/1‑93 (MQ=255)

                                                                                     cGAGCGTTTATCCATTACTCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTc   <  1:1083783/93‑1 (MQ=255)

                                                                                     cGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTc   <  1:2843813/93‑1 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATTCGTTTCt  >  2:4653403/1‑93 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  >  2:2759029/1‑93 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  >  2:2759030/1‑93 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  >  1:2552137/1‑93 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  >  2:4713064/1‑93 (MQ=255)

                                                                                      gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  <  1:843063/93‑1 (MQ=255)

GATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999683‑1999861

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