BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,496	+C	coding (301/408 nt)	phnQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,495	1	.	C	99.3%	1024.0 / ‑1.2	279	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (150/127); total (150/129)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 304 total aligned reads displayed.

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

 aGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑c                                                                                        <  1:2719260/93‑1 (MQ=255)

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTcc                                                                                   >  2:488312/1‑93 (MQ=255)

        aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCg                                                                                  >  2:1615805/1‑93 (MQ=255)

          tAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGc                                                                                >  1:4270473/1‑93 (MQ=255)

          tAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGc                                                                                >  2:567112/1‑93 (MQ=255)

              gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAg                                                                            >  2:5354168/1‑93 (MQ=255)

               tgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAgg                                                                           >  1:5790403/1‑93 (MQ=255)

                ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGc                                                                          <  1:1036451/93‑1 (MQ=255)

                 gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGcc                                                                         <  1:6233339/93‑1 (MQ=255)

                  aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCg                                                                        <  2:3039649/93‑1 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGa                                                                       <  2:4763378/93‑1 (MQ=255)

                    agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGat                                                                      >  2:5893903/1‑93 (MQ=255)

                     gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGata                                                                     <  1:5268511/93‑1 (MQ=255)

                     gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGata                                                                     <  1:569656/93‑1 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGatat                                                                    >  2:4160655/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  <  1:488312/93‑1 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  <  2:1822317/93‑1 (MQ=255)

                         ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCGCCCTGCTCACGCTGTCA‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 >  1:6556237/1‑93 (MQ=255)

                            tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcga                                                              <  2:5772774/93‑1 (MQ=255)

                             ttttGAGGTTGGCTAATGAAACGCTTCCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  1:5360005/93‑1 (MQ=255)

                             ttttGAGGTTGCCTAATGAACCGCTTCCCGCTATTTCTTTTATTCCCCCTGCTCCCGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  1:3475916/93‑1 (MQ=255)

                               ttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgacga                                                           <  2:1085850/93‑1 (MQ=255)

                                  aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGacac                                                        <  1:3903098/93‑1 (MQ=255)

                                   ggTTGTCTAAGGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACAcc                                                       <  1:258352/93‑1 (MQ=255)

                                   ggTTGGCTAATGAAACGCTTTCCGCTCTTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTACGGCACAGGCCGATATTATCGACGACAcc                                                       >  2:2253380/1‑93 (MQ=255)

                                     ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                                     >  1:28992/1‑93 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  2:5810917/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  1:1813452/93‑1 (MQ=255)

                                       ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCGCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  2:193170/93‑1 (MQ=255)

                                        gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  <  2:4377461/93‑1 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTTGCaa                                               >  2:6337392/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  1:1540784/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  1:1206669/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  2:4953703/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               <  2:2472050/93‑1 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTGCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               <  2:4508750/93‑1 (MQ=255)

                                            aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              >  1:3966587/1‑93 (MQ=255)

                                            aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              >  1:6352036/1‑93 (MQ=255)

                                            aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              >  1:4608365/1‑93 (MQ=255)

                                             tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAata                                             <  1:82561/93‑1 (MQ=255)

                                             tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAata                                             <  1:560109/93‑1 (MQ=255)

                                              gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            <  1:2974707/93‑1 (MQ=255)

                                              gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            <  1:2154121/93‑1 (MQ=255)

                                                aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          >  2:6486993/1‑93 (MQ=255)

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCGTTGGCAATATTcag                                        <  1:5053902/93‑1 (MQ=255)

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                        <  1:5955327/93‑1 (MQ=255)

                                                    cTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagca                                      <  2:545741/93‑1 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     <  2:3503287/93‑1 (MQ=255)

                                                      ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                    >  1:2113694/1‑93 (MQ=255)

                                                      ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                    >  2:1729470/1‑93 (MQ=255)

                                                      ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                    <  1:3073260/93‑1 (MQ=255)

                                                        ccGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                                  >  2:3910667/1‑93 (MQ=255)

                                                        ccGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                                  <  1:2011708/93‑1 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 >  1:1418637/1‑93 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 >  2:2512482/1‑93 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 <  2:4815266/93‑1 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  1:2354253/1‑93 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                <  2:1246152/93‑1 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  2:4651821/1‑93 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                <  1:3616490/93‑1 (MQ=255)

                                                           cTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGact                               <  1:4190319/93‑3 (MQ=255)

                                                            tATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTa                              >  1:6269508/1‑93 (MQ=255)

                                                              tttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt                            >  1:5027207/1‑93 (MQ=255)

                                                                tCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACTACACCATTGGCAATATTCAGCAGGCGATTAATGa                          >  1:2191736/1‑93 (MQ=255)

                                                                tCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  2:1340145/93‑1 (MQ=255)

                                                                  ttttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACg                        <  1:3551212/93‑1 (MQ=255)

                                                                  ttttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACg                        >  2:3552141/1‑93 (MQ=255)

                                                                   tttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       >  2:6582205/1‑93 (MQ=255)

                                                                    ttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                      >  1:931511/1‑93 (MQ=255)

                                                                     tATTCACCCTGCTCACGCTGTCCTC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     <  2:2161211/93‑1 (MQ=255)

                                                                     tATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  2:783197/1‑93 (MQ=255)

                                                                      aTTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                    >  2:2943889/1‑93 (MQ=255)

                                                                      aTTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGAAATATTCAGCAGGCGATTAATGACGCTta                    >  1:1005400/1‑93 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   >  2:3440872/1‑93 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   <  1:4416096/93‑1 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACAACATTGGCAATATTCAGTATGCGCTTAATGACGCTtat                   >  1:34172/1‑93 (MQ=255)

                                                                         cACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 >  1:2639377/1‑93 (MQ=255)

                                                                          aCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                <  1:4126551/93‑1 (MQ=255)

                                                                          aCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                >  2:406748/1‑93 (MQ=255)

                                                                          aCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                >  2:3749397/1‑93 (MQ=255)

                                                                            ccTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              <  1:1700448/93‑1 (MQ=255)

                                                                             cTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt             <  2:4203795/93‑1 (MQ=255)

                                                                              tGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg            >  1:4452872/1‑93 (MQ=255)

                                                                               gCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGa           >  1:254198/1‑93 (MQ=255)

                                                                               gCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGa           <  2:5047419/93‑1 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  2:1408223/93‑1 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          >  2:6000953/1‑93 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          >  1:1758124/1‑93 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          >  2:4298266/1‑93 (MQ=255)

                                                                                 ccccGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         <  2:2743611/90‑1 (MQ=255)

                                                                                  cACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCg        >  2:277886/1‑93 (MQ=255)

                                                                                   aCGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       <  1:3745329/93‑1 (MQ=255)

                                                                                    cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      >  1:1601676/1‑93 (MQ=255)

                                                                                    cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      >  2:5932572/1‑93 (MQ=255)

                                                                                     gCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg     >  2:3279770/1‑93 (MQ=255)

                                                                                      cTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    <  1:6456926/93‑1 (MQ=255)

                                                                                      cTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    <  2:1154794/93‑1 (MQ=255)

                                                                                      cTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    >  1:4807363/1‑93 (MQ=255)

                                                                                       tGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGc   >  2:3104642/1‑93 (MQ=255)

                                                                                        gTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCg  >  1:1336053/1‑93 (MQ=255)

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

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