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	96.7%	322.4 / 5.1	123	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (63/55); total (65/58)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 136 total aligned reads displayed.

TCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCG  >  NC_000913/4006406‑4006550

tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c                                                                        <  2:3697759/76‑1 (MQ=255)

tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c                                                                        <  1:1154921/76‑1 (MQ=255)

 cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc                                                                       >  2:3234143/1‑76 (MQ=255)

   cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc‑cc                                                                    <  1:4304845/76‑2 (MQ=255)

    aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc‑cca                                                                   <  1:1584784/76‑3 (MQ=255)

    aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc‑cca                                                                   >  1:2693767/1‑74 (MQ=255)

     gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGCAAAATAGCG‑T‑cc‑ccaa                                                                  >  2:538123/1‑73 (MQ=255)

         aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAAcccc                                                               >  2:2659472/1‑76 (MQ=255)

         aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAAcccc                                                               >  1:3731848/1‑76 (MQ=255)

         aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAAcccc                                                               >  1:3545053/1‑76 (MQ=255)

           tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCAc                                                             <  2:5487269/76‑1 (MQ=255)

           gcAAGGAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCAc                                                             <  2:258826/75‑1 (MQ=255)

           acAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCAc                                                             <  1:4297200/75‑1 (MQ=255)

            cAAGTAAGCGGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACg                                                            <  2:2948869/76‑1 (MQ=255)

            cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACg                                                            >  2:1939841/1‑76 (MQ=255)

            cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACg                                                            <  2:1997150/76‑1 (MQ=255)

             aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGc                                                           >  1:3728705/1‑76 (MQ=255)

             aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGc                                                           <  2:3576215/76‑1 (MQ=255)

              aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGcc                                                          <  2:2720290/76‑1 (MQ=255)

               gTAAGCCGCGCCTGACCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CCACCCCACGCCa                                                         <  2:1957984/76‑1 (MQ=255)

                tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCaa                                                        <  2:3667545/76‑1 (MQ=255)

                tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCaa                                                        <  1:4810632/76‑1 (MQ=255)

                 aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAAc                                                       >  1:456328/1‑76 (MQ=255)

                 aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAAc                                                       <  1:1018619/76‑1 (MQ=255)

                  aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAAcc                                                      >  1:5187539/1‑76 (MQ=255)

                  aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAAcc                                                      >  2:4885373/1‑76 (MQ=255)

                   gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCa                                                     <  1:3277843/76‑1 (MQ=255)

                   gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc                                                                       >  2:5798458/1‑58 (MQ=255)

                    ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAg                                                    >  2:4518377/1‑76 (MQ=255)

                     cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGc                                                   <  1:1855431/76‑1 (MQ=255)

                     cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGc                                                   >  2:1506203/1‑76 (MQ=255)

                      gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGca                                                  >  1:1791476/1‑76 (MQ=255)

                      gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGca                                                  <  2:3727301/76‑1 (MQ=255)

                      gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGca                                                  >  1:2611942/1‑76 (MQ=255)

                      gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGca                                                  >  1:1438335/1‑76 (MQ=255)

                       cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGcac                                                 >  1:3940642/1‑76 (MQ=255)

                       cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGcac                                                 <  2:2096368/76‑1 (MQ=255)

                        gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGcaca                                                <  1:4156736/76‑1 (MQ=255)

                        gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGcaca                                                <  1:1470113/76‑1 (MQ=255)

                         ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACaa                                               >  1:3304008/1‑76 (MQ=255)

                         ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACaa                                               >  1:644440/1‑76 (MQ=255)

                          cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAat                                              >  2:3618884/1‑76 (MQ=255)

                            gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAatat                                            <  2:2842427/76‑1 (MQ=255)

                            gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAatat                                            >  2:3247551/1‑76 (MQ=255)

                              ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGa                                          <  1:799954/76‑1 (MQ=255)

                              ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGa                                          <  2:2452794/76‑1 (MQ=255)

                              ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGa                                          >  1:2897881/1‑76 (MQ=255)

                              ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGa                                          >  1:4694194/1‑76 (MQ=255)

                                acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAAt                                        >  2:2753527/1‑76 (MQ=255)

                                 cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATc                                       <  1:5458757/76‑1 (MQ=255)

                                       cgcttcgaACGCAAGCACTGGGGGGGAAAGTAGCG‑TCCC‑CACCCCCAGGCCACGCAGCACAATATGAATCGCCAGt                                 <  2:150612/74‑1 (MQ=255)

                                       cacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                 >  1:2609099/1‑76 (MQ=255)

                                          ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                              <  1:5675142/76‑1 (MQ=255)

                                          ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                              <  1:2659472/76‑1 (MQ=255)

                                          ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                              <  1:4233461/76‑1 (MQ=255)

                                          ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                              <  2:4478068/76‑1 (MQ=255)

                                           tcaaACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCGCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                             <  1:5327778/76‑1 (MQ=255)

                                            caaACGCAATCACTGAGCGGGGAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCCGTCCggc                            >  2:1204274/1‑75 (MQ=255)

                                            caaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                            <  2:3390787/76‑1 (MQ=255)

                                            caaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                            >  2:5001835/1‑76 (MQ=255)

                                            caaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                            >  2:433333/1‑76 (MQ=255)

                                            caaACGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGCCCGGt                            >  2:702679/1‑76 (MQ=39)

                                             aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAGCCCCACGCCAGCCAGCACAATATGAATCGCCAGTCCGGTc                           <  2:1523131/76‑1 (MQ=255)

                                                cGCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGGCTGa                        >  2:3356767/1‑76 (MQ=255)

                                                 gCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                       <  1:3935403/76‑1 (MQ=255)

                                                 gCAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                       >  1:2691663/1‑76 (MQ=255)

                                                  cAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                      >  1:3702751/1‑76 (MQ=255)

                                                  cAATCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                      >  1:872959/1‑76 (MQ=255)

                                                   aaTCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCACCCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     >  2:4177237/1‑76 (MQ=39)

                                                   aaTCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     >  1:1037148/1‑76 (MQ=255)

                                                   aaTCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     >  2:4144738/1‑76 (MQ=255)

                                                   aaTCACGGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     <  2:1199335/76‑1 (MQ=255)

                                                     tCACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                   >  2:3474216/1‑76 (MQ=255)

                                                      cACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcccc                  >  2:2817467/1‑75 (MQ=255)

                                                      cACTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                  >  1:171555/1‑76 (MQ=255)

                                                       aCTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                 >  2:3274782/1‑76 (MQ=255)

                                                        cggAGCGGGAAAATAGCG‑TCCC‑CAGCCCCACGCCAGCCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                <  2:5072854/74‑1 (MQ=255)

                                                        cTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                >  1:881642/1‑76 (MQ=255)

                                                        cTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                >  1:3250307/1‑76 (MQ=255)

                                                        cTGAGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                >  1:2781612/1‑76 (MQ=255)

                                                         tGAGCGGGAAAATAGCG‑TACC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa               >  1:1582469/1‑76 (MQ=255)

                                                           aGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt             >  1:4086310/1‑76 (MQ=255)

                                                           aGCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt             <  1:331006/76‑1 (MQ=255)

                                                           aGCGGAAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt             <  1:5820067/76‑1 (MQ=255)

                                                            gCGGGAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTCAAGCCCAGCACTa            >  2:1580404/1‑76 (MQ=255)

                                                             cGGGAAAATAGCG‑TCCC‑CACCCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg           <  2:772199/76‑1 (MQ=255)

                                                              gggAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          >  2:4828733/1‑76 (MQ=255)

                                                              gggAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          <  1:989802/76‑1 (MQ=255)

                                                              gggAAAAGAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          <  2:1832035/76‑1 (MQ=255)

                                                               ggAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTTAAGCCCAGCAATAGAc         >  2:1822265/1‑76 (MQ=255)

                                                               ggAAAATAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc         <  2:5378014/76‑1 (MQ=255)

                                                                 aaaaTAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc       <  2:3728705/76‑1 (MQ=255)

                                                                 aaaaTAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc       <  1:1576588/76‑1 (MQ=255)

                                                                  aaaTAGCGCC‑CC‑CATCCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc      <  2:4438997/76‑1 (MQ=255)

                                                                   aaTAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa     <  2:3165182/76‑1 (MQ=255)

                                                                    aTAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc    <  1:4709901/76‑1 (MQ=255)

                                                                    aTAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc    >  1:1391998/1‑76 (MQ=255)

                                                                     tAGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGACAGTCCGGTCTGAAGCCCAGCAATAGACGc       >  1:574163/1‑72 (MQ=255)

                                                                     aaGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc   <  1:4238632/75‑1 (MQ=255)

                                                                      aGCG‑TCCC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg  >  1:3391041/1‑76 (MQ=255)

TCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CC‑CAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCG  >  NC_000913/4006406‑4006550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`