BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,805	T→C	W222R (TGG→CGG)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,805	0	T	C	98.6%	445.6 / ‑4.7	139	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (58/79); total (59/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCA  >  NC_000913/3822715‑3822893

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                      <  1:1105408/93‑1 (MQ=255)

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                      >  1:64383/1‑93 (MQ=255)

   ccGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGct                                                                                     <  2:587141/93‑1 (MQ=255)

     ggTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctac                                                                                   <  1:3381027/93‑1 (MQ=255)

     ggTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctac                                                                                   <  1:135019/93‑1 (MQ=255)

     ggTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctac                                                                                   <  1:739194/93‑1 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                  >  1:3684405/1‑93 (MQ=255)

       tGCTGGCGAAAGTCGCAAACCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                 <  1:1815126/93‑1 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                 <  1:865595/93‑1 (MQ=255)

       tGCTGACGAAATTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                 <  1:1815124/93‑1 (MQ=255)

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCa                                                                             >  2:4655044/1‑93 (MQ=255)

            gCGAAAGTCTCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                            <  1:2719529/93‑1 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGccg                                                                            >  2:401179/1‑91 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                            <  1:4567520/93‑1 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                            <  2:4671985/93‑1 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                            <  1:2264960/93‑1 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                          <  2:1576257/93‑1 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGGTACTATGCAGGCTg                                                                        >  2:4637090/1‑93 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTg                                                                        >  1:636779/1‑93 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTg                                                                        >  2:1144440/1‑93 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACCATGCAGGCTg                                                                        >  1:1249000/1‑93 (MQ=255)

                 aGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGc                                                                       <  1:4749721/93‑1 (MQ=255)

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCa                                                                      <  1:1964648/93‑1 (MQ=255)

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCa                                                                      <  2:2195310/93‑1 (MQ=255)

                   tCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAg                                                                     <  1:4395514/93‑1 (MQ=255)

                     gCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAggg                                                                   <  1:1387964/93‑1 (MQ=255)

                      cAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAggga                                                                  >  1:3321880/1‑92 (MQ=255)

                      cAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGt                                                                  <  1:2061909/93‑1 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTaa                                                                <  2:1237361/93‑1 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTCCAGGGTaa                                                                >  2:3299663/1‑93 (MQ=255)

                         aaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAg                                                               <  2:2230580/93‑1 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGt                                                              <  2:2766903/93‑1 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCa                                                          <  1:2834508/93‑1 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGCCa                                                          >  2:4313950/1‑93 (MQ=255)

                                  gcctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaa                                                      <  2:3131213/93‑1 (MQ=255)

                                  gcctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaa                                                      <  2:3357467/93‑1 (MQ=255)

                                   cctgcatTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaaa                                                     <  2:2291843/93‑1 (MQ=255)

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAg                                                    <  2:3931643/93‑1 (MQ=255)

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAg                                                    >  1:2402847/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAgg                                                   >  1:2026246/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGAAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGGCATGAAAgg                                                   >  1:4091812/1‑93 (MQ=255)

                                      tcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                  <  1:1687528/92‑1 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                  <  2:4557832/93‑1 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                  <  1:64339/93‑1 (MQ=255)

                                       cctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                 <  2:439974/93‑1 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCCTGAAAGGGCGcc                                              >  2:3634138/1‑93 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                              <  1:3629237/93‑1 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATACCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGGCGcc                                              >  1:3404980/1‑93 (MQ=255)

                                           ttCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCt                                             <  2:1998363/92‑1 (MQ=255)

                                           aTCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCt                                             <  2:2969762/93‑1 (MQ=255)

                                           aTCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGCCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCt                                             <  2:2774693/93‑1 (MQ=255)

                                            tCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATTGAGGCTGCAGGGGAAGTGCCATGAAAGGGCGCCTg                                            >  2:255426/1‑93 (MQ=255)

                                            tCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGCCATGAAAGGTCGCCTg                                            >  1:3491978/1‑93 (MQ=255)

                                               ccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTACGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTa                                         <  1:2260013/93‑1 (MQ=255)

                                                cACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCCACTATGCAGGCTGCAGGGTAAGAGCCATGAAAGGTCGCCTGTTAg                                        >  2:1397386/1‑93 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                       >  2:2477622/1‑93 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                       <  2:774243/93‑1 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                       <  2:1925589/93‑1 (MQ=255)

                                                  cGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGAt                                      <  1:4666891/93‑1 (MQ=255)

                                                  cGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGAt                                      <  2:2394330/93‑1 (MQ=255)

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTGCTATGCAGGCTGCAGGCGAAGAGCCAGGAAAGGTCGCTTGGTAGAttc                                    >  1:4031245/1‑91 (MQ=255)

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGc                                    <  1:3012376/93‑1 (MQ=255)

                                                     cAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCt                                   <  2:382895/93‑1 (MQ=255)

                                                      aaTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGGTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTg                                  >  1:3231209/1‑93 (MQ=255)

                                                       aTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGt                                 >  1:3145378/1‑93 (MQ=255)

                                                        tCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTc                                >  2:3168033/1‑93 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTccc                              >  2:1459614/1‑93 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTccc                              >  1:3054130/1‑93 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                             <  2:3899727/93‑1 (MQ=255)

                                                           gCAAGGAGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                             >  2:233113/1‑93 (MQ=255)

                                                            cAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCAc                            <  1:673660/93‑1 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                          >  1:3495706/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                          <  1:3896082/93‑1 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCGATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                         >  2:2180526/1‑93 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                         <  1:2043432/93‑1 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                         >  1:2227746/1‑93 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGtt                      <  2:4136039/93‑1 (MQ=255)

                                                                   agaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTc                     <  1:1291328/93‑1 (MQ=255)

                                                                    gaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTcc                    <  1:4080273/93‑1 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                   <  2:1241997/93‑1 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                   <  2:3250782/93‑1 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                   <  2:3717830/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCGGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                  <  1:1774323/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                  <  1:1774324/93‑1 (MQ=255)

                                                                        gAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTaa                <  1:722844/93‑1 (MQ=255)

                                                                         aaGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAAc               >  1:1019355/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg              >  1:1129861/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg              <  2:3492142/93‑1 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg             <  2:4053464/93‑1 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg             >  1:4298168/1‑93 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGGTCCCTAACgg             >  1:22045/1‑93 (MQ=255)

                                                                            gtGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACggg            <  1:1958320/92‑1 (MQ=255)

                                                                             tGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgggg           >  1:1576870/1‑92 (MQ=255)

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCCCCCCGTTCCCTAACGGGCg          >  1:170475/1‑93 (MQ=255)

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGt         <  1:1509233/93‑1 (MQ=255)

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGAAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTCACGGGCGtg        >  1:471653/1‑92 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg       >  2:2677816/1‑93 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg       <  1:280510/93‑1 (MQ=255)

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgc     >  2:2264743/1‑93 (MQ=255)

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca  <  2:2241136/93‑1 (MQ=255)

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCA  >  NC_000913/3822715‑3822893

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