BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,044,033	(C)_5→6	coding (1149/7152 nt)	b1978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,044,028	1	.	C	96.7%	683.7 / 13.6	215	coding (1144/7152 nt)	b1978	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/5); new base C (118/90); total (120/95)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 236 total aligned reads displayed.

CTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑CCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043937‑2044115

cTATGAACACTATTATGGCGATAAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAACCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑c                                                                                         <  2:860182/93‑1 (MQ=255)

 tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cc                                                                                        >  1:2362517/1‑93 (MQ=255)

 tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cc                                                                                        <  2:2398807/93‑1 (MQ=255)

      aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATGCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                   >  2:4490746/1‑91 (MQ=255)

      aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                   <  1:5082458/93‑3 (MQ=255)

       cAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccctt                                                                                  >  1:2733654/1‑90 (MQ=255)

       cAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccctt                                                                                  >  2:1870533/1‑90 (MQ=255)

          tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTccc                                                                                <  2:665310/93‑1 (MQ=255)

          tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTccc                                                                                >  1:3754416/1‑93 (MQ=255)

          tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTccc                                                                                <  1:1718632/93‑1 (MQ=255)

            ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑Gc                                                                             >  1:1200112/1‑93 (MQ=255)

             tatGGCGATGAAGTGCCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCt                                                                            <  1:1340796/93‑1 (MQ=255)

              atGGCGATGAAGTGGCCTTGTTCTATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCtg                                                                           <  2:3974108/93‑1 (MQ=255)

                   gATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGAc                                                                      >  1:2080751/1‑93 (MQ=255)

                     tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCt                                                                    <  1:2448577/93‑1 (MQ=255)

                     tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCt                                                                    <  1:5325191/93‑1 (MQ=255)

                       aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTc                                                                  <  1:2203952/93‑1 (MQ=255)

                        aGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCa                                                                 <  1:6417425/93‑1 (MQ=255)

                         gTGGCCCTGTTCGATAAAGACGATCGGCAAATTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAg                                                                >  2:4825645/1‑93 (MQ=255)

                           ggCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAgcg                                                              >  1:6616878/1‑93 (MQ=255)

                              ccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCgg                                                           >  1:2489258/1‑93 (MQ=255)

                               cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGa                                                          >  1:826668/1‑93 (MQ=255)

                                 gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGc                                                        <  2:167949/93‑1 (MQ=255)

                                   tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCaa                                                      >  1:3403948/1‑93 (MQ=255)

                                   tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCaa                                                      <  1:1093265/93‑1 (MQ=255)

                                     gATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACg                                                    >  2:3016347/1‑93 (MQ=255)

                                      aTAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGc                                                   <  1:4001711/93‑1 (MQ=255)

                                       tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGcc                                                  >  1:3900921/1‑93 (MQ=255)

                                       tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGcc                                                  >  1:88381/1‑93 (MQ=255)

                                       tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGcc                                                  <  1:5548638/93‑1 (MQ=255)

                                         aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                <  1:684133/93‑1 (MQ=255)

                                         aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                <  1:6052850/93‑1 (MQ=255)

                                         aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                >  2:5152958/1‑93 (MQ=255)

                                         aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                >  2:6151859/1‑93 (MQ=255)

                                         aaGACGATCGGCAAAGTAACCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                <  2:3854702/93‑1 (MQ=255)

                                          aGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAgg                                               <  1:3344546/93‑1 (MQ=255)

                                          aGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACATAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAgg                                               >  2:4016452/1‑93 (MQ=255)

                                           gACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAggg                                              >  1:2601830/1‑93 (MQ=255)

                                             cGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTa                                            >  1:2894685/1‑93 (MQ=255)

                                               aTCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTaaa                                          <  2:3187691/93‑1 (MQ=255)

                                                tCGGCAAAGTACTCCTCATGCCATAACCGCTGGCCTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAAc                                         <  2:251896/93‑1 (MQ=255)

                                                tCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAAc                                         <  2:5982987/93‑1 (MQ=255)

                                                 cGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACa                                        <  2:6408269/93‑1 (MQ=255)

                                                 cGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACa                                        <  1:3556661/93‑1 (MQ=255)

                                                  ggCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       >  1:5383102/1‑93 (MQ=255)

                                                   gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:6496072/93‑1 (MQ=255)

                                                   gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      >  2:6049983/1‑93 (MQ=255)

                                                   gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  2:4239702/93‑1 (MQ=255)

                                                   gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAACGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:474307/93‑1 (MQ=255)

                                                     aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGc                                    >  2:5177943/1‑93 (MQ=255)

                                                     aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGc                                    >  2:6213927/1‑93 (MQ=255)

                                                      aaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCg                                   >  2:1896371/1‑93 (MQ=255)

                                                      aaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCg                                   <  1:3599190/93‑1 (MQ=255)

                                                      aaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCg                                   >  1:5226415/1‑93 (MQ=255)

                                                       aGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGa                                  >  2:1532310/1‑93 (MQ=255)

                                                        gTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaa                                 >  1:4580589/1‑93 (MQ=255)

                                                          aaTCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaaa                               <  2:3918794/93‑1 (MQ=255)

                                                          aaTCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaaa                               >  1:2527818/1‑93 (MQ=255)

                                                          aaTCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaaa                               <  2:382788/93‑1 (MQ=255)

                                                             ccTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGa                            <  2:4738102/93‑1 (MQ=255)

                                                              cTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGac                           >  1:5652725/1‑93 (MQ=255)

                                                               tCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGaca                          >  2:4530475/1‑93 (MQ=255)

                                                               tCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGaca                          >  2:5853721/1‑93 (MQ=255)

                                                                cATGCCATAACCGCTGGACTTAACTATCC‑CCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGacac                         <  1:1897894/93‑1 (MQ=255)

                                                                 aTGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACAcc                        <  2:2075656/93‑1 (MQ=255)

                                                                    ccATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGCAAATGACACCCGt                     >  1:581110/1‑93 (MQ=255)

                                                                     cATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtt                    >  1:119072/1‑93 (MQ=255)

                                                                     cATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtt                    >  2:3500248/1‑93 (MQ=255)

                                                                     cATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtt                    >  2:101744/1‑93 (MQ=255)

                                                                      aTAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGttt                   <  1:4625284/93‑1 (MQ=255)

                                                                      aTAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGttt                   >  2:4160044/1‑93 (MQ=255)

                                                                      aTAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGttt                   >  1:4793379/1‑93 (MQ=255)

                                                                       tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  >  2:6475474/1‑93 (MQ=255)

                                                                       tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  >  1:2243618/1‑93 (MQ=255)

                                                                       tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  <  1:5714672/93‑1 (MQ=255)

                                                                        aaCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTg                 >  1:1319278/1‑93 (MQ=255)

                                                                        aaCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTg                 >  1:5410590/1‑93 (MQ=255)

                                                                         aCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                >  1:5924112/1‑93 (MQ=255)

                                                                         aCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                >  1:2836164/1‑93 (MQ=255)

                                                                         aCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                <  2:4944853/93‑1 (MQ=255)

                                                                         aCCGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                >  2:4763374/1‑93 (MQ=255)

                                                                          ccGCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               >  2:656619/1‑93 (MQ=255)

                                                                            gCTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGt             <  2:3267333/93‑1 (MQ=255)

                                                                             cTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTc            >  2:5629421/1‑93 (MQ=255)

                                                                             cTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTc            <  1:667435/93‑1 (MQ=255)

                                                                             cTGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTc            >  1:1515840/1‑93 (MQ=255)

                                                                              tGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCg           <  1:4058657/93‑1 (MQ=255)

                                                                              tGGACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCg           >  2:164348/1‑93 (MQ=255)

                                                                               ggACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGa          >  2:2263036/1‑93 (MQ=255)

                                                                                gACTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAt         >  2:6268916/1‑93 (MQ=255)

                                                                                gACTTAACTAT‑A‑CCCCCCTCCCGGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAt         <  1:4891658/93‑1 (MQ=255)

                                                                                 aCTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtt        >  2:370843/1‑93 (MQ=255)

                                                                                  cTTAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAttt       >  1:4357523/1‑93 (MQ=255)

                                                                                   ttAACTATGC‑CCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtttt      <  2:2343087/93‑1 (MQ=255)

                                                                                   ttAACTATCC‑CCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtttt      <  2:3056429/93‑1 (MQ=255)

                                                                                   ttAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtttt      <  2:1343205/93‑1 (MQ=255)

                                                                                    tAACTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTa     <  2:5569432/93‑1 (MQ=255)

                                                                                     aaCTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTac    >  2:4345126/1‑93 (MQ=255)

                                                                                      aCTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacc   >  1:2658838/1‑93 (MQ=255)

                                                                                       cTAT‑ACCCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacct  <  2:4586434/93‑1 (MQ=255)

CTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑CCCCCTTCCC‑GCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043937‑2044115

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