BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,183,028	T→C	V70A (GTC→GCC)	b3042 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,183,028	0	T	C	99.3%	840.2 / ‑3.3	272	V70A (GTC→GCC)	b3042	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (116/154); total (117/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 287 total aligned reads displayed.

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117

 tgtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCt                                                                                        >  2:4688661/1‑93 (MQ=255)

 tgtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCt                                                                                        <  2:7648638/93‑1 (MQ=255)

  gtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                       >  2:2703272/1‑93 (MQ=255)

  ggggAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                       <  1:6300082/91‑1 (MQ=255)

     gagaAAAAAATCCGCCAAGCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                    <  1:5022242/93‑1 (MQ=255)

     gagaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                    <  2:3673316/93‑1 (MQ=255)

      agaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   <  1:6830254/93‑1 (MQ=255)

      agaAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   <  1:1764531/93‑1 (MQ=255)

       gaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  >  2:4249784/1‑93 (MQ=255)

       gaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  <  2:7799692/93‑1 (MQ=255)

        aaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACgg                                                                                 <  1:1170173/93‑1 (MQ=255)

         aaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                                <  2:5190968/93‑1 (MQ=255)

         aaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                                >  1:7890181/1‑93 (MQ=255)

         aaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                                >  1:4731195/1‑93 (MQ=255)

          aaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGAc                                                                               >  1:6000009/1‑93 (MQ=255)

          aaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGAc                                                                               >  1:3130281/1‑93 (MQ=255)

           aaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACt                                                                              >  2:2765762/1‑93 (MQ=255)

              aTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                           >  2:1848777/1‑93 (MQ=255)

               cccGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTg                                                                          <  2:6134911/92‑1 (MQ=255)

                ccGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGa                                                                         <  2:3356989/93‑1 (MQ=255)

                ccGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGa                                                                         >  1:9073/1‑93 (MQ=255)

                 cGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaa                                                                        <  2:7338421/93‑1 (MQ=255)

                 cGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaa                                                                        <  2:6797042/93‑1 (MQ=255)

                   ccAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaa                                                                      >  2:1192090/1‑93 (MQ=255)

                    cAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaaa                                                                     <  1:734526/93‑1 (MQ=255)

                    cAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaaa                                                                     <  2:6480130/93‑1 (MQ=255)

                    cAAACCCGACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaaa                                                                     <  2:167743/93‑1 (MQ=255)

                     aaaCCCTACAGGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAAc                                                                    <  1:2913649/93‑1 (MQ=255)

                     aaaCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAAc                                                                    <  1:1071267/93‑1 (MQ=255)

                       aCCCTACAAGCGCAGCTGACGTACCTCGACCTTGCAAGCCGTGAAGAGTTCGACGTCCACACGCAAGCCCTGTTACGCACTCGTGAAAAACTg                                                                  >  2:951032/1‑93 (MQ=255)

                       aCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                  >  1:2495387/1‑93 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 >  1:4247960/1‑93 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 <  2:6611511/93‑1 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                                >  1:4983862/1‑93 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                                >  2:488702/1‑93 (MQ=255)

                          cTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCGTGTTACGGACTCGTGAAAAACTGgcg                                                               >  1:2601987/1‑93 (MQ=255)

                          cTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcg                                                               >  2:56325/1‑93 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  1:744249/1‑93 (MQ=255)

                           caCAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  2:5806199/92‑1 (MQ=255)

                           aaCAAGCGCAGCGGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  1:1585324/92‑1 (MQ=255)

                                 cgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTg                                                        <  1:2755681/93‑1 (MQ=255)

                                 cgcAGCGGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTg                                                        <  2:507839/93‑1 (MQ=255)

                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                       >  1:4927164/1‑93 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  2:7231168/93‑1 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  2:4992519/93‑1 (MQ=255)

                                     gCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACa                                                    <  1:3816965/93‑1 (MQ=255)

                                      cTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa                                                   <  2:474503/93‑1 (MQ=255)

                                        gACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAg                                                 <  1:1267693/93‑1 (MQ=255)

                                         aCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa                                                <  2:687938/93‑1 (MQ=255)

                                          cgcgCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGaa                                               <  2:2336858/93‑1 (MQ=255)

                                          cgcgCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGaa                                               <  1:1671607/93‑1 (MQ=255)

                                           gcgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAAt                                              <  1:7668573/93‑1 (MQ=255)

                                           gccccTCGACCTTGTAAGCCGTAAAGACTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAAt                                              <  1:6408949/90‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  2:4500246/93‑1 (MQ=255)

                                               cTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGc                                          <  1:2434849/93‑1 (MQ=255)

                                                tCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCg                                         >  2:3511283/1‑93 (MQ=255)

                                                 cGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGa                                        >  1:5728604/1‑93 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      >  1:2810566/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  1:1846919/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  1:7297633/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  2:7470591/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGGGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     <  1:1474055/93‑1 (MQ=255)

                                                     cTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTg                                    >  1:6779321/1‑93 (MQ=255)

                                                     cTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTg                                    >  2:6054038/1‑93 (MQ=255)

                                                       tGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGa                                  >  1:2191549/1‑93 (MQ=255)

                                                         aaaGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                >  2:5864633/2‑93 (MQ=255)

                                                          aaGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaaa                               >  2:1313618/1‑93 (MQ=255)

                                                          aaGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaaa                               <  1:2258784/93‑1 (MQ=255)

                                                           aGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAAt                              >  2:7077025/1‑93 (MQ=255)

                                                           aGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAAt                              <  2:4453555/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  1:3631290/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  1:4096493/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:6642632/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  1:7158977/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:3916515/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:388209/93‑1 (MQ=255)

                                                                 gAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGt                        <  1:4019237/93‑1 (MQ=255)

                                                                 gAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGt                        <  1:5562612/93‑1 (MQ=255)

                                                                   agagTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAc                      >  1:6526685/1‑93 (MQ=255)

                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                   <  2:4650618/93‑1 (MQ=255)

                                                                       ttCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgag                  <  2:280383/93‑1 (MQ=255)

                                                                       ttCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgag                  <  2:598903/93‑1 (MQ=255)

                                                                       ttCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgag                  <  1:458615/93‑1 (MQ=255)

                                                                         cGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGAt                <  1:7775844/93‑1 (MQ=255)

                                                                            cGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaa             <  1:4139710/93‑1 (MQ=255)

                                                                              tCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaa           >  2:3061751/1‑93 (MQ=255)

                                                                              tCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaa           >  1:307018/1‑93 (MQ=255)

                                                                               ccAAACGCAAGCCGTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaaa          <  2:2868882/93‑1 (MQ=255)

                                                                               ccAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaaa          <  2:3157668/93‑1 (MQ=255)

                                                                               ccAAACGCAAGCCCTGTGACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaaa          <  2:893647/93‑1 (MQ=255)

                                                                                cAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAag         >  1:6153167/1‑93 (MQ=255)

                                                                                cAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAag         >  2:1991550/1‑93 (MQ=255)

                                                                                  aaCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagca       <  2:4035505/93‑1 (MQ=255)

                                                                                  aaCCCAAGCACTGTTACGGAATCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagca       <  1:2092922/93‑1 (MQ=255)

                                                                                   aCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcag      <  1:1102693/93‑1 (MQ=255)

                                                                                   aCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcag      <  1:7015644/93‑1 (MQ=255)

                                                                                   aCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTCAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcag      >  2:5691630/1‑93 (MQ=255)

                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAACAAGCAGCCa   >  2:439319/1‑93 (MQ=255)

                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   >  2:5336167/1‑93 (MQ=255)

                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   <  2:2549320/93‑1 (MQ=255)

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117

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