BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	97.8%	1280.9 / 16.5	405	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (9/0); new base C (192/204); total (201/204)
Fisher's exact test for biased strand distribution p-value = 1.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 447 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                             <  2:8365426/93‑1 (MQ=255)

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                             <  1:517442/93‑1 (MQ=255)

 gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                           <  1:3290531/93‑1 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑c                                                                                         <  2:9598157/93‑1 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑c                                                                                         <  2:9051845/93‑1 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATCGCGT‑c‑c‑ccaa                                                                                    >  1:6534654/1‑90 (MQ=255)

       cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑c‑ccaac                                                                                   <  1:4178875/93‑5 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGc                                                                            >  1:1642541/1‑93 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGc                                                                            >  2:6490902/1‑93 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGc                                                                            >  1:3837777/1‑93 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CAC‑C‑AACCCCACGc                                                                            >  1:2297286/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGcc                                                                           >  2:7474107/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGcc                                                                           >  2:1710888/1‑93 (MQ=255)

                aTTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAGATAGCGT‑CTC‑C‑AACCCCACGCCa                                                                          >  2:5525578/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCaa                                                                         <  2:1809137/93‑1 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCaa                                                                         <  1:8341714/93‑1 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCaa                                                                         <  2:3523304/93‑1 (MQ=255)

                   cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAAcc                                                                       <  2:2475134/93‑1 (MQ=255)

                   cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAAcc                                                                       <  1:4465415/93‑1 (MQ=255)

                   cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAAcc                                                                       <  2:4903296/93‑1 (MQ=255)

                     cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAg                                                                     >  1:4315942/1‑93 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGc                                                                    <  1:5998996/93‑1 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAATACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑C‑CGAACCCCACGCCAACCAGca                                                                   >  1:7193793/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                                               <  1:1896296/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAata                                                              <  1:1133613/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAata                                                              <  1:9899381/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCAATGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAata                                                              <  1:7456314/93‑1 (MQ=255)

                             tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑CTC‑AACCCCACGCCAACCAGCACAatat                                                             >  1:1466092/1‑93 (MQ=255)

                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATg                                                            <  2:6650167/93‑1 (MQ=255)

                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATg                                                            >  2:512956/1‑93 (MQ=255)

                                aGTAAGCCTCACCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGaa                                                          <  2:2250978/93‑1 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGaa                                                          <  2:6087068/93‑1 (MQ=255)

                                aGTAAGCCACTCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGaa                                                          <  2:9483123/93‑1 (MQ=255)

                                  tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATc                                                        <  1:2863869/93‑1 (MQ=255)

                                   aaGCCGCGCCTACCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCg                                                       <  1:8778686/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  1:6973296/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  1:5503054/93‑1 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  2:1365520/93‑1 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CAC‑C‑AACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  1:8536287/1‑93 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  2:7199832/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  1:1570576/93‑1 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  2:4040038/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:10031478/93‑1 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  2:9233169/93‑1 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             <  1:6828849/93‑1 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:6870871/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:676018/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATAGCCAGTCCGGt                                             >  1:6475830/1‑93 (MQ=255)

                                              gCCCCCTTCAACGCTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            <  2:792255/93‑1 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  1:6736523/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  1:9808287/93‑1 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:7813067/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:477543/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  2:8515952/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       <  1:5368248/93‑1 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  2:5990163/1‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     <  2:4360190/93‑1 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  1:7376696/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  1:6629293/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:5784397/93‑1 (MQ=255)

                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  <  2:2170885/93‑1 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:4490372/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:4980047/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGCAGCCCCGCaa                               >  2:8188337/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  1:36661/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               <  1:2777868/93‑1 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:3958006/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  2:7891849/93‑1 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  1:2678751/1‑93 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             >  2:1465799/1‑93 (MQ=255)

                                                              ccaaCGCAAACACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg                            <  1:4566172/91‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  2:3763610/93‑1 (MQ=255)

                                                                 aCGCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACg                         <  2:2109656/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  1:9511408/1‑93 (MQ=255)

                                                                    cAATCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa                      >  1:8202332/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  2:2610861/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:4031555/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  1:2024240/93‑1 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  1:2350952/1‑93 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  >  1:5815728/1‑93 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGACAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  >  1:5619789/1‑93 (MQ=255)

                                                                         aCTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:8849174/1‑93 (MQ=255)

                                                                         aCTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:1250666/1‑93 (MQ=255)

                                                                          cTGAGCGGTAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                <  1:6950849/93‑1 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                >  1:2430366/1‑93 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                <  2:8787845/93‑1 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  1:625633/1‑93 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             <  1:4885949/93‑1 (MQ=255)

                                                                                gggAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  2:1136072/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          <  1:3167050/93‑1 (MQ=255)

                                                                                gggAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          <  2:2938173/93‑1 (MQ=255)

                                                                                 ggAAAATAGCGTCC‑C‑C‑AACCCCACGCCAAACAGCACAATATGAATCGCCAGTTAGAACTAAAGCCAAGCAATATACGCCACCCCTCCGCGaaa         >  2:3302209/1‑91 (MQ=255)

                                                                                  gAAAATAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  2:8637648/93‑1 (MQ=255)

                                                                                   aaaaTAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       <  2:3061654/93‑1 (MQ=255)

                                                                                     aaTAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  2:5309380/1‑93 (MQ=255)

                                                                                      aTAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    >  1:9329018/1‑93 (MQ=255)

                                                                                      aTAGCGT‑C‑C‑CTAGCCCCACGCCACCCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  1:574319/93‑1 (MQ=255)

                                                                                       tAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   >  2:6188087/1‑93 (MQ=255)

                                                                                       tAGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   >  2:278077/1‑93 (MQ=255)

                                                                                        aGCGTCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  >  1:9041295/1‑93 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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