BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,907,388	T→C	E213E (GAA→GAG)	pstC ←	PstC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,907,388	0	T	C	99.7%	1147.5 / ‑6.1	363	E213E (GAA→GAG)	pstC	PstC
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (175/186); total (177/186)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 375 total aligned reads displayed.

CATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTA  >  NC_000913/3907298‑3907477

 aTGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTcc                                                                                        <  2:5329398/93‑1 (MQ=255)

     tttCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAgg                                                                                    <  1:2877602/93‑1 (MQ=255)

          ccGACCGCGCGCCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGgtggtg                                                                               <  2:1008482/93‑1 (MQ=255)

             aGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAg                                                                            >  1:2227983/1‑93 (MQ=255)

                gcgcgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCa                                                                         <  1:6105615/93‑1 (MQ=255)

                 cgcgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCaa                                                                        >  1:6831650/1‑93 (MQ=255)

                  gcgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAAt                                                                       <  1:5289853/93‑1 (MQ=255)

                  gcgGCCAAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAAt                                                                       <  2:7122236/93‑1 (MQ=255)

                   cgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATa                                                                      <  2:4043983/93‑1 (MQ=255)

                   cgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATa                                                                      >  1:1970257/1‑93 (MQ=255)

                     gCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATAcc                                                                    <  1:5939627/93‑1 (MQ=255)

                      cccagccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCg                                                                   <  2:2698231/93‑1 (MQ=255)

                        cagccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  1:7907708/93‑1 (MQ=255)

                        cagccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  1:4179877/93‑1 (MQ=255)

                        cagccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  2:7343090/93‑1 (MQ=255)

                        cagccccagcATGATGCCGACGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  1:5069702/93‑1 (MQ=255)

                        cagccccagaATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  1:7306161/93‑1 (MQ=255)

                         agccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAg                                                                <  2:3035046/93‑1 (MQ=255)

                         agccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAg                                                                <  1:641621/93‑1 (MQ=255)

                          gccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAgg                                                               <  1:6785351/93‑1 (MQ=255)

                           ccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGc                                                              <  1:483200/93‑1 (MQ=255)

                            cccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGcc                                                             >  1:6926445/1‑93 (MQ=255)

                             ccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCg                                                            >  2:6809174/1‑93 (MQ=255)

                             ccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCg                                                            <  2:4975364/93‑1 (MQ=255)

                              cagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCGCCCAGGTGGTGCAGCCAATACCGTAGGCCGa                                                           >  2:7699059/1‑93 (MQ=255)

                               agcATGATGCCGCCGATAACACGATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGAc                                                          >  2:5736054/1‑93 (MQ=255)

                               agcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGAc                                                          >  2:6617311/1‑93 (MQ=255)

                                gcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGAct                                                         <  1:135346/93‑1 (MQ=255)

                                 cATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActc                                                        <  2:2881800/93‑1 (MQ=255)

                                 cATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACATCCCAGGTGGTGCAGCCAATACCGTAGGCCGActc                                                        >  2:4322463/1‑93 (MQ=255)

                                  atgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActct                                                       <  2:5188861/93‑1 (MQ=255)

                                   tgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCtt                                                      >  2:6414184/1‑93 (MQ=255)

                                   tgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCtt                                                      <  2:6093768/93‑1 (MQ=255)

                                    gatgCCGCCGATAACCCCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCttt                                                     <  1:2376424/93‑1 (MQ=255)

                                    gatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCttt                                                     >  1:4962487/1‑93 (MQ=255)

                                     atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtc                                                    <  2:7818269/93‑1 (MQ=255)

                                     atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtc                                                    <  1:2668074/93‑1 (MQ=255)

                                     atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtc                                                    <  1:198831/93‑1 (MQ=255)

                                      tgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtca                                                   >  1:2842559/1‑93 (MQ=255)

                                       gccgccGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcat                                                  >  2:3149481/1‑93 (MQ=255)

                                       gccgccGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcat                                                  >  1:7736801/1‑93 (MQ=255)

                                          gccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatcat                                               <  2:6858897/93‑1 (MQ=255)

                                            cgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatcatca                                             >  1:2129829/1‑93 (MQ=255)

                                             gATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCAc                                            >  2:5619959/1‑93 (MQ=255)

                                              aTAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCAcc                                           <  2:758993/93‑1 (MQ=255)

                                               tAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCg                                          >  2:4801252/1‑93 (MQ=255)

                                                aaCACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCgg                                         <  2:2022857/93‑1 (MQ=255)

                                                 acacCATTTTTGTTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  1:1261889/93‑1 (MQ=255)

                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        >  2:6246484/1‑93 (MQ=255)

                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  1:1094802/93‑1 (MQ=255)

                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  1:317436/93‑1 (MQ=255)

                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        >  1:5186675/1‑93 (MQ=255)

                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  1:4677558/93‑1 (MQ=255)

                                                  cacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCgggg                                       >  2:370769/1‑93 (MQ=255)

                                                  cacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCgggg                                       <  1:4785708/93‑1 (MQ=255)

                                                   acCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGt                                      <  2:2179742/93‑1 (MQ=255)

                                                   acCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGt                                      >  2:1438005/1‑93 (MQ=255)

                                                     cATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGttt                                    <  1:1189874/93‑1 (MQ=255)

                                                     cATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGttt                                    >  1:6569463/1‑93 (MQ=255)

                                                      aTTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTg                                   <  1:3357484/93‑1 (MQ=255)

                                                       tttttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGt                                  >  2:6565689/1‑93 (MQ=255)

                                                       tttttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGt                                  <  1:1447107/93‑1 (MQ=255)

                                                       tttttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGt                                  >  2:5840766/1‑93 (MQ=255)

                                                        ttttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGtt                                 <  1:4457746/93‑1 (MQ=255)

                                                          ttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCg                               <  1:1011289/93‑1 (MQ=255)

                                                             gTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAac                            <  2:5932499/93‑1 (MQ=255)

                                                               gAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAacaa                          >  1:2026172/1‑92 (MQ=255)

                                                                aaCGGAAGAACGATACGCCAGATAGCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAAcaca                         <  2:3270873/93‑1 (MQ=255)

                                                                 aCGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACAt                        <  1:3233005/93‑1 (MQ=255)

                                                                 aCGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACAt                        <  1:5475132/93‑1 (MQ=255)

                                                                  cGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATc                       <  2:1125004/93‑1 (MQ=255)

                                                                  cGGAAGAACGATACGCCAGATAACCCCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATc                       <  1:1777903/93‑1 (MQ=255)

                                                                   ggAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCa                      <  2:3509549/93‑1 (MQ=255)

                                                                   ggAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCa                      <  2:2377735/93‑1 (MQ=255)

                                                                    gaagaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCAc                     >  2:1786990/1‑93 (MQ=255)

                                                                    gaagaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCAc                     <  1:3998810/93‑1 (MQ=255)

                                                                     aagaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACg                    <  1:349150/93‑1 (MQ=255)

                                                                       gaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATATCGTATGCCGACTTTTTCTTCATCACAGGGGTTTGTTCGAACACATCACGCa                  >  1:5724708/1‑93 (MQ=255)

                                                                       gaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCa                  <  2:2227982/93‑1 (MQ=255)

                                                                       gaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCa                  <  2:5367340/93‑1 (MQ=255)

                                                                        aaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCAt                 <  1:7388018/93‑1 (MQ=255)

                                                                        aaCGATACGCCAGATAACCTCCCAGGTGGTGCAACCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCAt                 <  1:7625336/93‑1 (MQ=255)

                                                                         aCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCAtt                <  1:84207/93‑1 (MQ=255)

                                                                         aCGATACGCCAGATAACCTCCCAGGTGGGGGAGCCAATACCGCAGGGCGACTTTTTAATCTTCACCCGGGTTTGTTTGTCCACATCACGCact                >  2:1597324/1‑91 (MQ=255)

                                                                          cGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCTACTCTTGCATCATCACCGGGGTTTGTTCGAACACATCACGCAttt               <  2:3927721/93‑2 (MQ=255)

                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  1:5603646/93‑1 (MQ=255)

                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  2:2502467/93‑1 (MQ=255)

                                                                                 ccAGATAACCTCCCAGGTGGTGCAGCCAATTCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTCCcgccgc        <  2:4765769/93‑1 (MQ=255)

                                                                                  cAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCa       >  1:1704950/1‑93 (MQ=255)

                                                                                   aGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCaa      >  1:4250225/1‑93 (MQ=255)

                                                                                   aGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCaa      >  1:2589318/1‑93 (MQ=255)

                                                                                    gATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTATCGCCGCAAt     >  2:2537776/1‑93 (MQ=255)

                                                                                    gATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAAt     <  2:458381/93‑1 (MQ=255)

                                                                                     aTAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATg    >  2:1322793/1‑93 (MQ=255)

                                                                                      tAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGt   <  2:1555393/93‑1 (MQ=255)

                                                                                      tAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGt   >  2:815592/1‑93 (MQ=255)

                                                                                       aaCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa  <  1:3749864/93‑1 (MQ=255)

                                                                                       aaCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa  <  2:5445696/93‑1 (MQ=255)

                                                                                       aaCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa  <  1:6222016/93‑1 (MQ=255)

                                                                                       aaCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa  >  1:7165408/1‑93 (MQ=255)

CATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTA  >  NC_000913/3907298‑3907477

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